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Preface 

Beginning with Volume XX, the Deep Spi.ee Network Progress Report changed 
Irom the Technical Report 32- series to the Progress Report 42- series. The volume 
number continues the sequence of the preceding issues. Thus. Progress Report 
12-20 is llit* twentieth volume of tin* Dee p Space Network series, and is an uninter- 
rupted follow-on to Technical Report 32-1526. Volume XIX. 

This report presents DSN progress in flight project support, tracking and data 
acquisition (TDA) research and technology, network engineering, hardware and 
software implementation, and operations. Each issue presents material in some, 
but not all. of the following categories in the order indicated. 

Description of the DSN 

Mission Support 

Ongoing Planetary/ Interplanetary Flight Projects 
Advanced Flight Projects 

Radio Science 

Special Projects 

Supporting Research and Technology 
Tracking and Ground-Based Navigation 
Communications — Spacecraft/Ground 
Station Control and Operations Technology' 

Network Control and Data Processing 

Network and Facility Engineering and Implementation 
Network 

Network Operations Control Center 
Ground Communications 
Deep Space Stations 

Operations 

Network Operations 
Network Operations Control Center 
G r oun d Com m un i ca t i o n s 
Deep Space Stations 

Program Planning 
TDA Planning 
Quality' Assurance 

In each issue, the part entitled “Description of the DSN” describes the functions 
and facilities o‘. the DSN and may report the current configuration of one of the 
five DSN systems (Tracking, Telemetry, Command, Monitor & Control, and Test 
6: Training). 

The work described in this report series is either performed or managed by’ the 
Tracking and Data Acquisition organization of JPL for NASA. 
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Network Functions and Facilities 

N. A. Rer^etti 

Office of Tracking and Data Acquisition 


The objectives, functions, and organization of the Deep Space Network are 
summarized; deep space station, ground communication , and network 
operations control capabilities are described. 


The Deep Space Network (DSN), established by the 
National Aeronautics and Space Administration (NASA) 
Office of Tracking and Data Acquisition under the system 
management and technical direction of the Jet Propulsion 
Laboratory (JPL), is designed for two-way communications 
with unmanned spacecraft traveling approximately 16,000 
km (10,000 miles) from Earth to the farthest planets of our 
solar system. It has provided tracking and data acquisition 
support for the following NASA deep space exploration 
projects: Banger, Surveyor, Mariner Venus 1962, Mariner 
Mars 1964, Mariner Venus 1967, Mariner Mars 1969, 
Mariner Mars 1971, and Mariner Venus Mercury 1973, for 
which JPL has been responsible for the project manage- 
ment, the development of the spacecraft, and the conduct 
of mission operations; Lunar Orbiter, for which the 
Langley Research Center carried out the project manage- 
ment, spacecraft development, and conduct of mission 
operations; Pioneer, for which Ames Pesearch Center 
carried out the project management, spacecraft develop- 
ment, and conduct of mission operations; and Apollo, for 
which the Lyndon B. Johnson Space Center was the 
project center and the Deep Space Network supple- 
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mented the Manned Space flight Network (MSFN), 
which was managed by the Goddard Space Flight Center 
(GSFC). It is providing tracking and data acquisition 
support for Helios, a joint U.S./West German project; and 
Viking, for which Langley Research Center provides the 
project management, the Lander spacecraft, and conducts 
mission operations, and for which JPL also provides the 
Orbiter spacecraft.. 

The Deep Space Network is one of two NASA 
networks. The other, the Spaceflight Tracking and Data 
Network, is under the system management and technical 
direction of the Goddard Space Flight Center. Its function 
is to support manned and unmanned Earth-orbiting 
satellites. The Deep Space Network supports lunar, 
planetary, and interplanetary flight projects. 

From its inception, NASA has had the objective of 
conducting scientific investigations throughout the solar 
system. It was recognized that in order to meet this 
objective, significant supporting research and advanced 
technology development must be conducted in order to 
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provide deep space telecommunications for science data 
return in a cost effective manner. Therefore, the Network 
is continually evolved to keep pace with the state of the 
art of telecommunications and data handling. It was also 
recognized early that close coordination would he needed 
between the requirements of the flight projects for data 
return and the capabilities needed in the Network. This 
close collaboration was effected by the appointment of a 
I racking and Data Systems Manager as part of the flight 
project team from the initiation of the project to the end 
ol the mission. By this process, requirements were 
identified early enough to provide funding and implemen- 
tation in time for use by the flight project in its flight 
phase. 

As of July 19/2, NASA undertook a change in the 
interface between the Network and the flight projects. 
Prior to that time, since 1 January 1964, in addition to 
consisting of the Deep Space Stations and the Ground 
Communications Facility, the Network had also included 
the mission control and computing facilities and provided 
the equipment in the mission support areas for the 
conduct of mission operations. The latter facilities were 
housed in a building at [PL known as the Space Flight 
Operations Facility (SFOF), The interface change was'to 
accommodate a hardware interface between the support 
of the network operations control functions and those of 
the mission control and computing functions. This resulted 
in the flight projects assuming the cognizance of the large 
general-purpose digital computers which were used for 
both network processing and mission data processing. 
They also assumed cognizance of all of the equipment in 
the flight operations facility for display and communica- 
tions necessary i or the conduct of mission operations. The 
Network then undertook the development of hardware 
and computer software necessary to do its network 
operations control and monitor functions in separate 
computers. This activity has been known as the Network 
Control System Implementation Project. A characteristic 
of the new interface is that the Network provides direct 
data flow to and from the stations; namely, metric data, 
science and engineering telemetry, and such network- 
monitor data as are useful to the flight project. This is 
done via appropriate ground communication equipment 
to mission operations centers, wherever they may be. 

4 he principal deliverables to the users of the Network- 
are carried out by data system configurations as follows; 

• The DSN Tracking System generates radio metric 
data; i.e., angles, one- and two-way doppler and 
range, and transmits raw data to Mission Control. 


• I he DSN Iclemetry System receives, decodes, 
records, and retransmits engineering and scientific 
data generated in the spacecraft to Mission Control. 

• I he DSN Command System accepts coded signals 
from Mission Control via the Ground Communica- 
tions facility and transmits them to the spacecraft in 
order to initiate spacecraft functions in flight. 

I he data system configurations supporting testing, 

training, and network operations control functions are as 

follows: 

• I he DSN Monitor and Control System instruments, 
transmits, records, and displays those parameters of 
the D$_ J necessary to verify configuration and 
validate the Network. It provides operational 
direction and configuration control of the Network, 
and provides primary interface with flight project 
Mission Control personnel. 

• The DSN Test and Training System generates and 
controls siinulared data to support development, 
test, training and fault isolation within the DSN. It 
participates mission simulation with flight pro- 
jects. 

The capabilities needed to carry out the above functions 
have evolved in three technical areas: 

(1) The Deep Space Stations, which are distributed 
around Earth and which, prior to 1964, formed part 
of the Deep Space Instrumentation Facility. The 
technology involved in equipping these stations is 
strongly related to the state of the art of telecommu- 
nications and flight-ground design considerations, 
and is almost completely multi mission in character. 

(2) 'The Ground Communications Facility provides the 
capability required for the transmission, reception, 
and monitoring of Earth-based, point-to-point com- 
munications between the stations and the Network 
Operations Control Center at JPL, Pasadena, and to 
the inis; ion operations centers, wherever they may 
be. Four communications disciplines arc provided: 
teletype, voice, high-speed, and wideband. 'The 
Ground Communications Facility uses the capabili- 
ties provided by common carriers throughout the 
world, engineered into an integrated system by 
Goddard Space Flight Center, and controlled from 
the communications Center located in the Space 
Flight Operations Facility {Building 230) at JPL. 
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Cl) The Network Operations Control Center is tin* 
functional entity for centralized operational control 
of the Network and interfaces with the users, It has 
two separable functional elements; namely, Network 
Operations Control and Network Data Processing. 
The functions of the Network Operations Control 
are; 

• Control and coordination of Network support to 
meet commitments to Network users. 

• Utilization of the Network data processing 
computing capability to generate all standards 
and limits required for Network operations. 

• Utilization of Network data processing comput- 
ing capability to analyze and validate the 
performance of all Network systems. 

The personnel who carry out the above functions arc- 
located in the Space Plight Operations Facility, 
where mission operations functions are carried out 
by certain flight projects. Network personnel are 
directed by an Operations Control Chief. 


The functions of the Network Data Processing arc: 

• Processing of data used In Network Operations 
Control for control and analysis of the Network. 

• Display in the- Network Operations Control Area 
of data processed in the Network Data Process 
ing Area. 

• Interface with communications circuits for input 
to and output from the Network Data Processing 
Area. 

• Data logging and production of the intermediate 
data records, 

The personnel who carry out these functions are 
located approximately — <K) meters from the Space 
Flight Operations Facility. The equipment consists 
of mic.icomputcrs for real-time data system monitor- 
ing, two XDS Sigma 5s, display, magnetic tape 
recorders, and appropriate interface equipment with 
the ground data communications. 
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DSN Ground Communications Facility 

M. $. Glenn 

TDA Engineering 

The Ground Communications iaeilih) has hern designed to prot ide icliuhlc 
Earth- based, point-to-point voire mid data communications as part of the IX S.\ 


Tracking and Data Acquisition System. 


I. Introduction 

The Ground Coin in indent ions Facility (CCP ) is one of 
the three functional elements of the Deep Space Network 
and provides the capability required for the transmission, 
reception, and monitoring of voice and data communica- 
tions between the various locations of the DSN. The CCF 
uses common carrier circuits engineered by the NASA 
Communications (NASCOM) Division of the Goddard 
Space Flight Center, interconnected to specialized 
switching, terminal, and monitor equipment, integrated 
into a world-wide system, and operated in support of 
space flight missions. The CCF is composed of five 
subsystems: Teletype, Voice, High-Speed, Wideband, and 
Monitor. The Office of Tracking and Data Acquisition of 
the Jet Propulsion Laboratory provides the technical 
direction and systems management of the CCF and acts as 
the representative of NASCOM for communications 
switching functions on the west coast. 


II. GCF-NASCOM Interrelationships 

The interrelationships at the programmatic level 
between the Jet Propulsion Laboratory-developed Ground 


Communications Facility and the Coddard Space Might 
Center NASCOM are characterized as follows: 

(1) NASCOM 

(a) Provides ground communications ol all opera 
tional data for all NASA missions, including those 
supported by the DSN, 

(b) Accepts and supports communications require- 
ments established by the DSN. 

vc) Establishes basic characteristics of NASA opera- 
tional ground communications subsystems, such 
as teletype line rate, high-speed block size, and 
wideband circuit capability. 

(2) CCF 

(a) Provides ground communications lor all DSN 
missions using the services of NASCOM. 

(h) Establishes additional characteristics of all CCF 
subsystems, such as block multiplexing and error 
correction. 


I 

V. 
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III. Objectives and Goals 

The GCr design is based on standardized communica- 
tion techniques to provide more efficient transmission of 
user data and simple user and NASCOM interfaces. These 
objectives are met by: 

( 1 ) Providing message switching and routing. 

(2) Transmitting data which are essentially transparent, 
U\, user data are accepted and delivered in 
estab' shed formats and without additional errors. 

(3) Minimizing project-dependent equipment within the 

ccf. 

(4) Providing centralized data records. 

The continuing goals of the GCF are to provide highly 
reliable and cost-effective data transmission while main- 
taining a capability balance between the DSN and users, 
and include: 

(1) Equipment and routing redundancy to minimize 
single-point-of-failure impact. 

(2) Error performance which does not degrade data 
beyond KF-link error performance. 

(.3) Design coordinated with NASCOM Development 
Program. 

IV. Configuration and Functional Subsystems 

The current GCF configuration, including the related 
NASCOM functions, is shown in Fig. 1. This configuration 
is functionally organized into High-Speed Data, Wideband 
Data, Voice, Teletype, and Monitor Subsystems. 

A. High-Speed Data Subsystem 

The High-Speed Data (USD) Subsystem consists o ( 
assemblies that switch, transmit, receive, record, process, 
distribute, test, and monitor digital data and is used for the 
transmission of: 

( 1 ) All digital data of the DSN command, tracking, and 
monitor and control systems. 

{2> All low- or medium-rate data of the DSN telemetry 
system and the DSN test and training system. 

(3) All GCF monitor and control data. 

The High-Speed Data Subsystem provides a capability 
for transmitting and receiving serial bit stream formatted 
data over four-wire circuits having a 3.0-kl!z bandwidth. 
This serial bit stream is impressed on communication 
circuits at a continuous line bit rate divided into message 


segments referred to as high-speed blocks. The two key 
characteristics are: 

(1) Data blocks containing user data hits to be 
transmitted. 

(2) Filler blocks containing filler bits provided ny GCl- 
equipment when the user data bit rate is insufficient 
to maintain the fixed line bit rate required by design 
specifications. 

Each block is divided into three parts, header, text, and 
ending. Formats for the data blocks and filler blocks are 
illustrated in Tallies 1 and 2. 

The current plans are to provide the functional 
capabilities illustrated in Fig. 2 and to standardize at a 
1200-hit message segment and a line bit rate of 72(X) b/s. 
Other planned changes include conversion hem a 33-bit to 
a 22-bit error polynomi <1 and increasing the number of 
hits reserved in the data block ending from 36 to 40 bits to 
provide error correction, by re-transmission, for short 
outages or errors in GCF data transmission. 3 he purpose 
of these changes is to significantly reduce the error for 
nonreal-time replay of data. Figure 3 illustrates USD 
Subsystem configuration which is planned for the CY 1977 
and CY 1978 time period. The two configurations are 
required to provide continuous project support during the 
period of conversion from the existing Ground Data 
System to the new one for support of the Mariner 
Jupiter-Saturn and Pioneer Venus Projects. 

B. Wideband Data Subsystem 

The Wideband Data Subsystem consists of assemblies 
that switch, transmit, receive, process, distribute, test, and 
monitor data requiring the use of bandwidths greater than 
those provided by standard high-speed channels. The 
Wideband Data Subsystem illustrated in Fig. *1 includes 
standard wideband circuits as well as inter, site-microwave 
capabilities. The Wideband Subsystem is used for: 

(1) Telemetry data. 

(2) Simulation data. 

(3) Test and training data. 

(4) Data interchange within the DSN including opera- 
tional control. 

(5) Intrasite communications and timing signals. 

'Hie wideband data circuits to the deep space stations 
contain serial hit streams impressed on communication 
circuits at a continuous line bit rate typically 27.6, 28.5, 
50, 168, or 230.4 kilobits per second divided into 24(K)- or 
4800-bit message segments. Similar to the high-speed data. 
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tho message segments or data blocks contain user bits to 
be transmitted and filler bits provided by the CCF 
e<|uipment when the user data bit rate is insufficient to 
maintain the fixed line bit rate required by design 
specifications. I he data blocks are also divided into three 
parts: head, text, and ending, as illustrated in Tables 3 and 
4. 

C. Voice Subsystem 

The Voice Subsystem consists of assemblies that switch, 
transmit, receive, distribute, test, and monitor transmis- 
sions originally generated in vocal form and includes 
capabilities between the facilities of the Deep Space 
Network and to the Mission Control Centers. The 
functional capabilities and key characteristics include: 

(1) Standard voice-data grade circuits for all traffic. 

(2) Conferencing capability on one intercontinental 
circuit during noncritical periods for all deep space 
stations supporting a single project; individual 
circuits for each DSS during critical periods, 
resou rces perm i 1 1 ing. 

(3) User-controlled intercom switching. 

(4) Circuits used for high-speed data transmission 
(backup) if required. 

(5) Voice traffic recording in the central communica- 
tions terminal upon request. 

D. Teletype Subsystem 

I his subsystem consists of assemblies that switch, 
traismit, receive, distribute, test, and monitor digital 


signals originally generated in Baudot format at a teletype 
(I I\) rate of 100 words per minute. The operational use 
of teletype continues to be de-ernphasiz.ed and is used 
primarily for emergency, backup operational transmissions 
and administrative communications. Services and key 
characteristics include: 

(1) Handling Air Force Eastern Test Range (AFETR)- 
generated predicts for DSN initial acquisiton. 

(2) I ransmitting non-operational messages between the 
JPL Message Center and other locations. 

(3) Use of standard NASCOM format and the NASCOM 
communications processor for message switching. 

E. Monitor and Control Subsystem 

The Monitor and Control Subsystem consists of 
assemblies that gather, calculate, record, display, and 
report the operational configurations, status, and perform- 
ance of the GCF subsystems. A central monitor processor 
has been designed to receive inputs from other GCF 
subsystems as necessary to permit internal assessment of 
performance, problem detection, isolation, and correction. 

V. Typical Configuration 

The Viking Project represented one of the most 
extensive users of the GCF because of the two orbiting 
and two lander spacecraft and extended critical phases of 
the mission. The project requirements were largely met 
with standard GCF circuit configuration as was the 
intended goal. Illustrated in Fig. 5 is the GCF configura- 
tion used during the primary mission phase following 
touchdown. 
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Table 1. High-speed data bloek composition 


Subdivision 

22-bit Error polynomial 


33-bit Error polynomial 


Word 

Bits 

Word 

Bits 

Header 

1 thru 7 (16 bits each) 

112 

1 thru 7 (10 bits each) 

112 


8 

8 

8 

8 

Text 

8 (Bit 9) thru 73 (Bit 8) 

1040 

8 (Bit 9) thru 73 (Bit 12) 

1044 

Ending 

73 (Bit 9) thru 75 (Bit 16) 

40 

73 (Bit 13) thru 75 (Bit 10) 

36 

Totals 

75 

1200 

75 

1200 


Table 2. High-speed filler block composition 

Subdivision 

22-bit Error polynomial 


33-bit Error polynomial 


Word 

Bits 

Word 

Bits 

Header 

1 thru 3 (16 bits each) 

48 

1 thru 3 (10 bits each) 

48 

Text 

4 (Bit 1) thru 73 (Bit 8) 

1112 

4 (Bit 1) thru 73 (Bit 12) 

1116 

Ending 

73 (Bit 9) thru 75 (Bit 18) 

40 

73 (Bit 13) thru 75 (Bit 16) 

36 

Totals 

75 

1200 

75 

1200 


Table 3. Wideband data block composition 


Subdivision 

Word 

Bits 

Header 

1 thru 7 (16 bits each) 

112 


8 

8 

Text 

8 (Bit 9) thru 298 (Bit 12) 

4644 

Ending 

298 (Bit 13) thru 300 (Bit 16) 

36 

Totals 

300 

4600 


Table 4. Wideband filler block composition 


Subdivision 

Word 

Bits 

Header 

1 thru 3 (16 bits each) 

48 

Text 

4 (Bit 1) thru 298 (Bit 13) 

4716 

Ending 

298 (Bit 13) thru 300 (Bit 16) 

30 

Totals 

300 

4800 
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Viking Mission Support 

D. W. Johnston and T. Howe 

DSN Network Operations Section 


This report covers the most significant Viking M,ss, on events 
DS,V during August and September 1976. Intermediate Data product, on and DSS 
support !- alto summarized for this peri od. Viking DSN Discrepancy Report 
activity tor the period January 1975 through September 1976 » also included. 


I. Introduction 

The previous report in this series described the final 
preparation and testing prior to the Mars orbit insertion of 
Viking 1, the Mars orbit insertion, separation and Landing 
of Viking Lander 1 (VL-1), and concluded with the Viking 
2 approach maneuver on 27 July 1976. This report 
continues from that point and describes the sequence o 
events pertaining to the insertion through landing ol 
Viking 2. 


II. DSN Mission Operations Activities 

A. Viking Operations Activities 

Table 1 lists the significant Viking 1 and Viking 2 
activities supported by the DSN during this reporting 
period. 

The DSN supported a total of 18 significant events 
during this two-month period. On August 1, Viking 
Orbiter 1 (VO- 1) was already in orbit, with Viking Lander 
1 on the surface of Mars sending data daily via direct and 
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relay links. Viking 2 was approaching Mars, 7 days away 
from Mars orbit insertion (MOI). 

B. DSS Support 

Table 2 lists the tracking hours, per station, of the 
Viking spacecraft and the number of commands transmit- 
ted from each station during this reporting period. 

C. Intermediate Data Record Status 

During this reporting period the DSN data record 
capability provided the Viking Project Data Support 
Group with 394 original telemetry Intermediate Data 
Records (IDR), of which 378 IDRs were delivered within 
the 24-hour after-loss-of-signal requirement. Ihe average 
data throughput for this reporting period was 99.99 
percent of the available data. 

Table 3 is a breakdown of IDR production per DSS. 

D. Viking Discrepancy Reports 

A summary of Viking Discrepancy Reports (DR) 
covering the period from 1 January 1975 to 26 September 
1976 is shown in Table 4. 


13 


I 


III. DSN Special Activity Support 

A. Dual Orbiter Tracking 

The tracking of both Viking Orbiters by a single 64-m 
DSS began on August 6. Both Orbiters had been within 
the beamwidth of 64-m antennas for the three days prior 
to Mars orbit insertion of Orbiter 2, however the Viking 
Project chose the option of pointing the antenna at a 
single spacecraft in order to maximize data return, DSN 
predicts had shown that on August 3, a signal loss of 24 dB 
would result on one spacecraft while the antenna was 
pointed directly at the other. This loss would be reduced 
to just 4 dB for both spacecraft by midpoint tracking. True 
midpoint tracking was accomplished by using the (right 
ascension /declination) planetary mode of antenna drive in 
which hour angle (HA) and declination (dec) position 
coordinates are computed by inputting HA and dec for 
the three days beginning with the day of track. Offsets 
were then used to obtain midpoint tracking. 

B. Viking 2 Mars Orbit Insertion 

The Viking 2 MOI was successfully supported on 
August 7 by Madrid DSS 63 with DSS 62 as backup. The 
same series of events associated with MOI of Viking 1 
were repeated for Viking 2. The spacecraft science 
subcarrier was turned off, then the spacecraft was rolled, 
the low-gain antenna selected, the spacecraft yawed, and 
then rolled again so that the high-gain antenna pointed at 
Earth. The high-gain antenna was then selected and the 
high-rate subcarrier turned on and modulated by 2 kb/s 
data. A loss of uplink and downlink lock occurred during 
the yaw turn. ° 

A DSS exciter ramp was started a one-way light time 
prior to the predicted uplink loss of lock. A short 
insurance sweep was then executed to agree with the time 
at which the spacecraft was commanded back to the high- 
gain antenna. This tuning assured a two-way lock for MOI 
burn. 

The MOI burn lasted approximately 39 minutes, 35 
seconds, and resulted in an orbit of 27.4 hours, just 0.2 
hour greater than the desired orbit. The periapsis altitude 
of 1519 kilometers was ach eved and exceeded the desired 
altitude by only 19 kilometers. 

Following the motor burn, the reverse order of roll, 
low-gain antenna, yaw, roll, and high-gain antenna took 
place, and an exciter ramp was accomplished to reacquire 
the uplink. 

C. Mars Orbit Trim Number 1 (VO-2) 

On August 9, Mars orbit trim number 1 (MOT-1) was 
successfully performed during the Goldstone DSS 14 view 


period. This trim refined the Orbiter 2 orbit and resulted 
in an orbit of 27.3 hours duration. 'The motor burn lasted 7 
seconds. The bum was accomplished without the benefit 
of a roll, yaw, or roll turn by using the Sun-line reference 
technique. That is, the burn was timed for the point in the 
orbit in which the attitude of the Orbiter was corrected 
without a maneuver. This method preserved the attitude 
control gas supply. Although telemetry later indicated that 
the high-gain antenna had moved 0.4 degrees during the 
motor burn, due to backlash, no changes were observed in 
the “S” or “X” band downlinks. 

D. Mars Orbit Trim Number 2 (VO-2) 

On August 14, Australia DSS 43 supported MOT-2 for 
Viking Orbiter 2. DSS 63 provided backup support for the 
motor burn. This trim maneuver was to refine the orbit of 
Orbiter 2 with a motor burn of 3 seconds. No changes in 
link quantities were observed during this Sun-line 
maneuver. These first two maneuvers provided a new 
concept in site certification. Instead of synchronizing the 
spacecraft’s periapsis with a precise point on the planet, 
the spacecraft was left out of sync so that its periapsis 
walked around the planet at 40-deg steps per day. This 
permitted a good examination of the 40- and 50-deg north 
latitude band and all proposed landing sites. 

E. Mars Orbit Trim Number 3 (VO-2) 

MOT-3 for Orbiter 2 took place during the DSS 14 pass 
on August 25, with motor burn at 17:08:45 GMT, lasting 
for 72 seconds. This was a zero-turn maneuver, and the 
telecommunications link remained in the high-rate mode, 

(2 kb/s and 8-1/3 b/s) during the maneuver. The only 
observed change during the burn w?.s a decrease of about 
0.25 dB on X-band, which would correspond to a 0.25-deg 
change in the high-gain antenna pointing angle. This 
maneuver stopped Viking 2s westward walk, changed the 
spacecraft’s speed by 155 km/h (96 mi/h) and initiated a 
drift in the easterly direction toward precise coordinates 
for the landing site. 

F. Mars Orbit Trim Number 4 (VO-2) 

The resynchronization of Orbiter 2 over the Utopia 
landing site occurred on August 27 during the DSS 14 
pass. MOT-4 again used the Sun-line maneuver procedure 
and no changes in the downlink were observed. 

G Viking Lander 2 Presoparation Checkout, Descent, 
and Landing 

Network support of activities associated with the 
landing of Viking Lander 2 began on September 2 with 
DSS 43 supporting the separation minus 39-hour com- 
mand update. This command load prepared the mated 
lander for the preseparation checkout. This command load 
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was supported using configuration code 15. Australia DSS 
44 provided backup command support for DSS 43. 

Approximately nine hours later, preseparation checkout 
began. Spain DSS 63 supported this checkout. A special 
configuration was designed for the DSS 63 pass In as 
mud) as a Lander 1 direct link was to occur during the 
first part of DSS 63’s view period, two configurations were 
required. The first configuration was the standard three- 
spacecraft configuration for Orbiter/ Lander/ Orbiter, 
Code 30. Following the Lander direct link, a configuration 
was needed that would assure the receipt of the Lander 2 
checkout data at JPL. Figure l shows the configuration 
used for the second half of the DSS 63 pass. This 
configuration provided three processing channels for 
Under checkout science data. Two of these data streams 
were output via the wideband data line while the third 
was output on the high-sneed data line, providing dual 
transmission paths. In parallel with the Under 2 checkout 
data, Orbiter 1 was outputting 4 kb/s science data. Since 
no data channel was available at DSS 63 to process these 
data they were recorded on analog tape for postpass 
playback if required. During the checkout, the relay 
subsystems were turned on and performed normally. The 
10-watt UHF transmitter output power measured 10.72 

watts, which wa, 0.44 dB above the data base value of 
9.68 watts. 

ncc he OT Lilnder 2 1 checkout proceeded normally during the 
DSS 63 pass. DSS 62 provided an uplink for Viking 
Orbiter 2. 5 

s P ain °SS 63 supported the separation minus 9.5- and 
3.5-hour command updates and the separation “GO” 
command during their pass on the following dav, 
September 3. Configuration code 24 was used to support 
his command activity. This configuration code provided 
two high-rate science data streams with channel 2 of 
telemetry and command processor alpha outputting data 
to the high-speed data line while channel 3 of the 
lelemetry and Command Processor Assembly (TCP) beta 
output data for wideband data line transmission. This same 
configuration code had been used for support of presepa- 
ration checkout of Lander 1. 


The purpose of the “GO” command was explained in 
the previous article of this series. 

All commands scheduled during the DSS 63 view period 
were successfully transmitted. 

DSS P 63 DSS ^ pr ° Vided backl, P COI ”mand support for 


Goldstone DSS 14 was prime for separation, descent 
and landing support. DSS 11 acted as the backup 
command station for DSS 14. 

Since all of the descent data were to be at 4000 b/s, 
configuration code 15 was specified. Both telemetry and 
command processor strings were initialized for Orbiter 2, 
providing dual processing channels with their outputs 
multiplexed onto one wideband data line. 

Separation occurred at 19:39:59 GMT Earth received 
time on September 3. The landing sequence had been 
perfect up to this point; however, at 1 9:46:27 GMT, the 
X-band receiver showed a drop in received signal level 
and then went out of lock at 19:47:02 GMT. At the same 
time the S-band receivers were showing decreases in 

in°?i V rn ^ GC ’ These receivers all dropped lock at 
19:47:52 GMT. Goldstone DSS 14 responded to the 
situation and regained S-band lock within 2 minutes at a 
received signal level of -169 dBm. Since this signal lev, <4 
was consistent with the low-gain antenna with the high- 
rate subcarrier on it was assumed that the signal was not 
coming from a mispointed high-gain antenna. At 20:08:34 
GMT, approximately 19 minutes after the downlink 
reacquisition, commands were sent to assure the space- 
craft was on the low-gain antenna, and to select the cruise 
mode at 8-1/3 b/s. A round trip light time later (41 
minutes) telemetry lock was established and indicated that 
a switchover of the attitude control system had occurred 
and that the spacecraft had rolled off of Vega reference by 
some unknown amount. Engineering data also indicated 
that the descending Lander-to-Orbiter real-time relay link 
was working and that the Lander appeared to be in good 
health. This mode of operation was to continue through 
touchdown (22:58:20 GMT) on September 3 and until 
approximately 07:00:00 GMT the following day. 


When it was determined the spacecraft had rolled off 
the Vega reference, a command designated “fly back and 
sweep” was sent. This command would cause the 
spacecraft to reacquire Vega if the star was within 5 
degrees of the star tracker. Vega was not acquired, 
indicating a roll of more than 5 degrees. Shortly afterward, 
commands were sent to roll the spacecraft 360 degrees to 
obtain a roll attitude reading. The Canopus loss contin- 
gency plan identified in the previous article of this series 
was used to determine the spacecraft orientation in roll 
axis based on the point at which the X-band signal level 
reached a peak. Australia DSS 43 supported this effort and 
obtained a rapid X-band lock, making the star mapping 
possible. A peak was observed of -149.06 dBm proving 
two things: one, that the high-gain antenna was still 
working, and two, that the spacecrait had rolled about 22 
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degrees off Vega. The Orbitcr was then rolled to Vega, 
the high-gain antenna selected, the high-rate subcarrier 
turned on, and the playback of the first two Lander 
pictures initiated. The pictures were received without 
further problems. 

The second received picture, which was a 310-degree 
panoramic survey, showed a blemish on the Lander’s high- 
gain X band antenna that could have been caused during 
landing. Since the antenna dish is a very thin honeycomb 
design, it was vulnerable to damage while in the stowed 
position during landing if a substantial shock was 
experienced. 

H. First Lander 2 Direct Link 

The first Lander 2 direct S-band link occurred on 
September 5 with DSS 43 supporting. Configuration code 
61 was used by DSS 43 for this pass and was to be used for 
the first 20 Lander 2 direct links. As identified in the 
previous article of this series, Code 61 provides 2 
processing channels for Orbiter data and 4 processing 
channels for Lander direct link data. 

During the DSS 43 precalibrations it was determined 
that channel 3 of 1CP beta was inoperative. This channel 
was assigned to process Lander 250-b/s backup data and 
output via the wideband data line. A decision was made to 
use channel 2 of the DSS 42 TCP string for this 250 b/s 
backup data stream. The Viking configuration permits 
receiver number 2 located at DSS 43 to interface with 
Subcarrier Demodulator Assemblies 7 and 8 at DSS 42. 
Two hours prior to the scheduled Lander acquisition, a 
problem developed with a Helios spacecraft. Australia 
DSS 42 was assigned to track Helios with the 250-b/s 
backup stream being given up. Lander acquisition 
occurred at 03:30:00 CMT with a signal level of -151.5 
dBm, very close to the predicted level. An examination of 
the telemetry data indicated that the uplink sweep had 
acquired both Lander transponders and that the uplink 
signal levels f*>r low- and high-gain antennas were nominal. 
Since all data streams appeared to be nominal, an 
assumption was made that the damage to the high-gain 
antenna revealed by the second lander picture was either 
not severe enough to degrade the communication 
performance or that the assumption of damage was 
incorrect. Because the antenna dish faces upward while 
stowed, it was possible that some surface material could 
have been thrown up onto the dish during landing and that 
the apparent damage was really a discoloration caused by 
Martian dirt. This last theory seems to be the most 
popular. Since the high-gain antenna surface would never 
again be in a position in which it could be seen by the 
Lander cameras, a closer examination was not possible. 


I. Mars Orbit Trim Number 7 (*. J* 1) 

MOT-7 took place on September 11. This trim 
maneuver was supported by DSS 14. It broke the orbital 
synchronization of Orbiter 1 and caused it to be a walk 
around the planet. Lander 1 began a reduced-mission 
phase with continuing biological experiments, weather 
reporting, and picture acquisitions, but with both real- 
time and recorded data being transmitted via the Lander 
to Earth by direct link only. Viking Orbiter 1 begain its 
walk at 12:24:34 on 11 September following an engine 
burn that lasted 16 seconds. The resulting orbit moved the 
spacecraft approximately 40 degrees east each day. The 
orbit was to carry VO-1 completely around the planet in 9 
days continuing the walk to nearly half-way around the 
planet again until it reached a plane that would pass over 
the VL-2 location. The maneuver was successfully 
completed us'ing the standard roll, yaw, roll, and burn 
sequence. 

J. Mars Orbit Trim Number 8 (VO-1) 

This trim took place on September 20 and was 
supported successfully by DSS 14. The burn occurred at 
22:36:19 GMT Earth received time, ft was a sunline 
maneuver with no changes observed in the downlink signal 
levels. Tins burn accomplished a fine adjustment to orbital 
statistics in preparation for the synchronization burn to 
occur on September 24. 

K. Mars Orbit Trim Number 9 (VO-1) 

MOT-9 synchronized VO-1 over VL-2 to permit VO-1 
to become the relay station for Lander ? data. The 
maneuver was completed using the Sun-lmc technique. 
One anomaly occurred during this maneuver. Upon 
switching to the low-gain antenna the downlink was 
observed to be about 1.5 diB low, with the uplink being at 
least 6 dB lower than predicted. Spain DSS 61 had not 
turned off range modulation a* the required time and the 
observed degradation was the result of this modulation on 
the uplink. The trim was supported by DSS 63 and DSS 14 
during an overlapping view period. 

L. Pre-MOT-5 Test Burn Maneuver (VO-2) 

This trim was accomplished for the primary purpose of 
testing VO-2’s motors prior to MOT-5 to make absolutely 
sure everything was operational prior to MOT-5. The test 
trim was supported by DSS 43. Near the beginning of 
their pass, a hydraulic pump failure occurred preventing 
movement of the DSS 43 antenna. Goldstone’s DSSs 12 
and 14 passes were extended as long as possible, and, in 
addition, commands were sent to switch the spacecraft to 
cruise engineering only mode in the event DSS 43 could 
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not make repairs in time for the maneuver and DSS 44 
would be required to track the downlink. Australia DSS 43 
managed to recover prior to the main maneuver; however, 
the cruise mode of telemetry was maintained. 

M. Mars Orbit Trim Number 5 (VO-2) 

MOT-5’s purpose was to modify VO-2's orbit plane so 
that it could begin observations of the Martian north polar 
cap. This maneuver differed from the standard in that a 
command must be transmitted and received at the 
spacecraft before the motor burn would take place. The 
spacecraft was first rolled, yawed, and rolled on the low- 
gain antenna, and then the high-gain antenna was selected. 
A measurement of the X-band received signal level was 
made. A measurement below the predicted level would 
have indicated that the spacecraft high-gain antenna was 
not properly aligned and that the burn should not be 
attempted. The spacecraft had been previously program- 
med to recover from the maneuver automatically should a 
“GO" command not be received. The X-band downlink 
was found to be within tolerance and the “GO” command 
was transmitted. Figure 2 s'bows the time line associated 
with MOT-5. The top line indicates the ground time for 
DSS initiated events. The middle line shows the events 
when they arrive at the spacecraft or events that originate 
at the spacecraft. The lower line shows events occurring a 
round-trip light time after a ground originated event or a 
one-way light time after a spacecraft event. 


MOT-5 was performed flawlessly. The motor burn 
lasted 5 minutes and produced a velocity change of 343 
m/s. It changed the orbit inclination from 55.3 degrees to 
75 degrees. The walk will change 30 degrees per day with 
a planned resynchronization over VL-2 to take place after 
18 revolutions on October 18. VO-1 will then be released 
to start a second walk during solar conjunction. It will 
resync over VL-1 sometime in mid-November. 

IV. Conclusion 

With two spacecraft in orbit and two spacecraft landed, 
all returning extremely large amounts of data daily, the 
Viking Mission would have to be considered one of the 
most successful on record. It is also gratifying to note that 
during this exceptionally high activity period, with the 
DSN being utilized to near 100 percent capability 
continuously, the loss of data due to network hardware 
failures has been insignificant. This has been made 
possible, to a large degree, by utilizing 64-meter DSS 
“failure mode" configurations that wer< lesigned 'o 
optimize the data processing capabilities at a station in 
the event of a single-point failure in any telemetry stream. 
These configurat.ons were described in detail in previous 
articles in this series. Another factor that contributed to 
the excellent network performance was the use of the DSS 
Computer Assisted Countdown, which will be described in 
the next report. 
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Table l. Viking Operations Activities 


Date 

Spacecraft 

Activity 

1. Aug. 3 

Viking l 

Station keeping trim maneuver 

2. Aug. 6 

Vikings 1 & 2 

First simultaneous tracking of 2 
VOs by one DSS 

3. Aug. 7 

Viking 2 

MOI 

4. Aug. 9 

Viking 2 

MOT-1 

5. Aug. 14 

Viking 2 

MOT-2 

0. Aug. 25 

Viking 2 

MOT-3 

7. Aug. 27 

Viking 2 

MOT-4 

8. Sept. 2 

Viking 2 

Preseparation checkout 

9, Sept. 3 

Viking 2 

Separation, descent fit touchdown 

10. Sept. 4 

Viking 2 

First direct link 

11. Sept. 11 

Viking 1 

MOT-7 

12. Sept. 18 

Viking 1 

First solar oceultation 

13. Sept. 20 

Viking 1 

MOT-8 

14. Sept. 24 

Viking 1 

MOT-9 

15. Sept. 25 , 

Viking 1 

Lunar oceultation 

16. Sept. 28 

Viking 1 

Ground Communications 
Facility 7.2 kh/s high-speed data 
line test 

17. Sept. 29 

Viking 2 

Pre-MOT 5 test burn maneuver 

18. Sept 30 

Viking 2 

MOT-5 



Table 3. 

IDR production 


DSS 

Total 

Delivered 

Delivered 

original 

in 24 hours 

late 

14 

127 

126 

1 

43 

128 

122 

6 

63 

139 

130 

9 


Table 2. Station support for Viking 


DSS 

Tracks' 1 

Hours tracked 

Commands 

transmitted 

August 

11 

30 

229:35 

2 

12 

6 

42:55 

0 

14 

69 

494:34 

2248 

42 

27 

242:53 

1440 

43 

69 

571:21 

3094 

44 

7 

56:55 

0 

61 

31 

305:24 

3511 

62 

9 

83:45 

438 

63 

62 

541 :52 

2318 

Total (monthly) 

310 

2569:14 

13,051 

September 

11 

29 

216:14 

1430 

12 

6 

45:24 

6 

14 

72 

486:59 

1205 

12 

34 

304:02 

1532 

13 

69 

584:59 

2685 

44 

11 

91:29 

4 

61 

29 

269:36 

1028 

62 

7 

98:29 

332 

63 

81 

536:10 

1495 

Total (monthly) 

338 

2633:22 

9717 

Report total 

648 

5202:30 

22,768 

•‘The number of tracks includ. the number of passes, each of 
which includes one, two, or three spacecraft simultaneously. 
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Table 4. Viking Discrepancy Report summary, 1 January 1975-26 September 1976 


Resolution 





Deep Space Station 





Deep 

Space 

Network 

Ground 

Communi- 

cation 

Facility 

Network 

Data 

Processing 

Area 

Nr work 
Operations 
Control 
A i ea 

Total 

r* 

11 

12 

14 

42 

43 

44 

61 

62 

63 

71 



Station 

Dependent 

55 

47 

188 

45 

118 

22 

3Y 

25 

141 

16 

’-r 

t 

34 

133 

72 

940 

79.8 

Station 

Independent 

6 

i 

21 

8 

10 

1 

o 

4 

18 
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Pioneer Venus 1978 Mission Support 

R. B. Miller 

TDA Mission Support 


The Tracking and Data Acquisition portion of the Ground Data System, which 
will support the Differential Long Baseline Interferometry Wind Measurement 
Experiment, is described. 


1. introduction 

i The Pioneer Venus Multiprobe Mission includes a 

Differential Long Baseline Interferometry (DLBI) Experi- 
ment, which will attempt to measure the wind velocities 
in the atmosphere of Venus as four probes descend 
through the atmosphere. The fundamentals of the 
experiment are described in Ref. 1. The experiment will 
be using interferometry techniques to measure the 
components of the wind velocity perpendicular to the 
line-of-sight (Earth-spacecraft direction), and established 
doppler techniques to measure the velocity components 
along the line-of-sight. Each pair of stations which 
comprise the interferometer resolve only one component 
of the velocity. In order to resolve both components of the 
wind velocity perpendicular to the line-of-sight and to 
provide some measure of redundancy, four stations will be 
equipped to support this experiment. The 64-meter DSN 
stations located at Goldstone, California, and Canberra, 
Australia, will be utilized, along with the 9-meter Space 
Flight Tracking and Data Network (STDN) stations 
located at Guam and Santiago, Chile. 

The end product at each tracking station is a wideband- 
width digital recording covering the entire bandwidth of 

» 

'j 22 

5 

i 

;v 


interest, which includes five spacecraft signals. Plans to 
summer 1976 had assumed that the deliverable to the 
experimenter, located at the Massachusetts Institute of 
Technology (MIT), would be the actual wideband 
recordings produced at these tracking stations. This plan 
involved the Project procuring for the experimenter 
recorders capable of playing back these wideband 
recordings, and required that the experimenter implement 
and operate a system to convert the wideband recordings 
into a useable computer-compatible form. Estimates by 
MIT of the resources required to implement and operate 
this conversion process turned out to be a very significant 
portion of the total processing costs at MIT for the 
experiment. This problem, coupled with a similar 
expensive processing cost for the initial bandwidth 
reduction of the Pioneer Venus Orbiter occultation radio 
science data, led the DSN to form a special team in May 
1976 to see what the DSN might be capable of doing to 
help solve both of these problems for the Project. This 
activity culminated in a functional design and estimated 
resources required, which were approved by both the 
Project and the Office of Tracking and Data Acquisition in 
June 1976. The result will be a new occultation subsystem 
to be implemented in the 64-meter subnet (which will be 
described in a subsequent Pioneer Venus Progress Report 
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, irti de) and the agreement of the DSN to provide with 
Project’ financial assistance) equipment in Compatibility 
Test Area 21 (CTA 21) at JPL, Pasadena, necessary for 
reducing the bandwidth of the wideband digital recordings 
to produce standard computer-compatible recordings 
which will be the deliverable to the experimenter. 

The total Ground Data System to support the DLB1 
Wind Measurement Experiment, therefore, now com- 
prises four tracking stations located around the Pacific 
Basin in Australia, Guam, Chile, and California a tape 
conversion facility located at JPL, Pasadena, Bnd - fil j a * 
the actual science data processing to be performed at the 
Massachusetts Institute of Technology. The remainder of 
this article will describe the TDA-provided portions of 
this system; the tracking station configurations and the 
tape conversion facility at CTA 21. 


II. Tracking Station DLBI Configuration 


The key requirements for the tracking stations to 
support the DLBI experiment are described in Ref. 1. ine 
details of the requirements had been conveyed in written 
material to the DSN through the flight project. In 
November 1975, a meeting was held at JPL with the 
experimenters and representatives of DSN Engineering 
and STDN Engineering for the purpose of clarifying the 
requirements face-to-face with the experimenters, After 
this meeting, the DSN then proceeded to produce a first 
cut at a functional design for the DSN stations. T is 
functional design at a block diagram level was presented 
to the experimenter in mid-April 1976 in a telecon erence. 
This teleconference and other informal exchanges e 
tween DSN Engineering and the experimenter then 
resulted in a final design for the DSN configuration which 
was formally reviewed and accepted by the experimenter 
on 22 June 1976 in a meeting at JPL. The resulting DS 
design is shown at a functional level in the top part of rig- 
1 (A more detailed description of the receiver and 
calibrator design will be provided in a subsequent Progress 
Report article.) A key decision which was made in the 
design process was that the requirements for the DLBI 
experiment were sufficiently stringent and unique o 
warrant the DLBI requirement being met by a separate 
receiver dedicated to the experiment. The original 
concept had involved one of the open-loop receivers 
which will be provided for precarrier detection telemetry 
recovery including a second port 2 MHz wide to support 
the DLBI experiment. The final design involves a separate 
DLBI receiver which brings the total complement of 
open-loop receivers at Goldstone {DSS 14) and Australia 
(DSS 43) up to five per station for the Multiprobe entry. 


The configuration involves a calibration tone generator 
which produces two tones, one above and one below the 
five spacecraft signals. The DSN configuration includes a 
switching capability to insert the calibration tones ei her 
before or after the low-noise amplifier. 1 he tones will be 
inserted before the low-noise amplifier for testing 
purposes and for pre-pass checkout, but will have to be 
Inserted after the low-noise amplifier to avoid possible 
interference with the recovery of the telemetry data from 
the probes during the actual entry event The tones are 
also fed to a monitoring system which will strip the tones 
off the output of the receiver and compare them against 
the tones direct from the^one generator. 

One of the most difficult parts of the actual receiver is 
the final mixer which must be flat from essentia ly dc to 
1.8 MHz and have the images down at least 10 dB at the 
1.8-MHz point. In order to meet the stringent relative 
phase stability requirements of the experiment, certain 
critical modules will be contained in temperature- 
stabilizing ovens. In addition, the cables carrying the tones 
from the control room up to the antenna wi 1 be actively 
stabilized as well as the cable carrying the reference 
frequency into the first mixer of the receiver Redundant 
widebandwidth digital recorders will be provided at each 
station. For this experiment, the recorders will be 
operated at a 12-megabit per second rate operating at 76 
cm (30 inches) per second. 


It has been planned from the beginning that the STDN 
design would be derived by STDN Engineering based on 
whatever was the finalized design, for the DSN stations. 
With the finalization of the DSN design on 22 June, 
STDN Engineering proceeded with the detailed design or 
the STDN stations, which culminated in a formal review 
of the STDN design on 16 September at Goddard wherein 
the experimenter did concur in the STDN design. 

Notice in Fig. 1 that the functional design of the STDN 
stations is essentially identical to the DSN except that it is 
not required to have the calibration tones switchable 
between the input and output of the low-noise junp 
since there will be no telemetry recovery from the SI UN 
stations. STDN Engineering was able to take advantage of 
the functional commonality and produce a detailed design, 
which involves use of many common modules with the 
DSN configuration. The calibration tone generator will be 
a common design using common modules as well as the 
critical last mixer mentioned above. The plan is to have 
the DSN procure, with Goddard funding, the common 
modules for the STDN stations. The STDN receiver 
involves using the first local oscillator (first mixer) in the 
standard STDN operational receiver, the multiple- 


JPL DEEP SPACE NETWORK PROGRESS REPORT 42-36 


23 


I 


frequency receiver. The analog-to-digital converter and 
recorder will be loaned to the STDN by the DSN as 
originally planned. The DSN is responsible for providing 
the necessary interface electronics to hook the recording 
subsystem up to the STDN Frequency and Timing 
Subsystem and power equipment. The interface of the 
STDN receiver to the recorder was to be a Goddard 
responsibility, but this engineering interface problem has 
been eliminated by STDN use of the same final mixer as 
in the DSN configuration, 

ill. Station Performance Validation 

There will be several means of validating the proper 
performance of elements of the station configuration for 
testing purposes and for pre-pass checkout as well as 
during the actual tyent. The calibration tone monitor, 
comparing the tones a* they go up to the input or output 
of the low-noise amplifier against the tones coming out of 
the rear-end of the receiver, will enable detection of 
several possible problems in the receiver chain. The 
recorders themselves are 18-track and have three indicator 
lights per track which are capable of detecting cata- 
strophic failures in a particular track as well as the 
capability to do a bit error computation on a single track 
at a time. Tlie micro controller on the recorder will be 
able to switch the bit error computation every 15 seconds 
from one track to another. This monitoring is possible on 
only one recorder at a time; however, the monitoring can 
be switched between the recorders without affecting the 
recording process during the period of interest. 

For a total end-to-end system check at the DSN 
stations, it is planned at the current time to provide a 
spectrum analyzer type of device to operate on a full 
read-after-write output from the recorder in real-time. 
This will enable the detection of all signals of interest, or 
at least a single signal of interest, to show that the total 
end-to-end system is configured properly and that actual 
data are going to the recorder. The possibility of 
employing a similar process at fhe STDN stations has 
been discussed, and it appears feasible to provide a 
commercially available spectrum analyzer capable of 
detecting the bus signal using a digital-to-analog converter 
after the full read-after-write on the recorder. The 
advantage of this kind of approach for the STDN stations 
is that it would also be assuring that the station was 
actually on point. This is known automat* -ally at the 64- 
meter stations because of the capabii»t> detecting the 
telemetry in real-time. Other alternatives are being 
discussed for the STDN stations, such as the feasibility of 
tracking the Orbiter spacecraft using a regular closed-loop 
receiver just for the purpose of validating that the station 


is indeed on point. These monitor alternatives are still 
being evaluated by the STDN. 

IV. The DLBI Tape Conversion Facility 

The bottom part of Fig. 1 portrays the functional 
elements in the tape conversion facility to be provided at 
the Compatibility Test Area {CTA 21) at JPL, Pasadena. 
The process involves recovering one spacecraft signal at a 
time from the wideband recording by digitally mixing an 
estimate of that particular spacecraft signal with the 
output of the wideband tape recording and then digitally 
filtering the product so that the output is of sufficiently 
narrow bandwidth to be accommodated on a standard 
9-track computer-compatible tape recording. The digital 
bandwidth reduction mixer and filter is envisioned to be a 
piece of special-purpose hardware (involving firmware), 
while the controller and formatter will be software 
written for a standard DSN minicomputer already 
available in CTA 21. The planned resulting bandwidth will 
be on the order of 1 kHz. Each of the five spacecraft 
signals (four probes and the bus) and the two calibration 
tones will each be independently processed. It will be 
Project responsibility to provide the DSN with the 
estimate of the frequency history of each of the five 
spacecraft signals. 

V. DLBI Ground Data System Integration 
Plan 

Figure 2 shows the overall major milestone schedule for 
the DLBI Wind Measurement Experiment of Pioneer 
Venus. For integrating the DSN and STDN elements, it is 
planned to do a trial installation in an STDN station 
located at Goldstone of all equipment required for the 
STDN configuration in mid-November 1977. This should 
enable solving any final integration and interface problems 
utilizing local engineering support. For interfacing with 
the experimenter, it is planned to provfJe on 1 October 
1977 a computer-compatible tape that is correct in 
format, followed on 1 January 1978 by a computer- 
compatible tape which contains actual tones which have 
passed through the entire TDA portion of the Ground 
Data System. It is anticipated that the first computer- 
compatible recording containing actual signals (from the 
ALSEP or some other suitable spacecraft signal) will be 
available for the experimenter at the end of the first week 
of February 1978. At least two of the four stations should 
be ready in mid-February 1978 to start supporting some 
end-to-end system checkouts, with all four stations ready 
by mid-March 1978. It is planned to be fully operational 
by 1 April 1978, and to start supporting regular 
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interferometry observations (currently envisioned as 
ALSEP tracking) from that time until the actual entry 
event the first of December 1978. This will give at least 


eight months to shake out any remaining problems in the 
total system, including the equipment and software at 
MIT. 


Reference 


1 Miller, R. B., “Pioneer Venus 1978 Mission Support,” in The Deep Space 
Network Progress Report 42-31, pp. 11-14, Jet Propulsion Laboratory, 
Pasadena, Calif., Feb. 15, 1976. 
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Ffg. 2. Pioneer Venus 1978 major milestone schedule for DLBI Wind Measurement Experiment 
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Helios Mission Support 

P. S. Goodwin 
TDA Mission Support 

E. S. Burke and R. E. Morris 
DSN Network Operations Section 


Operating beyond its 18-month design lifetime, the Helios- 1 spacecraft 
continues to gather valuable scientific data. Having achieved its fourth 
perihelion, the Helios-1 spacecraft is being closely followed by Helios-2, as the 
latter emerges from a lengthy solar conjunction period, in which considerable 
scientific data were collected on Experiments 11 and 12 ( Celestial Mechanics and 
raradatj Rotation). The scientific opportunities to compare the data from the 
two Helios spacecraft continue to add to mans knowledge of the suns influence 
upon the inner solar system. 


1. Introduction 

This is the twelfth article in a series that discusses 
Helios-1 and -2 mission support. The previous article (Ref. 
1) reported on a Helios-1 spacecraft power anomaly, 
Helios-2 occultation, STDN-OSN telemetry cross-support 
engineering-level tests, tracking coverage and DSN 
performance. This article covers the Helios-1 third solar 
occultation and fourth perihelion. A Helios-2 spacecraft 
emergency and third solar conjunction (occultation) are 
also discussed. Additional topics include STDN-DSN 
telemetry and command cross-support. Ground Data 
System tests, tracking coverage, plus DSN system 
performance for August and September 1976. 
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II. Mission Operations and Status 

A. Helios-1 Operations 

The Helios-1 spacecraft occulted the sun for a third 
time on September 23, 1976. The spacecraft was 
transmitting on traveling-wave tube amplifier 2 (TWTA- 
2), high-power mode, and stored experiment data on-board 
the spacecraft for memory readout later at a higher bit 
rate. As the probe cleared the “blackout” region, heading 
towards its fourth perihelion, the spacecraft team 
commanded TWTA-2 to medium power. The medium 
power mode was selected to reduce the risk of excessive 
traveling-wave tube (TWT) assembly temperature as the 
spacecraft approached the sun. 
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At 1907 GMT on Tuesday, October 5, 1976, the Helios- 
1 spacecraft successfully passed its fourth perihelion at a 
minimum distance of 0.309 AU (approximately 45 million 
km) from the sun. The maximum solar impact was equal 
to 10.445 solar constants. The overall performance of the 
spacecraft and its 10 scientific instruments, which had 
exceeded the design lifetime of 18 months, was excellent. 
All 10 experiments were fully configured for prune 
mission mode and delivering valuable science data. 


Several members of the scientific team observed the 
perihelion passage from the German Space Operations 
Center (GSOC). The scientists are finding Helios- 1 
perihelion data very interesting since it provides a second 
unique opportunity of comparing measurements with 
Helios-2 (presently approaching its second perihelion, 
which is centered around October 20, 1976). 


B. Helios-2 Operations 

At 17:44:52 GMT on September 4, 1978, after 
transmission of a “ranging-on” command (at 1715 GM1), 
the downlink of the Helios-2 spacecraft was ost This 
occurred during pass number 234 over DSS 11 (Gold 
stone). The failure was observed while attempting to 
acquire the ranging signal. Several commands were 
transmitted altering the spacecraft's RF configuration in 
an attempt to restore the downlink. Tins eventually 
resulted in reestablishing the Helios-2 spacecrafts down- 
link, but near threshold. To learn if the low signal was due 
to DSS 11 or the spacecraft, the services of Dbi> 4/ 
(Australia) were attained from a scheduled Viking pass. 
DSS 42 was brought up as a Helios follow-on pass a er 
DSS ll’s loss of signal. Subsequent troubleshooting 
(altering the spacecraft configuration by commands) over 
DSS 42 tentatively restored the downlink signal at 
GMT in the following configuration: TWTA-1, ig 
power, driver 1, and noncoherent mode. A spacecraft 
emergency was declared at 2254 GMT. After reranfigur- 
iny the spacecraft in the combination of TWTA-1, 
medium power, driver 1, coherent and medium-gain 
antenna the DSS 42 receiver obtained lock. Switching to 
the spacecraft’s high-gain antenna (HGA) re-sulted m a 
solid receiver and telemetry lock at 0121 GMT, Septem- 
ber 5. The spacecraft configuration and health were 
verified, and the spacecraft emergency terminated at 0500 
GMT on September 5. 

After a few days in which to analyze the spacecraft 
anomaly, it was theorized that the loss of downlink was 
probably caused as a result of two commands: (1) rangmg- 
on, and (2) TWTA-2 medium power. Pending further 


investigations it was recommended to (1) avoid ranging 
commands until after the second perihelion, and (2) stay 
on TWTA-1 in the medium-power mode. 

Helios-2 experienced its third “blackout” (solar occulta- 
tion) period during September 20-29, centered on 
September 26, 1976. All scientific data were read into 
memory on-board the spacecraft and later read-out when 
higher bit rates were available. 

The Helios-2 spacecraft will pass its second perihelion 
around October 20, 1976. The results of this perihelion 
will be covered in the next progress report. 


C. STDN-DSN Engineering Test and Support 

Encouraged by the DSN engineering tests regarding 
STDN real-time telemetry and command cross-support 
completed during May-June 1976 (Refs. 1 and 2), the U. 
Helios Project Manager requested support for Hehos-l 
-uid -2 by the STDN Goldstone Station (CDS). This cross- 
support was requested for the October perihelion* of 
Helios-1 and -2, thence continuing until November 15, 
when STDN GDS will temporarily be de-comrnitted from 
all flight support for a scheduled equipment modification 

program. 

Before the STDN-DSN real-time Helios cross-support 
could be considered operational, the link between S1DN 
GDS and the DSN Goldstone inter-site microwave system 
needed more permanency. The interface installed for the 
engineering tests in June was a one-way temporary coaxia 
cable. 


After considering available options, it was decided that 
the most suitable and economical PJ 0V ‘ .“J® “ 

full duplex interconnection between STDN GDS and DSS 
11 was to reestablish the old Project Apollo microwave 
link between STDN GDS and DSS 14. There it could join 
the DSN Goldstone inter-site microwave system. I he 
Apollo-era equipment was reinstalled and successful y 
tested by September 30, 1976. 

The plan was to first send real-time spacecraft 
telemetry over the newly established microwave link (Fig. 
1) from STDN GDS to DSS 11 (via DSS 14 and GCF 10, 
both at Goldstone) where it would be processed by the 
back-up TCP string and sent to JPL. The uplink-lock and 
command tests would be conducted a week later followed 
by a telemetry and command demonstration pass. Upon 
successful completion of this demonstration pass, the 
STDN-DSN Helios cross-support configuration would be 
committed for operational support. 
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The first real-time Helios spacecraft telemetry was 
processed over the new STDN-DSN configuration (see 
Figs. 1 through 3) on October 1, 1976, during a 
demonstration pass. Data from two other STDN-DSN 
Helios cross-support passes were accumulated on October 
15 and 16. The results are shown in Table i. 

The STDN-DSN Helios command configuration was 
successfully tested during the week of October 11, 1976. 
Several spacecraft uplink acquisitions and commanding 
sequences were conducted for STDN (GDS) operator 
training during this week. With the successful STDN- 
DSN Helios Telemetry and Command cross-support 
demonstration track performed on October 18, 1976, the 
STDN-DSN cross-support configuration was declared 
operational for Helios Project support. 

Although the data from this cross-support configuration 
were degraded (approximately 3.5 dB) when compared to 
a DSN 26-meter station, the data received from Helios at 
perihelion and distances less than 1 AU are excellent. 

D. Ground Data System Tests 

1. High-Speed Data Line Change to 7.2 kb/s. The Deep 
Space Network (DSN) plans to convert all high-speed data 
terminals to 7.2 kb/s on December 1, 1976. This 
conversion will standardize DSN high-speed data line rates 
with NASA Communications (NASCOM) standards. On 
September 27, 1976, during Helios-1 pass 658 over DSS 12 
(Goldstone), a special Ground Communications Facility 
(GCF) test was performed to verify that the DSN GCF 
conversion to a 7.2-kb/s high-speed data line (HSDL) rate 
was compatible with mission-dependent Mission Control 
and Computing Center (MCCC) software programs and 
hardware. The test consisted of sending parallel data 
streams on standard 4.8-kb/s HSDL and the new 7.2-kb/s 
HSDL to MCCC for Model 360-75 computer processing 
(monitor and tracking data were not tested at this time). 
Although one command anomaly did occur, there were no 
problems attributable to the 7.2-kb/s interface during the 
test. The planned changeover to the 7.2-kb/s HSD system 
on December 1, 1976 should be essentially “invisible” to 
the MCCC data system and to DSN-MCCC operations. 

2. New Telemetry Capability at the German Weilheim 
Station. In April 1976, the German Helios Project Office 
received word that their request to modify the Telecom- 
mand Station at Weilheim to include a telemetry data 
receiving capability had been approved. Conversion was 
planned for late 1976 (Ref. 2). By August the conversion 
was progressing quite satisfactorily. The telemetry 
receiving equipment had been transported from the 100- 
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meter Eflfelsberg Station to Weilheim. After internal 
(German) engineering tests were completed, internetwork 
interface tracking validation tests were scheduled involv- 
ing Weilheim and DSS 62 (Spain). The DSN supported 
two of these tests— the first on September 13 and the 
second on September 17, 1976. These tests were highly 
successful, and the Weilheim Telemctry/Cominand 
Station (67/68) became operational on September 20, 
1976. 

E. Actual Coverage Versus Scheduled Coverage 

This report covers tracking activities for a 63-day 
period from August 9 through October 10, 1976. 
Operations for this period include Helios-1 cruise, solar 
occultation, and fourth perihelion. Helios-2, still in solar 
conjunction phase, occulted the sun for the third time. 
Total tracks for both spacecraft were 101 for a total of 
862.7 hours. Although this represents only 82 percent of 
the passes taken last period, there were 61 more hours of 
coverage. The Helios spacecraft received 41.4 percent of 
the DSN tracking time allotted to both Pioneer and 
Helioy, after Viking requirements were satisfied. Helios- 1 
was tracked 50 times for a total of 441.4 hours; Helios-2 
was tracked 51 times for 421.6 hours. The average pass 
time was 8.8 hours for Helios-1 and 8.2 hours for Helios-2. 
This compares with the last report period of 7.03 and 5.5 
hours, respectively. Only 25 hours of 64-meter subnet 
support was allocated for the Helios spacecraft during this 
period due to Viking support requirements. As the Helios- 
1 and -2 spacecraft trajectories approach Earth, only the 
26-meter subnet will be required for tracking, with the 
exception of ranging at DSS 14. DSN tracking support is 
expected to continue at its present level throughout the 
next report period. 

III. DSN System Performance 

A. Command System 

The Helios extended mission command activity contin- 
ued at a low level during August and September for both 
spacecraft. A total of 2461 commands were radiated to 
both Helios-1 and -2 spacenraf,*. This compares to 2910 
commands radiated during the last .eporting period. This 
low activity is due to two reasons: (1) the Helios-2 
spacecraft spent the greater portion of this time in solar 
conjunction (sun-Earth-probe (SEP) angle less than 5 
degrees), and (2) heavy Viking activities resulting in 
relatively few tracks for Helios. The cumulative command 
totals through September are 40,792 for Helios- 1 and 
14,345 for Helios-2. 
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There were three command system aborts during 
August and September— all with Helios-1. Two of the 
aborts were caused by DSS transmitter failures and the 
third by a fault in a DSS Command Modulator Assembly 
(CM A). Helios-2 experienced no command aborts. The 
cumulative command system abort count was raised to 13 
for Helios- 1 and remained at 3 for Helios-2. 

Total command system downtime for the months of 
August and September rose to 12 hours and 50 minutes 
(almost twice the figure for the previous 2 months). 
Approximately one-half (5.3 hours) of the reported 
downtime for this period was caused by one failure at DSS 
44 (a faulty circuit breaker). Of the 10 failures recorded, 5 
were transmitter associcted. 

B. Tracking System 

As the Helios spacecraft speed toward their even- 
numbered (in-bound) perihelions, their velocity with 
respect to Earth causes a negative doppler shift which 
exceeds the capability of the standard doppler extractor. 
To compensate for this condition, an additional 1-mega- 
hertz bias is added to the standard doppler bias at the 26- 
meter stations for a total bias of 2 megahertz. (Stations 
with Block IV receivers (64-meter'i are able to select a 
minus 1-MHz bias.) The bias shift is required to avoid 
ambiguous doppler counts when the negative dopier 
approaches or exceeds the amount of positive doppler 
fixed bias. 

The Network Operations Analysis Group (Tracking) 
computes the dates to add this special bias into the 
systems for each spacecraft. This doppler “offset” bias will 
be required for Helios-1 from September 29 through 
October 26 and for Helios-2 from October 6 through 
November 7, 1976. 

Also affected are the transmit and receive frequencies at 
the Helios-supporting DSN stations. Because of this 
negative doppler frequency, a ground station must 
transmit at a lower than normal frequency and receive a 
higher than normal frequency to maintain relative 
constant frequencies at the spacecraft receiver. On 
September 23 the DSN stations supporting Helios-1 began 
operating on voltage-controlled oscillator (VCO) transmit 
S-band channel 20-b and receive S-band channel 22-a, as 
opposed to channels 21-a and 21-b, which are normally 
used. This will continue until November 2, 1976. Helios-2 
VCO change to I4-b and lff-a for transmit and receive, 
respectively, will be required between September 27 and 
November 14, 1976. The times to switch to these 
frequencies are computed by the Network Operations 
Analysis Group. 
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The Helios-1 spacecraft’s ranging transponder began to 
operate once more on July 14, 1976, when its temperature 
again came within operable limits (Ref. 3). Good ranging 
points are now possible with Helios- 1 until about October 
20 when it is predicted that the temperature of the 
ranging unit will again slip into the non-operational range 
(5°C to 18°C). 

Ranging with the Helios-2 spacecraft has not been 
attempted since the spacecraft emergency on September 
4, 1976. Because a ranging command was thought to be 
partially to blame for the loss of spacecraft downlink, it 
was agreed that ranging would not be tried again with 
Helios-2 until after the second perihelion on October 17, 
1976, when better link conditions would exist. 

C. Telemetry System 

During the Helios-1 and -2 superior conjunction periods, 
two types of data have been routinely collected for 
internal DSN use. The first type of data consists of hourly 
readings from the system-noise-temperature (SNT) record- 
ings. The second type was spectral broadening test data. 
These data are being used to expand a data base that will 
be used later to improve telemetry predict values at small 
SEP angles. 

Telemetry support for Project Helios has been rather 
routine— broken only with special downlink signal predic- 
tions for the STDN-DSN cross-support configuration. 
DSN telemetry predicts have been within 0.6 dB in 
automatic gain control (AGC) and 0.8 dB in signal-to-noise 
ratio (SNR) for Helios-1. Helios-2 figures are not so readily 
computed since tracking has beer, minimal during this 
period, caused by extremely small SEF angles. 


IV. Conclusions 

Emerging from a relatively quiet cruise phase, the 
Helios-2 spacecraft experienced its third solar occultation 
on September 23, 1976. Spacecraft ranging was again 
possible as the spacecraft ranging assembly temperature 
rose above 18 degrees Celsius. On October 5, 1976, the 
Helios-1 spacecraft passed within 0.309 AU (approxi- 
mately 45 million km) of the sun as its instruments 
recorded the fourth perihelion. The spacecraft s health is 
good as it continues its trajectory which will bring it closer 
to Earth. 

Overcoming a spacecraft emergency on September 4, 
1976, the Helios-2 spacecraft passed its third solar 
occultation on September 26, 1976. Owing to the nature 
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of the spacecraft emergency, the spacecraft team is 
cautiously delaying, until after the second perihelion 
(October 20, 1976), any further attempt at spacecraft 
ranging. 

The awaited decision confirming STDN-DSN Helios 
cross-support for Project Helios came in mid-September 
1976. An STDN-DSN Helios cross-support configuration 
was verified and placed into operation on October 18, 
1976. This telemr'rv and command configuration utilizes 
the facilities shov* 1 1 in Figs. 1 through 3. 

The DSN is continuing plans to convert to the 7.2-kb/s 
high-speed data terminals. One of the tests performed to 
verify hardware/software compatibility was run during 
Helios-1 pass 658 in September. All equipment performed 
well. 

The German rieiios command station at Weilheim was 
modified to include telemetry and command capabilities 


(telemetry added) and began operations on September 20, 
1976, The German Effelsberg 100-meter radio telescope is 
no longer supporting Helios. 

The Helios-1 and -2 spacecraft tracking activities 
remained at about the same level as in the previous 
period. This was primarily because the Helios-2 spacecraft 
remained in solar conjunction for most of this time period. 
Tracking coverage should increase during the next period 
owing to both spacecraft passing perihelion and approach- 
ing closer to Earth. 

Although the number of commands sent to the Helios 
spacecraft decreased, there were three command system 
aborts during this report period (all on Helios-passes). 
Except for transmit and receive frequency changes and 
added doppler bias to compensate for the high doppler 
rate at this point in the trajectories of the spacecraft, the 
performance of the Telemetry and Tracking Systems 
remained fairly routine. 
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Table 1. STDN-DSN Hellos cross-support engineering test results 


Day of year 


Predicted 
(DSN), dl) 

STDN 


DSN vs STDN 

(1976) 

Bit rate 

Block III, dB 

MFR, dB 

Block III 
difference, dB 

MFK 

difference, dB 

275 

(Oct 1) 

32 

9.8 

6,2 

4.8 

3.6 

5.0 

64 

8.5 

4.0 

2.6 

4.5 

5.9 

289 & 290 
(Oct 15, 16) 

128 

8.1 

5.3 

N/A 

2.8 

N/A 
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Fig. 1. STDN (GDS)-DSS microwave configuration 
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Fig. 2. STDN (GDS) cross-support configuration 



Fig. 3. DSS cross-support configuration 
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Terminology of Ranging Measurements 
and DSS Calibrations 

T. Komarek 

Telecommunications Systems Section 
T. Otoshi 

Radio Frequency and Microwave Subsystems Section 


This article presents a proposed set of basic terminology related to deepn space 
ranging measurements. Calibration equations are derived for the “ dish-mounted 
zero delay device ” method for 26-m antenna systems and the “ translator ” method 
for 64-m antenna systems. 


I. Introduction 

Data produced by the ranging system can be used in 
many ways: spacecraft orbit determination, space plasma 
measurements, earthquake detection, etc. Whatever the 
objective is. the basic element is always the measurement 
of the time elapsed between transmission and reception of 
a ranging code by an Earth-bound Deep Space Station 
(DSS). The ranging code is modulated on a RF carrier 
transmitted from a DSS to a distant spacecraft (S/C), 
where it is demodulated and retransmitted back to Earth, 
Tlie time delay measured, therefore, necessarily includes 
the DSS and spacecraft hardware delays, as well as 
medium effects due to tropospheric refraction and “decel- 
eration” in the ionosphere and space plasma. 

This article addresses terminology of ranging measure- 
ments with obvious emphasis on DSS delay calibrations 


since this is the most complicated part of the ranging 
telecommunication link. To make the semantics meaning- 
ful. it is useful to recognize the following “axioms:” 

( 1) Range is always one-way; it is measured in units of 
length. 

(2) Delay is measured in time units (whatever they 
are).* 

(3) “Light time” corresponds to propagation with the 
velocity of light in vacuum. 

This article was originally a response to the immediate 
needs of the Vikirig Radio Science Team. It is suggested, 
however, that the terminology presented would be con- 


*The "range unit" is a time unit since it is an integer fraction of the 
time period of the DSS master oscillator, 
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sidered for acceptance by the Ranging Accuracy Team 
and, through this body, acceptance by other organizations 
involved, i.e., DSN, Spacecraft Telecommunications, 
Navigation, and Mission Control. 


II. Definitions and Explanations 

The following definitions apply to the general case 
ranging system geometry shown in Fig. 1: 

Topocentric range p, in meters, is the distance between 
the DSS and spacecraft reference locations. 

Reference locations are the spacecraft antenna aperture 
plane, DSS antenna axis crossing for the azimuth- 
elevation (az-el) mount, and the intersection of the polar 
axis with the hour angle (HA) plane for the hour angle- 
declination (HA-Dec) mount. 

Round trip light time , in seconds, is 
RTLT = 2 p/c 

where c is the velocity of light in vacuum. 

Round trip propagation time is the time for the ranging 
modulation to propagate from the DSS reference loca- 
tion through the propagation media (planetary atmo- 
sphere, space plasma, . . .) to the S/C reference location 
and back, as measured by a correlation-type ranging 
machine. 


Zcorrkctiox ~ 2-way ground station delay correction 
which must be made to account for the 

(a) difference in the DSS turnaround 
device path and the ranging-to- 
spacecraft path up to the DSS 
aperture plane. 

(b) difference between the DSS an- 
tenna aperture plane and the DSS 
reference location. 

888. DSS Calibrations 

A. Dish-Mounted Zero-Delay-Device Method 

Station delay calibrations by means of a dish-mounted 
zero delay device (ZDD) are presently used by the 26-m 
DSS subnets. This technique (Fig. 2) uses the geometric 
optics principle in which, for a large antenna with a 
parabolic geometry (including Cassegrain S/X reflex feed 
systems), the distance for any ray path from the primary 
feed phase center to the paraboloid aperture plane via the 
microwave optics airpath is a constant. Under this assump- 
tion, the time delay for any ray propagating along the 
microwave optics path from the feed horn phase center 
to the aperture plane is the same as that for any other ray 
path between these two points (going from the feed horn 
phase center to the aperture plane or from the aperture 
plane to the feed horn phase center via the microwave 
optics ray paths). 


RTPT ~ RTLT + BIAS uedia Assuming that the internal delay of the ZDD is negligi- 

ble, it can be seen from Fig. 2 that when ranging to a 
spacecraft 

Delay to spacecraft, D s/C , in seconds, is the total delay 
measured by a ranging machine when ranging to a 

spacecraft. It is related to RTPT by Ds/C ~ Dzdd + T * + t r.«p + BIAS 8/C + r Rt down + Tjr 


RTPT — D 8/c BIAS b/c ~ BIAS D8S + Zcorrection 

( 1 ) 

where all entries are in seconds, and 
BI ASs/c — spacecraft turnaround delay 


(2) 

where D zdd is the round trip delay from the ranging ma- 
chine to the ZDD located on the dish surface and back to 
the ranging machine. It includes transmitter and receiver 
delays as well as’ round trip airpath delays within the 
antenna structure. 


BIAS D ss — 2-way DSS delay pre- and post-track 
measured from the ranging machine to 
a DSS turnaround device (either ZDD 
or translator) and back. 


It can also be seen from Fig. 2 that 

RTPT — Tft + Th + r X + T R<U p + r Ri down + Tx + Tfc + T6 

( 3 ) 
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SoH »ng for + t*,, from (2) and substitution into 
(3) gives 

Rl P7 = Ds,c ~ BIAS,,,- — Dzdu + 2r/, + 2to (4) 

Comparison of (4) with (i) gives 


BIAS OSS ™ Dy.oo 

( 5 ) 

ZcORK KCTinS = 2 Tk + 2 TI, 

(6) 

where, from the geometry of Fig. 2, the 

one-way airpath 

delays ti, and to (in seconds) are 


h 

T h 

c 

(V 

b 


7 1 — cosS 

c 

(8) 


and where 


h — distance in meters between two planes parallel to 
the aperture plane. The first plane intersects the 
ZDD and the second plane contains the declina- 
tion (secondary) axis. 

b = perpendicular distance in meters between the hour- 
angle (primary) and declination (secondary) axes. 
This distance in JPL antenna drawings is shown to 
be 6.7056 meters (22.0 ft). 

c = speed of electromagnetic wave in vacuum (2.997925 
X lO'ins 1 ) 

S = spacecraft declination angle 

1 lie values of Dy nn given by (5) are reported for pre- 
and post-track on the “DSS Post-Track Report.” The 
ZroKHKmos given by (6) depends not only on the physical 
location of the ZDD on the dish surface but also on the 
spacecraft declination angle 8, Values of h and b are 
obtained at eaeh 26-m DSS from physical measurements 
and are the constants supplied by the DSN. These con- 
stants together with the declination angle information 
should be used by the Orbit Determination Software to 
determine the Zcorkei’tiox from (6) for the appropriate 
tracking pass. 

B. Translator Method 

The DSS calibration method using the dish-mounted 
ZDD suffers by multipath errors due to multiple reflec- 
tions from the antenna structure. For 64-m stations, where 


the multipath error can be several meters, another calibra- 
tion method was developed using the Block IV translator 
(Fig. 3). In this method, only a portion of the DSS hard- 
ware delays is pre- and post- track measured. The rest 
of the hardware delays (including antenna and cabling) 
is calculated from geometrical dimensions (in the case of 
waveguide and airpath) and precalibrated (in the case of 
complex components such as diplexers and filters). 

When ranging to a spacecraft, the ranging machine 
measures D.syr, the total two-way delay-to-spacecraft. 
D, /r contains the propagation delays and both DSS and 
S/C biases. From Fig. 3, using t to denote one-way delay, 
it can be seen that 


D s/r — T, + 7.! + 7 C , i,,, + T l{tUV -f BIAS,,,. + 


+ 


TRtdown * Tcv/owh * T4 T 7*2 


( 9 ) 


where 

all entries are in time units 

T i. 7 2 , 7., — DSS hardware delays shown in Fig. 3. 

7n = propagation delay (including medium effects) 
from tl\e DSS antenna aperture plane to the S/C 
reference location. 

T c — one-way delay from feedhorn phase center to 
aperture plane of the DSS antenna via the micro- 
wave optics airpath. The S-band airpath includes 
points on the ellipsoidal reflector, diehroic plate, 
hyperboloidal subreflector, and the main para- 
boloidal reflector surface. 

7/, = one-way delay between the aperture plane and 
DSS reference location 

Both r ( . and r {) are constants for a given DSS, independent 
of the point on the dish surface where reflection takes 
place, due to the parabolic geometry of DSS antenna 
dishes and Cassegrain design. 

Also from Fig. 3, it can be seen that the round trip 
propagation time is 

RTPT — TO + 7[{,up + 7 II, (town + TO (10) 

But from (9) 


7k, up "b Tlt,<loWn ~~ Ds/C 



BIASg/o 


7(', up 7 dot/i„ 


( 11 ) 
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Substituting (ii) into (10) gives 

RTPT = 2 "> + D - - £ r, - rr.., ~ 

I I 

- BIASa.c 

Defining one-way delays Ts and T . as 
A 

““ T r,tt}> ~ ~r> 


T r.riown 


( 12 ) 


1 11 Tc,ii<nrn " 7 


I), 


then substituting into (12) leads to the RTPT ex nr o« ion 
when ranging to spacecraft of ' P 

RTPT = Da,c - BIAS S/C ~y] t, 

ill 

''■hen ranging to translator, the ranging machine met 
lhe two-way delay-to-translator. 


(13) 


measures 


From Fig. 3 


Djclatok — 7, r xi.. n 


+ To 


(14) 


vvhe 


re 


\L.\roK = one 


-way delay from the uplink sample point 
"’tough the translator to the downlink 
sample point. This delay is calibrated by a 
portable ZDD. y 


From (13) and (14) 

RTPT = D„ c - BlAS s v - D xtAn , - ^ 

i :t 

as; 


BIASiuut - Dx,. \rt)tt 

o 

ZrouiiKtnox = Txt.xroit — T , 

or from substitutions of T , and T „ 


i-3 


(16 

(17) 


Z, 


OIUiKt'TtOX 


'XL.wnn + 2r/j — t., — - 

* T (\Up 7 (\(IotOn 

(17a 


where all terms are defined in Fig. 3. 

» avf i I (Dxutoh) IS a measured delay 

K ported by a station pro- and post-track. Values of 

for rj" ‘ mC <TrMn) f ° r the S/X ri ‘ fle * feed configuration 
tor 64-tn antenna syste.ns have been determtoed by 
Batelaan and are reported in Ref. 1. 


Reference 

Otoshi, T. Y., ed., *4 Collection of Articles on s/Y / r 

Ranging Tests, Technical Memorandum 33 747 ^ Experiment Zero Delay 

tory, Pasadena, Calif,, Nov. 1975 , p , 50 . ’ 4 ’ '° ' 1 Pro P ulsion Labors- 
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D s ^ c - MEASURED DELAY -TO-SPACECRAFT 
MEASURED BIAS D5S 




© 


DSS REFERENCE 
LOCATION 



CORRECTION 

DS5 ANTENNA 

APERTURE 
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TOPOCENTRIC RANGE p 


PROPAGATION VELOCITY = VELOCITY OF LIGHT IN VACUUM (C) 

— VW- PROPAGATION VELOCITY < C 
RTPT - D s/C - BIAS S/C - BIAS dss + ^CORRECTION 


Fig. 1. Ranging measurement definitions 



r = ONE-WAY DELAY 

D - ROUND TRIP DELAY TO ZDD MEASURED 
ZDD BY RANGING MACHINE 


Fig. 2. DSS calibration — dish-mounted ZDD method 
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Magnitude of 64-m Elevation Axis Movements 
Due to Alidade Temperature Changes 

N. t. Hung, H, Phillips, and R. Zanteson 

DSN Engineering Section 


In projected very dong baseline interferometry (VLBI) work a reference point 
for 64-m antennas is the intersection of the elevation and azimuth axes. This 
report describes a minimum-level effort at DSS 14 to determine the nuignitude 
of the effect of diurnal changes in the temperature of the alidade legs on the 
height of the elevation axis. The thermal expansion between the lowest recorded 
temperature —3°C (27°F) and the highest, 36°C (97 °F), over the period covered 
was 8.9 mm (0.35 in.). 


I. Introduction 

In projected VLBI work a reference point for 64-m 
antennas is the intersection of the elevation and azimuth 
axes. This report describes a minimum level effort at 
DSS 14 to determine the magnitude of the effect of 
diurnal changes in the temperature of the alidade legs 
on the height of the elevation axis. The results, over a 
period from March through September, 1976, show an 
average height change of 1.70 mm (0.07 in.) during a day, 
with a range from 0.1 mm (0.004 in.) to 3.5 mm (0.14 in.) 
for the period covered. The thermal expansion between 
the lowest recorded temperature — 3°C (27°F) and the 
highest, 36°C (97°F), over the period covered was 
8,9 mm (0.35 in.). 
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II. Test Instrumentation 

Alidade temperature measurements were made using 
surface thermometers placed on an inner surface (to 
avoid direct sunlight) of the front alidade legs (Figs. 1 
and 2). Initially, thermometers were placed on the right 
front leg, After one month’s data were taken, one of the 
thermometers was moved to the left leg to determine if 
the temperatures there were significantly different, 

The alidade leg and ambient temperature readings 
were recorded three times each day by maintenance per- 
sonnel, generally at 0600, 1200, and 2000 hours. The 
temperature readings for the right leg were averaged, 
and the differences between the high and low reading for 
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these data and for the left leg and ambient tempera- 
tures were tabulated. The results were then reduced 
statistically, and are shown for each month in Table 1. 
The dimensional change is based on a length of 20.96 m 
(835 in.), the distance from the elevation bearing to the 
machinery deck level, and a coefficient of thermal expan- 
sion of 1.08 X 10 V°C (6 X 10~V°F). 


III. Conclusions 

The maximum elevation change of the elevation axis 
measured during one day was 3.5 mm (0.14 in.). The 
maximum change over the period of record was 8.9 mm 
(0.35 in.), These numbers may be useful in evaluating 


methods for monitoring the position of the axis inter- 
section. 

The precision of the results obtained is seriously 
limited by the limited number of data points for each 
day. Tlie use of continuous recording equipment would 
provide much more accurate information, but at a con- 
siderable cost in instrumentation and data reduction. The 
need for this must he defined before such a program is 
undertaken. 

An accurate correlation between temperature measure- 
ments and elevation position was beyond the scope of 
this effort. The use of temperature measurements as an 
element in determining the position of the elevation axis 
would require considerable development effort, if 
practical. 
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Table 1. Changes in alidade leg temperature, length, and ambient temperature 


Month, 

1970 


Night log 


Loft log 


Ambient 


AT, C Ah, nun Average AT, C Ah, nun AT, °C 

(standard deviation) (for average Af) (standard deviation) (for average At) {standard deviation) 


Mar. 


Apr. 


May 


Jim, 


Jul. 


Aug. 


Sep. 


7.58 

2.88 


7.82 
3, *16 


8.5*1 

2.57 


7.82 

3.31 


7.51 

3.23 


.8.3 1 
3.07 


5.08 

2.81 


0.017 


0.018 


0.20 


0.17 


0.17 


0.19 


0.12 


0.50 

3.73 


7.82 


3.23 


5.90 

3.50 

5.09 

2.88 


6.91 


2.35 


•1.55 

2.44 


0.15 


0.18 


9.77 

4.70 


0.14 


0.13 


0.16 


0.10 


9.59 

3.94 


9.82 

4.46 


10.80 

3.28 


7.65 

3.27 
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An Improved Upper Bound on the Block Coding Error 
Exponent for Binary Input Discrete Memoryless Channels 

R. J. McEliece 

Communications Systems Research Section 

J. K. Omura 

University of California at Los Angeles 


For coded telemetry systems it is important to know the tradeoff between the 
error probability and the complexity of implementation. For systems using block 
codes, the block coding error exponent is a good way to estimate this tradeoff. In 
this article we show how the new upper bounds on the minimum distance of 
binary codes obtained by McEliece et al. result in improved upper bounds on the 
coding error exponents for binary input memoryless channels. 


Consider a binary-input memoryless channel with input 
alphabet A “ {0,1} output alphabet B , and transition 
probabilities {p(i/|x) : xe A,y eB). Let C = {x ( ,x L . ••• ,x, w } 
be a binary code of length n and rate R =■ n~' logs M for 
this channel, and assume that each of the M codewords 
is sent with probability 1/M. Let rf„,j n (C) denote the 
minimum Hamming distance between distinct code- 
words, and let P,(C) denote the probability of maximum- 
likelihood decoder error when the code C is used on the 
given channel. 

Now define 

8(ji,R) ~ —max d,,,n,(C) (1) 

n 

P t ,(n,R) — min P, .(C) (2) 


where the maximum and minimum in (1) and (2) are 
taken over the set of all codes of length n and rate > R. 
And finally define* 


S(R) — lim S(n,R) 
11 « 

(3) 

E(R) = lim -- log 2 P F (n,R). 

(4) 




The maximum minimum distance S(R) is unknown ex- 
cept at the points R — 0, 1 : 8(0) — 1/2, 8(1) 0. The block- 

coding error exponent E(R) is known only for R = 0 and 


‘These limits are not known to exist. However, in what follows 
every upper hound is an upper bound on the corresponding lim 
sup's, and the lower hounds are lower bounds on the lim inf’s, so 
there is no harm in pretending the limits exist. 
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ft > ftcri f> ftcri t being a number to be defined below. We 
shall now briefly survey the known upper and lower 
bounds on E(R), and indicate how the new upper bound 
S*(ft) on S(R) obtained in Ref. 1 can be used to improve 
the known upper bounds on £(R) for small values of R. 

First, the sphere-packing bound E„ P (R ) and the random 
coding bound E r (R), valid for all rates R less than 
channel capacity (Refs. 2, 3): 

Fr(ft) < E(R) < E„ t ,(R), (5) 

The two bounds in (5) are equal for sufficiently large 
ft, and in fact the number R cri , cited above is the point 
where these two bounds meet. (Formulas for E r and 
for binary symmetric and binary erasure channels are 
given in the Appendix.) 

Next, we have bounds which depend on the 
Bhattacharyya parameter (Refs. 4, 5) for the channel, 
which is defined by 

veil 

These bounds are 

aD < E(R) < aS(ft) (6) 

where 0 < D < 1/2 is defined implicitly by R = I - H,(D), 
where //,(%) is the binary entropy function. (The lower 
bound in (6) is called the expurgated bound E ei (R); 
it is only valid for 0 < R < R', where ft' is the rate at 
which the expurgated bound meets the random coding 


bound.) As mentioned, the function S(R) is unknown, so 
the upper bound in (6) is ineffective. However, by 
using the bound S(R) < 8*(R) obtained in Ref. 1 (for 
numerical values of 8*(ft), see Table 1 in Ref. 1), we 
obtain an upper bound 

E(R) < «S*(R) (7) 

which can be evaluated, and which is already better than 
any previously known upper bound for small values of ft. 

Finally, Shannon et al. (Ref. 3) have shown that if 
£ 0 (R) is any upper bound to E(R), then so is the convex 
hull of the curves E 0 (R) and E tp (R). In particular, by 
taking E a (R) = (l/2)a (from (6) and the fact that 
8(0) = 1/2), we see that E(R) is bounce i from above by 
the line passing through the point (0,(l/2)a) which is 
tangent to £ s „(ft). This bound is called the straight-line 
bound E„i(R). However by taking E 0 (R) = a8*(R) 
(cf, (7)) we can obtain an upper bound which is signifi- 
cantly better than min (E sl (R), E xp (R)) for a considerable 
range of ft’s. We illustrate this in Fig. 1 with a binary 
symmetric c hannel w ith crossover probability e = 0.01, 
a = -logo V4<1 - e) = 2.329, and in Fig. 2 with a 
binary erasure channel with erasure probability e = 0.01, 
« — ~ "log" c — 6.644, In both figures die unknown region 
for 0< ft < ft ( ri f in which E(R) lies is shaded. A final 
point worth mentioning is that the new upper bound 
(7) on E(R) always matches the expurgated bound 
£ tv (R) in slope at ft = 0. (Both slopes are — oo; this is 
well known for the expurgated bound, and follows for 
the bound (7) from the results of Ref. 1.) This fact 
supports the conjecture that E(R) — E CI (R) for ft < ft crj , 
for binary input channels. 
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Appendix 

E r (R) and E sp (R) for Binary Symmetric and Binary Erasure Channels 

For a binary symmetric channel with crossover probability e, the random coding exponent is given by 

n - r ~ . (l -I- V4e(l -7) 0 < R < 1 ~ tf=( V«/(V* + Vl “«)) 

^r(R) - _ h,(D ) 1 " H 2 (\/e7(Vi*+ Vl ~ <0) < H < 1 - H 2 (e) 

where T,(D) - -Dlog 2 e - (1 - D)log s (l - «), and D satisfies (7). The sphere packing exponent is 

E ap (R) - r«(D) - H 2 (D) 0 < R < 1 - H 2 (e). 

(Hence R„„ = 1 - HMVtf “ VT^» *< R > = E '< R > = E -< R > for R * R "*"> F ° r the binary ° h<m 

nel with erasure probability e: 

1 1 - R - log, (1 + e) 0 < R < 1 — 2e/(l + e) 

E ' (R) = | E„„(R) 1 - 2e/(l + e) < R < 1 - e. 

where 

E,p(R) ■= (1 S f - lo & (d - £ ) + <2P )' 

where p is determined by R “ 1 “ €2 p /(l — € 4- e2 p ). 
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A Radar Study of the Backup Martian 
Landing Sites 

G. S. Downs, R. R. Green, and P. E. Reichley 

Communications Systems Research Section 


The Goldstone radar system at DSS 14 was used to probe the Martian surface 
at 8495 MHz in a narrow strip between -6 deg and -2 deg latitude. The Viking C 
landing sites lie in this strip, and their altitudes, rms surface slope, and 
reflectivity are presented here. 


The latest Martian opposition (condition of closest 
approach) occurred on December 9, 1975. At that time 
the sub-Earth point on Mars was at a latitude of -1 deg, 
traveling south to a minimum of -6 deg in mid-January, 
1976. Turning northward, the sub-Earth point reached a 
latitude of 24 deg by July 1976 (the planned time of the 
first Viking landing). The backup landing sites, to be 
considered in the event of a failure on the first landing, 
were in a narrow strip just south of the Martian equator. 
The region contained possible landing sites that could be 
probed by radar at a time when the sensitivity of Earth- 
based radars was near a maximum (the sensitivity 
decreases as the fourth power of the distance). 

These backup sites, known as “C” sites, were probed at 
a frequency of 8495 MHz. The R&D X-band transmitter at 
DSS 14, Goldstone, California, was used as a signal source. 
The 64-m antenna at DSS 14 was used for transmission and 
reception. Small cells 25 km E-W by 150 km N-S were 
isolated on the Martian surface. This isolation was 
accomplished by (1) phase modulating a CW signal with a 
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pseudorandom sequence upon transmission, thereby 
facilitating differential range measurements; (2) Fourier 
analyzing the output of each of 32 range j.'ates to produce 
a 64-point spectrum every 65 ms. About 461 spectra were 
combined to produce an average spectrum every 30 
seconds. The basic bit length of the phase modulation was 
8 ps. The maximum doppler shift (due to the rotation of 
Mars) measurable was 984 Hz, so the maximum resolution 
was 15.6 Hz. This corresponds to an E-W resolution of 
0.07 deg in longitude along the doppler equator within 10 
deg of the sub-Earth point. The resolution wa s dropped to 
0.5 deg throughout the data analysis to decrease the 
effects of receiver noise. The N-S resolution, determined 
by the modulation bit length, was about 2.5 deg. These 
resolutions correspond to the cell size quoted above. 

The 30-second data frames collected during a particular 
observing run were combined to yield backscatter 
functions (reflected power versus the angle of incidence) 
for each 0.5-deg region located along the doppler equator. 
Each function was analyzed to determine residual range, 
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the roughness parameter C, and the reflectivity p. The 
parameters C and p are referred to as the scattering 
parameters. The manner in which the cells are isolated on 
the Martian surface and the techniques of determining 
range and the scattering parameters are described in 
detail in Refs. 1 and 2. 

Observation runs were made between 20 December 
1975 and 25 February 1976 for the purpose of studying 
the Viking C sites. The results of 12 runs are condensed 
here by combining results within three latitude bands, 
each 2 deg wide. The bands are centered at latitudes of -6 
deg, -4 deg, and -2 deg. All numerical results at a 
particular longitude, pertaining to a particular parameter 
and located within 1 deg in latitude of the center of the 
strip, were averaged together. The results are presented in 
Figs. 1 and 2. Figure 1 corresponds to the strip 2 deg wide 
N-S centered at -6 deg latitude, while Fig. 2 corresponds 
to the strip centered at -4 deg, with one exception: all 
data in Fig. 2 for longitudes greater than 300 deg apply to 
the strip centered at -2 deg latitude. The residual ranges 
have been converted to altitude in km. The reference 


surface corresponding to 0-km altitude is the ellipsoid of 
Standish (Ref. 3) with a radial scale of 3397.515 km. The 
Martian ephemeris used for the final analysis was obtained 
from the JPL ephemeris DE 96, modified by a small range 
polynomial determined by M. Standish using, in part, the 
observations described herein. The roughness parameter C 
has been converted to B r , the estimate of the rms slope of 
the surface, where B r = C~^ /z radians. The one-standard- 
deviation uncertainty in altitude is less than 0.5 km. The 
plot of each B r and p in Figs. 1 and 2 is a vertical bar. The 
middle of the bar represents the estimate of the 
parameter. The length of the bar represents one standard 
deviation above and below the estimate. (Any bar 
exceeding 8 deg or 12 percent is truncated). Hence there 
is a 68-percent probability that the true value of the 
parameter lies on the bar. 

The positions of the seven C sites are marked on the 
appropriate figures. The site nomenclature and position 
are those in use as of 22 April 1976. In some cases the site 
is actually a region and is marked at two longitudes or 
latitudes. 
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LONGITUDE, deg 


Fig. 2. Martian radar surface properties as seen by radar at 8495 MHz, -4 deg latitude, 
or —2 deg latituue (longitude > 300 deg): (a) altitude, (b) rms slope, and (c) reflectivity. 
Longitudes for the C sites are marked 
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The Effect of Direct Current Bias in the 
Computation of Power Spectra 

E, R. Rodemich 

Communications Systems Research Section 


We determine the effect of dc bias in the approximate computation of spectra 
of Gaussian processes by hard limiting. 


I. Introduction 

The spectral density S(«), |<d| < n, of a stationary 
Gaussian process {**}, — co<k<cc, of mean zero, 
can be expressed in terms of the covariances R x (k) = 
E{x n x nik ) by 

S(o») •= R*(0) 4 2 jp R*(k) cos foo. (1) 

■Vi 

if we put 

> 0 , 

V ‘ ' j~l,x, <0, (2) 


Then it is known (Ref. 1) that R u (k) = satisfies 

the relation 

R,(k) =B,(0)si n r|-n„w]. (3) 

In practice, the series in (1) is truncated, and the 
correlations Ry{k) arc estimated from samples of finite 
size, which leads to random errors in the evaluation of 
S{< u ). These errors will not be considered here. Neglect- 
ing them, wc use the variables {?/(,} to estimate R v (k), 
then apply (3) and (1) to get S(w). This leads to a 
saving in computation time over the direct estimation of 
R r (k) from {x k }, since these estimates require a large 
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number of arithmetic operations. The variance o-” ~ 
Zlr(O) •" E{xl), if needed, must he estimated separately. 
This variance, which enters S(m) as a scale factor, is often 
unimportant, because it can be affected by many extra- 
neous factors. 

We consider the following problem: Suppose that 
an unknown bias a is added to the Gaussian process {x*}, 
so that we get {x k -f (i),-zc < k < oo. The formula 
(2) cannot be applied if a is unknown. However, we 
can take 

( + 1, x,- 4- a > 0, 

1/i • \ , 

I - 1, .Vi -r a <0. (4) 

Then v/e can estimate the numbers E k = E(y„y tuk ) (no 
longer correlations) take 

R ?(k) - sin (jEk), (5) 

in analogy with (3), and form 

S*(») = R* (0) +2 2) R*( k) cos h>. (6) 

k-1 

If a ~ 0, S*( w ) is S(w)/<r-. If a^O, S*(w) approximates 
S(w) in some sense. We shall investigate how close this 
approximation is. 

To find S(®) exactly, (3) needs to be replaced by a 
fov T iula involving a non-elementary function which de- 
pends on u/ ( t. This is not practical for real-time compu- 
tations. Hence we consider the above approximation. 

II. Formulas 

Get the Gaussian random variables x k have variance 
rr-. Then the value of E k = E(y„y n+k ) depends only on 
pk — Rj(k)/ir'-, and can be written as 

= f(Pk)> (7) 


( + I,f > ~a/iT, 

!/W = ) 

This follows if we evaluate E k directly as the expectation 
of a function of x„ and x, by integrating over the 
bivariate Gaussian distribution (Ref. 2), 

If we differentiate (8) with respect to z, the result- 
ing integral can be simplified by integration by parts, 
and we get 

Integrate this equation from 0 to z. Using the variable 
of integration u = sin % 

i /“""K-TTahr)*,- <» 

To evaluate /( 0), note that at z = 1, the distribution 
in (8) is concentrated on the line U ~ t->. Then /(l) — 
E(l/(*i)*) = 1. Hence, putting % — 1 in (9), 

1 ' m = H ' ( - irisnr) rfu - do) 

By numerical integration of tlie integral in (9), going 
from ~ — 0 to z ~ 1 and to z — — 1, tables of f(z) were 
obtained for various values of fl/y. Some of these are 
plotted in Fig. I. 

If we plot [sin (*/2)/(s) j instead of f(z), the curve for 
a ~0 becomes a straight line, and the other curves are 
also modified, but they do not come significantly closer 
to the line for a — 0. The separation between the curves 
is the error in this method of approximating R T (k)/<P, so 
at first glance this method does not look too promising. 

However, instead of R*(k), we should be considering 

__ smfj £*] - sin f"|-/(0)l 

R,(k) = L- J L? — J 


* [* r u \ n N f t\ 4- tl - 2zt t t, 1 , 7 

j y J r y(E)y(tj) exp UU x Jt* 



7T _ 


sin 

j E "J _ sin 


sin 

J “ sin 

f/(0) 


- C, [R*(k) - C,], 
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where 


and 



0, — (i — c 3 )-«. 

R-Ak) is plotted as a function of Ii^k)/^ in Fig. 2. We 
see that R^(k) is a good approximation to R t (k)/<r if 
RAk)U is not too close to -1 (which is usually true). 


Forming the numbers R y (k) requires knowing some- 
thing about the function f(z). However, they are only 
wanted to construct the function 


S(m) — 1 -J- 2 2 R y (k) cos ka>. (12) 

ks i 


The error R y (k) -- Rj(k)/tr is negative for R x (k) > 0, 
with a minimum near R y (k) = 0.44^ for all values of a/ a 
between 0 and 1, For R r (k) < 0, the error is positive, 
inci easing as R y (k)/ <t~ goes from 0 to — 1, These errors 
are given at R#)/V - 0.44 ancl -0,30 for several values 
of a / it in Table 1. 


Another method which we can consider for comparison 
is to first compute the correlation 


then take 


R,,(k) 


^ - m 

1 - /( 0 ) ’ 


and 



This function can be found more directly as 

S(«)=C 1 [S*(-)~C s 8(«)]. (13) 

The piesence of a dc component in (j/ t ) causes a spike 
in the spectrum S*(<n) at w ~ 0. 1 his can be removed by 
inspection, which is the significance of the term — C..8(«) 
in (13). Then the multiplier C l can be adjusted to give 
whatever scale is desired for the spectrum. Thus we get 
S(<»), defined by (12), without any particular knowledge 
about a or the function f(z) 


Sfw)- — 1 2 ]jP Rj-(k) cos k<». 

k l 

Tins method is suggested by the original relation (3). 
Curves for R„(fc) as a function of R y (k)U- y similar to the 
curves of Fig. 2, can be plotted, These curves have the 
same general appearance, but they range farther from 
the line for a = 0. Values of the error /?,(£) - R r (k)/,r" 
are given in Table 1 for comparison. Wc see that this 
method is not as good, as well as being harder to 
implement. 
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Table 1. The errors in R^k) and RJk) for R x (k)/a* ■ 0.44 
(upper number) and -0.30 (lower number) 


(l/o 

0.2 

0.5 

1.0 

jyfc) - RjLky/o* 

-0.0001 
+ 0.0003 

-0.0031 
+ 0.0093 

-0.0342 

+0.0822 

fl>) - RM/** 

-0.0031 

+0.0064 

-0.0192 
+ 0.0378 

-0.0738 

+0.1261 
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Standard Random Number Generation for MBASIC 

R. C. Tausworthe 

DSN Data Systems Section 




I. introduction 

The first MBASIC random number generator (Ref, I) 
implemented on the Univac 1108, used a linear congru- 
ential method, *„ +1 = 5 1! %, (mod 2 as ), followed by nor- 

mobibhM ^ the ianS<? ThiS Senerator was used 
probabJy because it was already available in the U1108 

s atistics package. Empirical tests by users, however 
later proved that the generator possesses very nonrandom 
correlation properties indeed, especially unless great care 
is taken in specifying the initial “start” value. 

(Jt%ZTfu n def ? Ult S ' arting value ’ or “seed,” 

L" f 9 , \ iad J®? 11 s f ec ; ted because theoretical results 
(Ref. _) claimed that the longest nonrepeating sequences 
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are obtained when r„ is an odd power of S. Unfor 
nately, the particular x„ chosen, while perhaps generati 
a long sequence, does not seem to produce numbers w 
vei-y good randomness properties, 

Tins article describes an alternate generator of a tv 
'hose randomness has been theoretically shown to' 1 
vastly superior and which can be implemented on ai 
computer, despite word length restrictions (the UlU 
algorithm was tailored to 36-bit words). It is a machin 
independent algorithm. 

Tho generation method is almost as fast as the line- 
meth0d ’ bUt " 0t qUita The rati0 of speec 
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II. The Linear Recurrence Modulo 2 Method 

In 1965, the author (Kef. 3) showed that numbers 
produced as successive binary words of length s taken L 
bits apart {s < L) from a linear (shift -register) bit-stream 
recursion of the form 

a„ — 4 -f c v . 4 a,,.,, (modulo 2) 

form a pseudorandom sequence whenever the polynomial 
f(x) — x p 4- 4 ••• 4- C\X 4- 1 is primitive over 

GF( 2). The algebraic structure of these pseudorandom 
numbers provided a way of proving that, over randomly 
chosen starting values in the numeric sequence, 

the correlation between numbers is essentially zero, 
being ostensibly equal to — 2~’\ for all numbers in the 
sequence separated by less than (2 P — s — 1 )/L, subject 
to the restriction (L,2'‘ — 1) — 1, The author further 
showed that adjacent /c-tuples of such numbers were 
uniformly distributed for 1 < k < ( p/L ). 

The algorithm for computing the numbers is simple, 
especially when f(x) is a primitive trinomial, say 
x p 4 x Q -f 1, where q < p/2. An even greater simplifica- 
tion is possible, as is shown in this article, when p is an 
even multiple of L, as is the case for the trinomial 
t r,;u _l Y :<; + ^ (f rom Zierler and Brillhard (Ref. 4). The 
generator based on this polynomial has period 1.4 X lO 1 '*", 
has virtually no (average) correlation between any num- 
bers separated by less than 5 X 10 1SS , has 28-bit precision 
numbers available, and has adjacencies up to 19 dimen- 
sions uniformly distributed. Runs-up and -down statistics 
up to length 16 are impeccable. The period and maxi- 
mum correlation distance are, in fact, so great that the 
generator would have to produce numbers at a nano- 
second rate for more than 10“- years before nonrandom 
distribution or correlation effects would be noticeable 
as nonrandom. Almost 4 X 10" numbers would have to 
be examined to detect deviations in runs-up and runs- 
down statistics as nonrandom. 

These pseudorandom number generators have been 
widely studied (Refs. 5—5) since 1965, both theoretically 
and empirically, and have been “promoted to pride of 
place in the field of pseudorandom number generation 
(Ref. 7).” 

Tootill (Ref. 8) has even discovered generators of this 
type for which “there can exist no purely empirical tests 
of the sequence as it stands capable of distinguishing 
between it and [truly random sequences].” For reasons 
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having to do with computer storage and precision, the 
generator of this article is, unfortunately, not one of 
these. Nevertheless, the generator described is vastly 
superior to any linear congruential generator in existence. 

III. The Method 

The algorithm for producing the succeeding p bits 
from the current set of p bits in the x 1 ' 4 x' 1 4 1 generator 
is, 

(1) Left-shift the p bit string by q bits, inserting q 
zeros on the right, dropping q bits on the left. 

(2) Add modulo 2 (exclusive-or) the original and 
shifted p bits. 

(3) Right-shift this result by p — q bits, supplying 
p — q zeros on the left, and dropping p — q bits 
on the right. 

(4) Add the results of 2 and (3) above modulo 2 to 
give the next p bits. 

The proof that this algorithm works is very simple, and 
the reader is invited to apply the algorithm to the bit 
string a, a- • • ■ a,, and use the reduction formula 
(ip* m r a,» (try it with aia 2 a x aia r , with — 

J+W f t/ffl to see what is happening). 

Note in the method that the number of computations 
required for generating the next p bits grows at most 
linearly in p. Assuming p > > L, then partitioning the p 
generator bits into words of length L produces a number 
of words that also increases in proportion to p. Therefore, 
the number of computations for L bit precision random 
numbers, to first-order effects, is independent of the 
recursion degree p. Making p large, however, has advan- 
tages in increasing randomness properties. 

It is true that as p increases, more and more registers 
are required to compute each new set of p bits, and 
shifting many registers at once presents a small incon- 
venience in most computer languages. These factors 
place a small speed and storage overhead on the genera 
tion process; but as we shall see, even this is not extreme 
due to the particular trinomial chosen. 

The trinomial .v"" 1 - 4 Jt i7 4- 1 (p = 532, q = 37) has sev- 
eral things to recommend it: (1) 532 is factorable into 
28 * 19, which means that 19 words each having 28 hits 
precision can be generated all at once by the algorithm; 
(2) up to and including 19-tuples of adjacently produced 
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random numbers will be uniformly distributed; (3) 28 
and 19 are botli relatively prime to the period — 1), 
so no ill effects occur as a result of beginning new words 
every 28 bits; (4) the period and correlation distance are 
so great as never to be witnessed in the lifetime of the 
universe; and {5) the q — 37 value produces good runs-up 
and runs -down statistics (Ref, 8), up to runs of length 16. 

On converting 28-bit fixed-point mantissas to real num- 
bers, S-digit precision results. Some machines, such as 
the 1108, may have to reduce this precision in order to 
fit the floating-point exponent field into the word (the 
l '1108 has only 27-bit mantissa precision). In such cases, 
the most significant bits of the generated words may 
always be retained without losing randomness qualities, 
so that all implementations produce essentially identical 
results, within machine precision. This philosophy is 
present in the algorithm which follows in Section IV. 

The RANDOMIZE function which initializes the gen- 
erator uses a multiplicative linear eongruential algorithm 
to generate the first 19-number “seed,” from which the 
rest of the generated numbers grow. The particular value 
for the multiplicator a in the algorithm 

*n-i £=■ ax n (modulo 2 1 ') 

was chosen as 41,475,557 10 for the L — 28 case from 
theoretical results published in Ahrens and Dieter (Ref. 9), 
who show that good randomness results when 

a = 5 (modulo 8) 

where | = (5 1/1! — l)/2 (the “golden section number”). 

IV. The Algorithm 

Managing the 532-bit shift-register is the main trick in 
implementing the method of the previous section. The 


algorithm given in this section utilizes an array Wj, 
i — 0, • * • ,18 of 6-bit computer words (b > 28) sufficient 
to encompass the p — 532 span of bits to be operated on 
(and retained), and delivered in 28-bit chunks whenever 
the RANDOM function is invoked. 

On machines having word sizes smaller than 28 bits, 
double words will have to be used for each word in the 
algorithm below. Two values of the RANDOMIZE 
argument that round to the same internal fixed-point 
representation will produce the same random sequence; 
conversely, every unique fixed-point representation of 
the argument produces a unique random sequence. In 
addition, so long as the value of the argument is the 
same and stays within the precisions of two differing 
machines, the sequences produced on each will be the 
same, within machine precision. 

Counting the number of elemental operations (load, 
store, shift, etc.) for the algorithm, one finds about 10 4- / 
operations per number generated, where f is the number 
of operations in “floating” the fixed-point number. The 
linear eongruential algorithm requires only about 3 + / 
such operations, so the ratio of speeds is approximately 
3:1. 

Tlic RANDOM function is about 23+/ operations 
long, as compared to 3 + f for its linear eongruential 
form, and data storage is 21 words vs 2. However, even 
though the program requires perhaps 7 times as much 
core as the linear eongruential form, the total is still 
probably u ider 50 locations, of negligible concern in 
most installations. The asymptotically random sequence 
of Tootill (Ref. 8) requires 607 words to store the array 
W, alone. (This, coupled with the fact that only 23-bit 
precision was available in that generator, is why it was 
not considered here.) 

The two functions, described in CRISP-PDL syntax 
(Ref. 10), are as follows: 


Function: RANDOMlZE(starter:universai real) 

< * This function declares and initializes a 19-clement 

< * array, W,„ . . . ,W with random numbers generated 

< * by a linear eongruential method. An integer I is 

< * set to zero to enable RANDOM to select YV 0 as the 

< * first random number. W and 1 are permanent data 

< * structures, accessed only by RANDOMIZE and RANDOM 
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.1 constant multiplier; integer — 41475557 

.2 variable /:intcger, l\ integer, 

W : array [0 . . 18] of universal integer 
,3 if (starter <0) < * MBASIC manual specifications * > 

4 : starter: clocktime < * returns current time of day as integer* > 

.5 : - - >’(starter=0) 

.6 : starter; = multiplier 

>(else) < *starter>0* > 

7 : starter: = fix(starter) < * floating-to-integer conversion * > 

: , . endif 

.8 W„: _r starter 

.9 loop for /-■ 1 to 18 

10 | W; - (Wj ,* multiplier) modulo 2* 28 

.1.1 *— <— repeat with next j 

.12 Z 0 < * set up to pick W„ as first random number * 

endfunction 

Function: RAX DOM : real 

< * This algorithm makes use of a 19-word array, 

< * W,„ . . . ,W, s, each with b> 28 bits. Each word contains 

< * 28 hits of the generator, right justified. An 

< ’ integer variable 1 on entry contains the index of 

< the word array next to be returned as the random value. 

< « Roth \V and Z are permanent data structures, accessed 

< « only by RANDOM and RANDOMIZE. The latter of these 

< initializes / to zero and W to the seed. 

variable /: integer 

if (/ 19) < * all words have been used up * > 

: /; — 0 < * reset to first element in array * > 

: loop for / • 0 to 16 < * exclude last two words * > 

: | load W (+I ,\F;^ into A,„A| 

: 7 left shift A, by — 28 < * join bits in stream * > 

: 7 left shift A,„A, by 9 < * q ^ 37 is 28 i 9 * > 

: 7 \V r . (W, NOR A„) sion formula * > 

: <-«- repeat with next / 

: load W,„W 0 into A ll7 A, < * now compute W I7 : * > 

: left shift A, by b- 28 < * join \V,„ \F„ bit streams * > 

: loft shift A,„A, by 9 < * A„ now has final 19 bits of W,e> 

: U 7 ,-: --\V i7 XOR A,, < * and first 9 bits of stream shifted 495 * > 

: load into A„,A, < * do similarly for * > 

: left shift Aj by b 28 

: left shift A„, A, by 9 

: \YV: 1 - (W,* XOR A„) 

: . . endif 

RANDOM float(W,)/2 * < 28 < convert to real * > 

Z: Z i 1 

endfunction 


.1 

.2 

.3 

.4 

.5 

.6 

.7 

.8 

.9 

.10 

.11 

.12 

.13 

.14 

.15 

.16 

.17 

.18 

.19 
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V. Conclusion 

The algorithm given has been implemented as the 
RANDOM function in the MBASIC processor currently 
on the Caltech PDP-10 computer. AH tests so far run on 
it validates the randomness properties claimed by the 
theory. In that theory, by the way, the only factor left to 
chance is the specification of the initial “seed.” The stated 
uniformity, zero-correlation, and runs statistics are all 
based on the single assumption that the seed be chosen 


randomly. Of course, the default value canned in was 
not randomly chosen, but chosen specifically to look 
random except for the first word, and certainly, to the ex- 
tent of our tests, this appears to have worked beautifully. 

We have demonstrated that the generator is also capa- 
ble (as is every known random number generator) of 
producing numbers with Sit variations from random- 
ness over a few thousand samples when the wrong seed 
is supplied. 
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It is shown that Reed-Solomon (RS ) codes can he decoded by using a fast 
Fourier transform (FFT) algorithm over finite fields GF(F„), where F„ is a Fermat 
prime, and continued fractions. This new transform decoding method is simpler 
than the standard method for RS codes. The computing lime of this new decoding 
algorithm in software can be faster than the standard decoding method f for 
RS codes. 


I. Introduction 

Recently, Gore (Ref. 1) proposed the usage of a Finite 
field transform over GF(f/"), where q is a prime number 
and n is an integer, for decoding RS codes. Michelson 
(Ref. 2) has implemented Mandelbaum’s algorithm 
(Ref. 3) and showed that the decoder, using a transform 
over GF(q “ ), is faster than a more standard decoder 
(Ref. 4). The disadvantage of this transform method is 
that the code length is such that the most efficient 
FFT algorithms cannot be used to yield fast transform 
decoders. 

Rader (Ref. 5) proposed transforms over rings of inte- 
gers modulo Fermat numbers. Such transforms can be 
used to compute error-free convolutions of real integer 
sequences. Agarwal and Burins (Ref. 6) extended Rader’s 


Fermat number theoretic transform by using the gen- 
erator u -- \/2 for the transform rather than a — 2. If 
\/2 is the generator of the transform, the transform has an 
FFT algorithm which can be used to calculate transforms 
with as many as 2"‘- points of integer data. This trans- 
form was extended to residue classes of quadratic inte- 
gers Ir n (^ 2), where \/2 is a root of x~ — 2 r - 0, F„ is a 
Fermat number and I t (| denotes the set of integers mod 
F„ (Ref. 7). McClellan (Ref. 8) lias realized recently the 
hardware for the Fermat prime theoretic transforms. He 
showed that the arithmetic used to perform these trans- 
forms required only integer additions and circular shifts. 

Recently, Justcscn (Ref. 9) proposed that transforms 
over GF(F n ), where F„ - 2 2n 4- 1 for n - 1 ,2,3,4 is a 
Fermat prime, can be used to define RS codes and to 
improve the decoding t iency of these codes. 
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Recently, tlu* authors (Ref. 10) extended the transform 
to the finite field of type /, (| (\ s/ 2) (isomorphic to GF(F„)), 
w here a/ 2 is a root of the polynomial P(x) = x* - 2 over 
GF(b„), and l r>l denotes the set of integers modulo F„, 
Again the arithmetic used to perform this transform 
requires only integer additions, circular shifts, and a 
minimum number of integer multiplications by powers 
of \/2. An FFT over the finite field of type hjV 2) can 
be list’d to encode and decode RS codes of as many as 

1 symbols for n • 3,4. Encoding and decoding can be 
accomplished faster and more simply than any other 
know’n standard decoder for RS codes of the same 
symbol range. It was also shown (Ref. 10) that the FFT 
over G/'(K* 2" I 1), where K and n are integers, can be 
used to encode and decode a class of RS codes. A special 
case of the radix — S FFT over GF(q -), where r/ - 2" -- 1 
is a Mersen ue prime, was developed to encode and 
decode another class of RS codes. 

I he decoding of systematic ltecd— Solomon codes using 
the translonu over GF(F„) was composed of the follow- 
ing three steps (Ref. 10). 

(1) Compute the Ft V over GF(F„) of the received 
code .Y-tuple; i.o„ 


Theorem 1: Let q be a prime number. Also let A(x) be 
the formal power series of form 


S‘x) (1) 

i=0 

w'hcre a t eGF(q), the degree of A(x),cl, is a real integer, 
and ^ 0. Define the set 


V i It 


and rf in nn I 


n t 


such that addition is given by 


G(x) = B(x) + C(x) = £ hxr-' + £ c,x f >, e > f 

i -0 i o 

= X " >xe ~' 

i-a 


where 


g, 


j bu i <e - f 
I b, -I Cj n J >e — f 


S 


K 


E 


and multiplication is given by 


wlu're y ]u eGF(F„) and u is an element of order N. 

(2) I’se Berlekamps iterative algorithm (Ref. 11) to 
detf'rmine it, from the known Sj = E t for i -- 
I--- 1 " 4 and / --- 1,2, ■ ■ • ,2f. Then compute the 
remaining transform error E r 

(3) Compute the inverse of the transform over GF(F„) 
°f S\ E s to obtain the corrected code, 

An advantage of this transform decoding algorithm 
over other methods is that a FFT over GF(F„) can be 
used to compute the syndromes and error magnitudes, 
In this paper, Berlekamp’s iterative algorithm can be 
modified by using continued fractions in GF(F„). This 
modified Berlekamp’s algorithm can be easily imple- 
mented on a digital computer. 


H(x) = B(x) ■ C(x) = ( £ bx -^ (£ o,x’ A 

= X 

r u 

where 

lh ■ t, b;C, , 

) 0 

Then R is an infinite field. 

Proof: It is evident that R satisfies the postulates of a 
commutative ring with unity element. An additive iden- 
tity element and a multiplicative unit element in this 
ring are 


II. New Approach to Decode Reed-Solomon 
Code Using the Transform Over GF(F„) 

In this section, a new' approach is developed to define 
and decode RS Codes. The following theorem and defini- 
tions are needed. 


and 


D(x) = £ di - 0 for i = 0,1,2, 


F(x) 


0,/n = l,fj ~ 0 forj > 0 
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Every nonzero element of F, i.e., 


A(*) =£«•**'* 

I ~i\ 


that the field of type 7, H (\V2) is isomorphic to GF(F„) 
for n = 3,4 and that a - ^f%eGF(F») is an clement of 
order 2"M In these fields a systematic RS code can be 
specified in GF(F») as follows: 


has an inverse element B{x ) defined by 


B(x) = A(x)~' 




I -II 


Assume the code length for tlie RS code is A 2” 4 . 
Let a code word be represented by /(*), a polynomial 
of degree N 1 over GF(F„). The generator polynomial 
of g(.v) is defined as 


win ‘re 


g(*) 


f- 1 

Z> 


i i 


a*) 


K = 

b , = for (* = x > 2 . ‘ ‘ * ) 


where 


d = 2 K < A 7 — 2 ,l+l 


Hence F is a field. 

If the set R is composed of all power series not con- 
taining negative powers of x, i.e., 

R — i ' (l j | (i jCfJ/Vi i , il i mi jirtsl 1 1 1«> i lit oi-’t’r / 

(2) 

then it is evident that R is a subring of F . The integer 
part [A(.v)] of A(x) in (1) is defined by 

[A(x)} 


and a - i/ 2 ,a- - (>/ 2 )‘, ••• ,a rf = (^ 2 ) l/ are the roots of 
g(.v) in GF(F„). The resultant RS code* with A 7 symbols, 
which is a multiple of the generator polynomial, is com- 
posed of d - 1 parity check symbols and iV - (d 1) 
information symbol 0 , where d is the minimum distance 
of the RS code. If Ms the number of errors, the code will 
correct, then d 2 / 4 1 . 

Suppose that the code word f(x) = f„ I- /,* + • • • 
-f- /v-iA , v ' ! is transmitted over a noisy channel. The re- 
ceived code word R(x) “ yo + 71 * 4- • • • f y.v-i *'" 1 is 
composed of the original code word with the addition of 
possible errors, i.e., 

y(--v) - f(x) 4 c(x) 


Let p{x) be the ratio of elements in R, i.e., a rational 
fraction form of type, 

t 

p(x) - 

L'W-' 

i 0 

where 

Z °i- v ' £ ' 1 ' X) b t x f ~'eRcF 

i r| i cj 

By theorem 1, it can be proved that p{x) is an element in 
the field F. An element in F of form p(x) is called a 
rational element. 

Let GF(F„) be the finite field, where F„ - 2- n + 1 is a 
Fermat prime for n = 1,2, 3, 4. It was shown in Ref. 10 


where c(x ) — e 0 -i- e x x I e»x- -r ••• I e.v-i-v '' 1 is an error 
polynomial. 

Upon receiving the message y(.v), the first step in the 
decoding process is to take the FFT of the message in 
lr n (V S'). The transform is taken of the received N -tuple 
message (y Ml y lt • • ■ the coefficients of the polynomial 

y(x). This transform is 

s K =- Zy^) hn 

ii n 

it .0 

--- 

M-U 'I " 

- F k 4- E* for K “ 0,1, - • • ,N 1 
Since f(x) is a multiple of g(-t), /(a*) = 0 i • 1,2, ■ ■ ■ ,d ~ 1- 
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Hence, 


Multiplying the above equation by YjX', on*, gets. 


S» E, - 

r; - E ^[(^V]* for fe = 1,2, •■•//- 1 


( 3 ) 


where E K €GF(F„) is periodic with A\ Let Y* and X, be 
the /th error magnitudes and the ith error location, re- 
spectively. Then the syndrome in (3) becomes 


YiXN - a.Xf'-- 1 -1- <r 2 Y,X'" - 4- ••• + (-l)V,Y f XJ 


Summing on i for i — 1,2, •••,/, then 




(-iv<r,Ev,x{ = o 

i -1 


Using (4), wo have. 


s* - E, «£V.X,* fori- - 1,2, ,rf - 1 (4) 

i -1 

The generating function of the sequence ( E K ) is defined 
as a formal power series. That is. 


Sj +! — <nS/ +f -i 4- ••• + ( — l)VfSj ~ 0, for/ < t 

and 

E)< f ■■■ ViEj+t , 'h •“ 4- ( — l)'(r<Ej = 0 for /> t 


E(x) = Ei-t" 1 + £ 2 .r- 4 E r 1 
where 

Ea-cGFCF,,) 

Substituting (4) into (5), one gets 


E Ea-*" a ' 1 5) 


Thus 


LTv)----:££y,X?.r-‘ 

A i r t 

-Ev.f (x.*-r 

i ■■■ : a i 

> , Xi.r -1 

" i -x,*-' 


£(x) — E,*- 1 4- E 2 .r- -r £ t x- 3 + 


,L .. X 


^ m_ 

X — X, <r(x) 


( 6 ) 


where 


F(x) and o-(x) — IT (* — X ) efl 

i 1 


( 7 ) 

It will be shown in the next section that <t(x) in (6) can 
be calculated by using continued fraction approximations 
when only the first 2 1 coefficients of E(x), i,e., S,, 

S,m are known. If the coefficients of <r(x), i.e,, rr, for 
i — 1,2, -,t, are known, then Eq. (7) is used to obtain 

E„, E,/, E, m, •■*, E.v-i, and the transfonn of the N-tupIe 
error pattern i.e., (E,„ £,, JE a , • ■ •, E.v-i) is obtained. 
Thus, the A' -tuple eiTor pattern (<?„, e,, ■••,c v . 1 ) is found 
by taking the inverse transform over I F (\ s,/ 2) of E A - for 
/c — 0, 1, • ■ •, N — ■ 1. Finally, the original N-tuple code 
word can be computed by subtracting <?„ from the re- 
ceived code word r„. 


III. Implementing Berlekamp’s Algorithm by 
Using Continued Fraction Approximations 

It was shown in the previous section that E(x) = F(x)/ 
<r(x) in (6) is a rational element in the field of all formal 
power series F. Thus, using a procedure precisely similar 
to that used for rational elements in the real number 
field, described in Appendix A, it is possible to use con- 
tinued fractions to develop a finite sequence of rational 
approximations to E(x). That is, the recursive formula on 
convergents is given by. 


Note that E(x) is a rational element in F. Since 

fj(x) - x‘ — a-iX* -1 4- (t 2 x 4- ••• 4- ( — l)*o-/, then 

rr(Xj) = 0 - X{ ~ .r,XH + 

•}■ • • • + ( — l)*tff, for i “ 1,2, ’ • t 


E T (X) 


where 


qs(x)Ps-,(x) - P,s-- 2 (x) _ P s (x) 
( i *(x)<Ts-i(x) 4- oy?_ 2 (x) “ <rs(x) 


( 8 ) 


P } (x) - q t (x),P 0 {x) = l,t7 0 (x) = 0, cti(x) = 1. 
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The partial quotients q, N (x) in (8) can be computed by the 
following formula recursively, 

ft, -.,(*) = qs(x)R,. { {x) + R.s(x) (9) 

where ft-,(.r) — E{.v), R„(x) = 1 and q,y(:v) is obtained as 
the “integer part” of R„^(x)/Rs m(.y) and JL(.v) is the 
“remainder”; and 

( - 1 YRs(x) = a«(*)E(x) - Ps(-r) (10) 

where a,-(x) and P N (.v) also satisfy the recursion in Eq. (9) 
with the initial values given in Eq. (8). By applying 
Euclid’s algorithm to the rational element E(x) in F, 
observe that E(x) - P s (x)/v s (x) will be terminated when 
Rs(x) - 0. 


It follows from the norm properties that 

IM*)|| = |!c/n(x)o'.s.- 1 (x) -I- <Ts->(x)\\ 

= M - !! c/.(a-)|i (11) 


But, by (9), 

- vmx) - Rs.,(x) ii 
- 1! IW*)!! 

,i / mi irlW*)!| ^ , 


Thus, (11) becomes 


The norm of A(x)eF will be used below in the proofs of 
theorem 2. This norm is defined as follows: 

Definition 

The norm of A(x) — ^2 ii A( a*) b, is defined by 

||A(x)!! = 2 rf 

where d is a degree of A(x) 

The properties of norm ||A(x)!‘ are 

a) AB = :'A • B|l 

b) f A|| >0 and ■ A !! ~ 0 ifA = 0 


■ /r,s'(.v) I! > for 1 < s < n. 


To compute the norm of the difference E(x) — P s (x)/ <r.s(x), 
we observe that 


i E(x) - 


Ps(x) 


, "*(*) i 

Then, by lemma 1, 


||S(*M*) - Wll _ "Mx)\\ 
M*)ll ^cr»(«)H 




w 


<*) 


< 


E(x) - 


Ps.*(x) 

<ry.,(.x) 


for 0 < s < n — 1 

(12) 



d) || A = ft! 1 < max(;|A||,!!B : : ) if !'A|| - |IB!! andHA ± /3H 
— max( M AI:, B ) if |r A jj ^' ] B\\ 

Lemma 1: Let n be the smallest finite integer such that 
R n (x) — 0, where R,(x) is defined in (10); i.e., E(x) — 
P„(x)/(T n (x). Then |]R>-(x)|j is a monotone decreasing 
sequence for s = 0,1,2, n and |jo\s-(*)| is a monotone 
increasing sequence for 5 = 0,1,2, ■ n. 

Proof: 

By (9), one gets, 

W) = <Mx)R*-i(x) + Rs(x) 

where deg ft,(.t) < Ky_,(.v). This implies |jRs(x)|[ < jjR,--,(-v)|j 
for n = 1,2, -",17. Furthennorc, since tri(x) = 1 and <r„(x) 
= 0, then |M*)|l > ]|cr 0 (a:) ||. Assume ||ar s _,(x)|| > |]«N- 2 (x*)j| 
for all ^ < n. By (8), 


For decoding RS codes over GF f (F„), we only know the 
first 2 1 coefficients of E(x) in (6). That is, 

E(x) = E, .v 1 + E,.r- + • ■ • E^xr- 1 + X.y--'- 1 + • • ■ 

X(*) 

where X(.v) is an unknown element in F. The following 
theorem is developed to recover the rational element E(x) 
in F precisely when only the first 2f coefficients of E(x) in 
(6) are known. 

Theorem 2: Let E(x) = P(x)/a(x) in (6) be a rational 
element in F defined by theorem 1, where P(x) and 
o-(.v)eft are defined in (2) and j|E(x)|| < 1. Let X(x) be 
an unknown element in F such that deg X(x) < — 2 deg 
u(x). If the first deg X(x) + 1 coefficients of E(x) are 
known, i.e., 

E(x) + X(x) = E,x’ + E»x~~ + ••• + E l! ,«.v«x,M^ K - Y(r,+1 


Gs{x) — q s (x)<rx- } (x) -I- ir». a (x) 


+ Xx lh '* xu) + ■ 
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where 

X(x) ~ Xx' h ‘* XiT) 4 ••• 

1 lion E(x) cun be obtained by using the continued frac- 
tion algorithm operating on E(x) + X(x). 

Proof: By (10), we know that 

(-inUx) - cr«(x)fl2(.r) 4 X(.v)] -- P,(.r) 

“ ((j.s(.v)E(.v) — P.s(.v)) 4- (t.s-(.y)X(.v) 

where <r. s (x)A r (.r) indicates the location of the unknown 
coefficients in R s (x), 


Thus, from (15), (17), (18), we have 

!|<'n i(*) 11 ~ t| Rs---(x) ||-' (19) 

Substituting (19) into (14), one obtains 
!! Rs-i(x) |[- > |l X’(.r) I 1 

Hence q4 x )> which is independent of X(x), is obtained 
by the Euclidean division of R s -M by Rs~\(x) if and 
only if 

ii Rs-t(x) !l“ > |i X(.v) !' (20) 


We see that following the Euclidean division of R s Jx) 
by R s - i(.v), it follows immediately that ( j s (x ), which is 
independent of A’(.v), can be detenu ined if and only if 


deg Rs-i(x) - deg(rs- 1 (x)X(.v) > deg R s ,(x) 
- deg Rs-i(x) -- deg q„(x) forS > 2 


(13) 


(Note that the left side of (13) indicates the number 
of known coefficients in the divisor and the right side of 
(13) indicates the degree of the partial quotient f/«(.v).) 
It follows from (13) that 

i (.r) J ( J-I . «| ( <K iJn ,u-|A(Xt 

By the properties of noi*m, (13) becomes 

[I Rs-i(x) |[ 2 > I; 2 (.v )" ; <r,-,(.Y) ! ' jjX(:e)|; (14) 

But, by (11) in the proof of the lemma 1, one has, 

!!^,(x)!| (15) 

Since, by (9), 

II R«- 3 (x) II - !! q4x)R,^(x) + R,(x ) ; (16) 

Then, by the lemma 1, (16) becomes 

I! «*-=(*) !! = !l q*(x)i\.,(x) !! - 1! q s (x) |J • j| Ry-,(x) || 


From (8), one gets 


and 


li ff i(x) !l- 1 
If Ho(ar) I! = 1 


( 17 ) 


(18) 


Let q„(x) be the last partial quotients such that q»(x), 
which is independent of X(.y), can be determined by the 
Euclidean division algorithm. It follows from (20) that 


WR,,(x)^< 'X(.y) 1 1 
By (19), (21) becomes 

|! O'h-i(-Y) X(x) 

This implies 

||x(yL|-> < i;<r, l+ i(.v) [|- 

Since J| w(.y) |j" < j|X(.*) j r\ thus (22) becomes, 

M*) ir < K + i(*) !r 


Consider either 


E(x) 


PM 


or 


E(x) 


<r>t(x) 


PM 


£(.v) - ^ 

r(x) 


r »{x) 


> 


!fc»- 


m 


v(.V) 


If 


then 


,c,v) -0 
|| 1 1 a(x) ; 


_ P(x) _ PM 

n ( x ) <r„(x) ' 


( 21 ) 


( 22 ) 


(23) 
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If 


E(x) 


PM 

<r„(x) 


> 


E(x) - 


m 

tr{x) 


= 0, 


then 


Ef*) - 


P.{*) 




-li(w-feil) 


< max 


/ P...W P..M\| 

V«.i(i) <'«( r ) /I 

/Hr-' P "* ,( - r) !l || P ** l(5t) P "^|h («> 


By (12), (24) becomes either 


, v F„(.v) 


1 Eiin(x) 





or 


If 


i v ?«(*) I 


i ?..,(*) p..w 

l £(l) - »„(*) | 

< 

, a, ,(-V) 

, P„(.f) 1 


| r . , P...W 

- „„(*> ! 

i j 



this implies R„(x) — ft«ri(*) — 0. If 


E(x)~ 


W 

"«( x f 


< 


Pntl(x) P'i(x) 


ir„+i(x) <r»( x ) 

\P n . t (x)an(x) ~ <r^{x)P n (x) || 
I] U n+ i(x) |] • l| <T„(x) || 


which implies 


Then 


(25) 


By the same procedure used in the derivation of (A-10) 
in Appendix A, (25) becomes 


P(x) P m (x) 


|| ?(*)*»(*) ” P„(x)tr(x) || = 0 

P{x) _ Pn(x) 


< 


i! !! Ii li 


tr(x’) (t, X (x) 

Multiplying (26) by || <r(x) || ||o-*(ac) || gives 

[I cr(^) II 

I! P{x)(T„(x) - P n {x)<r(x) II = || ffji+ , (£)"][ 

By (23), tills yields, 

SI P(x)tr„(x) - Pn(x)<r(x) || < 1, 
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<r(.\‘) ~ <r n (x) 

Hence the theorem is proved. 

A simple example of theorem 2 for decoding a RS code 
in GF(F„) is now presented. 

Example: Let GF(2 L ’ 2 + 1) be the field of integers modulo 
the Fermat prime F, “ 17. We co .isider a 2-error correct- 
ing 8-tuple RS in GF(17). (Note that this example is the 
same example in Ref. 10.) 

Assume the information symbols are l,2,3,2eGF(17); 
i.e, l(x) = l.r + 2x° + 3r’ + 2x\ By the example in 
(Ref. 10), the encoding of /(*) is the polynomial 

b(x) = lr -I- 2x" + 3r* t- 2.x‘ + 15x 3 + 12x- 4- 2x* + 5 

Suppose that two errors occur in the received words, i.e., 

y ( x ) s= 5 + 2.v + 9x- -4 15.V 3 + 2.r‘ + lx fi + Sx* + lx' 

By the example in (Ref. 10), the syndrome can be calcu- 
lated, using a FFT over G/'(F«). That is, 

8-1 


S, = E„-=i;y„2"* for k = 1, 2,3,4 


Hence, 


S, - E, = -8 
S, = E, = -5 
S : , - E, = 11 
S 4 = E, = -1 


( 26 ) By (6) 


P(x) 

E(x) = — 8x _1 - 5.x- 2 + Hx- 3 - Jr 4 + Xx- 1 + ■ • ■ - 


where 


X(x) 


tr{x) = fl (X - X t ) = X 2 - </tX* + ( T -> , 

i-i 


(27) 
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r' i 0 t ' fl r c " t ,!» Si, f 1Wm 2. *) can bo determined by the use of a cc 

h'(-)< ■-^irixj -f - 4 m ,27) then, by turned fraction which is given below in tabular form. 


Rv-n(.v) — q s (x)Rs-i(x) 4- Rs(x) 


as(x) = q y (x)<T> «_,(x) 


1 0 

1 8.r* ■■■ 5.v ' 1I.V-' - .r‘ X.v • 4 

O.r 1 4- O.v- + Ox-' 4- Ox-' -(-••• 

- 8x-' - 5x-- 4- 11 a-* - x-4 + Xx- 5 4 ... 

2 2.v 4 . >3 

- 8x- - 5.r* [1 ' 

+ llx- 3 - x“ 4 1 1 - lO.r ' 4- 5 r 8 - o A -.a + x x ~> + ... 

Xx ~’ + • • • 10x-' ~ 5.r ! - 2x- a + Xir' + ■ • 

I O.v 1 — 15.v- — x" ' ! -I- • ■ - 


10x-- -}- 3x— 1 + Xx~ l 4- ... 


6.v — 4 


7x - 4- 3x- 3 - 8.v 1 - 5x~ 2 + llx- 3 - x~ 4 + Xx r ‘ 4- 
X.r ' -i- - • - - S.r 1 4- r » jl xx- 3 4- • ■ • 

-6x- 2 4- x.\- :t 4 - 

- 6.v-- 

- Ox - 2 - x.v -"' 


8x- 1 - 5x- 2 4- llx- 3 fx,(x‘) -0*0 4 1 - ] 


x-‘ 4- Xar* 4- 


2x -4 3 10.V-- + X- ’ 


-4 Xx~* 4- • • • 


6x — 4 0 -4 Ax :i 4- ... 


- 2x -4 3 


<r ; = (6x — 4) (2.v 4- 3) 
4- 1 

— x- — 2x 4- 9 


From the above tabular form, observe that 
K t “0 4- Xx~\ Hence, «(x) <T ,( X ) ~ x- - 2x 4- 9 where 

tr, — 2 and <r u = 9, By (7), one gets 

~ 2Fj* i -r 9E) — 0 for / > 2 (28) 

From (2(8), the rest of the transform Ej of the error 
pattern is E r , - 1, E„ = 11, E 7 = 13, £„ - £„ == 12. By 
example in Ref. 10, the inverse FFT over GF(2- 2 4 - 1) 0 f 
the Ej for / — 0,1, • • •, 7 is given by 


(eo,e u e,,e :<> e. h e rnGr „e 7 ) = (0,0,14,0,0,15,0.0) 
The corrected code word is 

b(x) = y (x) - e(x) = (5,2,9,15,2,1,2,1) 

- (0,0,14,0,0,15,0,0) 
- (5,2,12,15,2,3,2,1) 
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Appendix 

The Computation of Continued Fractions by Using Euclid’s Algorithm 


Lrt S be an irreducible rational element S in the field 
of real numbers. In this Appendix, it will be shown that 
a finite sequence of rational approximations to S can be 
constructed by using continued fractions. 

Let S = a/b, where a and b are integers, be an irre- 
ducible rational element in the field of real numbers. 
Using Euclid's algorithm, one gets, 

a = befi 4 7, (A -1 ) 

b - y iQ 2 + 

71 - 7 ^a + 73 


or 


- q\ + 


172 + 


+ . 


4 


q« 


(A-3) 


where q K for k = 1,2, * • ■, n are called the partial quo- 
tients. 

Let us define the convergents S K for k — 1,2, ••*,« as 
follows 

Si = Qi 


yk -2 d" 7 /t-i9a’ + 7 a 


7,1-:, — 7 n-zQ n-1 "f- y»-\ 
y n -2 “ 7 «-i ( 7» y>\ ~ ynQn 


or 


a 

7i 

lT~ q ' 

h 

h 

73 

— = q -2 4 


72 

7i 

7i 

73 

— = q . » d 


72 

72 


4 . J2L 

q * 

7 A-l 

7 A-l 

7»-3 

7 n-1 

— = f/n-i 

4 

7«-s 

7»-a 

7»i-2 

- — = q» 

7 «-l 



By (A-l), S = a/b can be developed by a continued frac- 
tion, as follows: 



q 3 + 


4 


C)k + «A 


, k < n 


(A-2) 


s = = q ' + T, 

1 

S., = 9i+ , 1 

q, H 

; ' <7s 

n 1 

S,i ~f) ~~ ^ ^ 


1 

-4 

<7» 

From (A-l), we observe that S A - is a finite sequence. In 
other words, S u S-, •••, S*. ••• will be terminated when 
y „ = 0. Thus, S„ = a/b, where n is a finite number. 

A recursive formula for convergents is generated as 
follows. Let = 1 and Qo = 0. Then set 


c ii^Z 1 

1 ' Q, 


q-q, 4 1 _ q»Pi 4 Bn __ 


S*j q 2 • 1 0 q-iQi + Qo 

_i (^ + i) p ' +p ° 

( ,h+ -t) Q,4 Q ° 


s 3 = q 1 + 


Q2 


+ F "(Z) 171 _ <7:. Pu + P. 

P, q ; (?2 + 

(q.g, 4- <?„) 4 


Zi 
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The recursive convergents are defined as 


Hence 


q A-i 4- Fa-- 2 _ Pa 
QkQk - l + Pa-2 ~ Qk 


(A-4) 


with 


Hat-2 — QkRk-i + Ra 


where P, = rj u P 0 = 1 , Q 0 - 0, p, = 1, for k = 2,3/ • •, n. 

(A-4) 

In order fro calculate Eq. (A-4), it is necessary to com- 
pute the partial quotients. To do this, by the same pro- 
cedure used in the derivation of (A-4), we can show first 
that S in (A-2) can be expressed in the following form: 


H-i — ( — l) -1 R!i = S 
R 0 = (-1)°R'„ = 1 

and by (A-5) 

(-l) K R h - = Q K S = P K 


(A-7) 


(A-8) 


(Qk + aK)Pii-t + Pa-2 
(Qk + a A -)Q A --i + Pa- 2 


where a h - is defined in (A-2) 


If S has this form, then 


To show that S A - ~ (Pa/Pa), which is computed by 
using a continued fraction, is an irreducible fraction, i.e., 
(Pa, Pa) = 1, consider the difference between Sa and Sa-j 
for K > 1. That is, 


[Qk 4- a a) (Pa-, ^Pa-,) ~ S Pa-2 Pa'-2 



PkQk-x ~ PaPa-i 
QkQk-i 


and, 


(A-9) 


I 

“ — = Qk T a K = 

CtA-1 


Pa - u — Qa'-'jS 

fvr- g A -.,s 


where Ra ~ P K ~ Pa-S. It follows that 


R'a--,. 

Ra-, 


(A-5) 



Finally, 


R'a - qM-i + Ra-2 (A-6) 


Let l s = P k Qk-i - QaPa-i. By (A-4), 


/a = PaQa-, ~ Pa-Pa--, - (c/a-Pa-, + Pa^Pa-, 

~ (QkQk~: + Pa-2)Pa-, 

- - (Pa-,Pa - 2 - Qk-\P A-z) 

= - 7*-i (A-10) 

Since /, — P,p 0 — PiPo — Qi * 0 — 1 • 1 = — 1, one has, 
by (A-10), /m — —7, = 1. With the above result, one has 
7 a = ( — 1) A \ It follows that 


By (5-A), the initial condition of Ra for K ~ 0,1, is 
given by 



(~1) A ~ 

QkQk-1 


for k > 1 


R' - P-, - p_,s - -S 
R' ( = P 0 — P 0 S — 1 — 0 • S = 1 

Define a new function R A in forms of Ra by Ra = 
( 1) A Ra. Then (A-6) becomes 

- (-1)^Ra-2 - q K (-iy-'R K ., - (-1)*Ra 

It is evident that 

Ra -2 = QkRk-i + Ra for both even and odd k 


or 

Pa-Pa-, - Pi. P a-, - (-1)* fork > 1 (A-ll) 

If (Pa, Pa) = d K , then, by (A-ll), d K \( — l) K . This im- 
plies that (Ik — 1. Hence (Pa,Pa) — 1. 

A simple example, showing how to compute the rational 
approximations to an irreducible rational number, is pre- 
sented in the following tabular form. For this example, 
S is the fraction 38/105. 


72 
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1 




0 

1 


0 




1 

0 


1 

38 „ 38 

105 105 

0 

38 

105 

P, = 0*1 + 0 = 0 = q. 

= 0*0+ 1 " 1 

Si — 1 

2 

38 29 

1 "105 105 

2 

29 

105 

P., = 2*0 + 1 = 1 

Q., = 2 • 1 + 0 = 2 

s. = T 

3 

38 29 9 

105 ~ 1 * 105 + 105 

1 

9 

105 

P 3 = 1 * 1 + 0 = 1 

Q 3 = l*2 + 1 = 3 

S ’=T 

4 

29 ^ 9 2 

T05 _3 ‘l05 + I05 

3 

2 

105 

P, = 3 • 1 + 1 = 4 

Q, =3-3 + 2-11 

4 

S, - n 

5 

9 2 i 

l05 ” 4 ‘T05 + 105 

4 

1 

105 

P, = 4 • 4 + 1 = 17 

= 4*11 + 3 = 47 

• 47 

6 

2 1 

- o 0 

105 “ 105 

2 

0 

P.i = 2 • 17 + 4 = 38 

Q, =2*47 + 11 = 105 

S - — 
11 105 


From the tabular form when s — n — 6, one observes 
Rti “ 0. By (A-8), 

_ F, ; _ 38 
S _ Sfi ” Q„ ~ 105 


For a more detailed discussion of the relation of Euclid’s 
algorithm to the continued fraction associated with a 
rational element in the field of real numbers, see Ref. 12, 
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The Fast Decoding of Reed-Solomon Codes Using 
High-Radix Fermat Theoretic Transforms 

K. Y. Liu and t. S. Reed 

University of Southern California 

T. K. Truong 

TDA Engineering Office 


Fourier-like transforms over GF(F„), where F„ = 2-'" *t- 1 is a Fermat prime, 
have found application in decoding Reed-Solomon codes. It is shown here that 
such transforms can be computed using high-radix fast Fourier transform (FFT) 
algorithms requiring considerably fewer multiplications than the more usual radix 
2 FFT algorithm. A special 256-symbol, 16-symhol-error-correcting, Reed- 
Solomon (RS) code for space communication-link applications can be encoded 
and decoded using this high-radix FFT algorithm over GF(F S ). 


t. Introduction 

Recently, Justesen (Ref. 1) and Reed, Truong, and 
Welch (Refs. 2, 3) proposed that transforms over GF(F„) 
(Refs. 4, 5) can he used to define Reed-Solomon (RS) 
codes (Ref. 6) and to improve the decoding efficiency of 
these codes. The transform over GF(F„) is of the form 

Z a{Qy" for 0 < / < <2 — 1 (I) 

Ml 

where F„ - 2-" + 1 is a Fermat prime for n < 4. In 
(1) the transform length d divides F„ — 1, a(t) e GF(F„), 
and y is a primitive dth root of unity which generates 
the d element cyclic subgroup 


G,= { 7 .yV>/'U} 

in the multiplicative group of GF(F„). The inverse trans- 
form of (1) is 

a(t) - ( f /)- 1 1: A(f) y-f for 0 < t < d - l ( 2 ) 

/ o 

where (d) denotes the residue of d modulo F„ and (d) 1 
is the inverse of (d) in GF(F„). 

To transform longer integer sequences over GF(F„), 
from Ref. 5, one can use the fact that y ~ 3 is a primi- 
tive element in GF(F„). Such a y gives a maximum trans- 
form length of 2-\ 
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In space communication links, it was shown in Hof. 7 
that in tho concatenated E ~ 16-error-corrocting, 255- 
symbol RS code, each symbol with 8 bits and a K - 7, 
rate L* or hi, Viterbi decoded convolutional code, can 
be used to reduce the value of E h /M„ required to meet a 
specified bit-error rate F,„ where E, t is the received 
energy lor each bit, and .V„ is the noise power spectral 
density at the receiver input. 

Figure I presents a curve of concatenated code bit 
probability of the error bound vs E t ,/N u for a K 7. 
II convolutional code with 8 bits per RS symbol. 

Since 3 is an element of order 2*" in GF(F„) (Ref. 5), 
an RS code of as many as 2" symbols off) bits or eh can 
be generated in GF[F). Hence, by Ref. 2, the Fermat 
theoretic transform over GF( F ;i ) can be used jo decode 
an RS code of 2' symbols. For a given 223 information 
symbols, each of 8 bits, as mentioned above by Ref. 2, 
221 information symbols in GF(F :f ). i.e., S t - 0, S.j, •• -.S .,,, 
can be represented in the range from 0 to 2- 1 1. After 

encoding the information symbols, the parity check sym- 
bols in the 2.58 -symbol MS code may occur in the range 
between 0 and 2 J . If 2 J ‘ is observed as a pantv check 
s\ mbol, deliberately change this value to 0. now an error. 
The transform decoder will correct this error automa- 
tically. lienee, the RS code generated in GF(F,) can be 
used to concatenate with a K 7, rate * s or ht convolu- 
tional code. 

rhe ai ithmetie used to perform these transforms over 
GF\Fj requires integer multiplications by powers of 3 
and integci additions modulo 1' However, integer mul- 
tiplications by powers of 3 modulo F„ are not as simple 
as multiplications by powers of \/2 modulo F„. which 
can be implemented by circular shifts (Ref. 5). To 
icnicdy such a problem, it is shown here tliat high-radix 
fast Fourier transform (FFT) algorithms can be list'd to 
reduct' the number of multiplications required for trans- 
forming integer sequences in GF(F„). 


II. High-Radix FFT Algorithms Over GF(F„), 
Where F„ Is a Fermat Prime 

In order to develop high-radix FFT algorithms over 
C-F(F„), it is desirable, as we shall s<>c, that multiplica- 
tions involving the 2'th roots of unity in GF(F„) be 
simple operations. This is made possible from the 
fact that tht- 2'th roots of unity over GF(F„), where 
2 < / < n •• 1 are plus or minus power of 2 mod F„. 


To see this, note that if 2|s, then 

(rh2*/ a )*s=2* mod F„ 

and 

- 2-" • 2" - - 2" mod F„ 

Hence, by theorem 2.20 of Ref. 8, the conf nonces 

**'--2* mod F„ (3) 

and 

.v- - —2'* modF„ (4) 

have exactly two solutions given by 

.v - ±2* - mod F„ (5) 

and 

.v . mod F„ (6) 

respectively. Xow let y be a primitive c/th root of unity 
in GF(F„), where d ~ 2' with 1 < t < 2". Then by 

theorem 1 of Ref. 9. 

'/ 2 -- —I. mod F„ 

Also by (6), 

7 1 1 : “ 1 mod F„ 

Combining (5) and (6), one obtain; 

'/ s ~ ±2^'* mod F, 

where k ■ 1, 3. By repeatedly applying (5) and (6) in 
this manner, one lias finally 

y‘ *’ ■’*' mod F„ (7) 

where 2 < i < n -f J and k - 1, 3, 5, • • •, 2' 1 — 1 . 

The high- radix FFT algorithms over GF(F„) are 
similar to those over the field of complex numbers 
(Reis, 10, 11). The following example illustrates the 
radix 16. decimation -in -frequency, twiddle-factor FFT 
over GF(F :[ ). 

Examp Let F, 2- 1 * 1 257, d 16 “ = 256 . The 

radix 16, decimation -in-frequenev, twiddle factor, FFT 
algorithm over GF(F S ) is described as follows. 
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Now the choice of y — 3 1 gives 


Let f and t in (1) be expressed ns 


/ = /. • 16 f /„ 

(8) 

/ — i, • 16 4 /„ 

(9) 


where 

f u t, =0,1, 2, 3, ••■,15 
Substituting (S) and (9) into (1), one has 


A(f) = £ E a ( tl ‘ 16 + l,i+/o) ( ' 1 ‘ lMo) 

tu u r, <i 

( 10 ) 

Since y J = y ,o! ^l mod F 3 , (10) becomes 


A(f)= £ jr «(/, • 16 i h,)y /,<# , ^ /utl 1<>,/uf0 
£ «(*i • 16 t- L)y'"' 1 ‘«1 y'“'“ 

.o ■■ J 


L r ' 


Jlto 1 « 


Let 


16 i U = 


!>(/,• 16 + *„)/•■• l,; y /oi ° 


Btf.-lfl + 16 f 

i„ ii 

The radix 16, 256-point, FFT algorithm over GF(F.,) 
is then composed of the following stages: 


Stage 1 ; 


B,tf..'16 + f (1 ) - 
Stage 2: 

Bj(fn • 16 + /,) = 

From (7), 


£«(*,• 16 f Uy'"' 1 1,1 


y/0^l 


EB.(/.,’16 I- L)y r ’'" " 


(ID 


( 12 ) 


y“> =n s ±2* mod F, 

where A' -- 1, 3, 5, 7. It is shown in Ref. 9 that if y is a 
primitive element in GF(q), where q is a prime, 
then w"‘ is also a primitive element in GF(q ), where 
m 3,5, •••jj — 2. It is well known (Ref. 5) that 3 is 
a primitive element in GF(F„). Thus 3 W is also a primitive 
element in GF(F„) for m -- 3, 5, ■■ ■ 2-" — 1. 


y" 1 - (3 a ) , “ ^ • (3 IM ) S = (■-2 3 ) a ^2 mod F ;i 

since 3"%,: - •2 :| (mod F.,) and 2'L :-. —2 (mod F ■). Hence, 
a: j n (li) can take on only the values r*-l or a power 

of 2, 

Since multiplications by ±:l involve only sign change 
and since multiplications involving powers of 2 mod F :i 
can he achieved by circular shifts, the 16-point discrete 
Foumci transform in the brackets of (11) can be evaluated 
without multiplications. These results are then referenced 
by multiplying by the so-called twiddle factor 
Using a similar argument, (12) can also be evaluated 
without multiplications. 

The number of operations required to perform a FFT 
of 256 points using the radix 2, the radix 4, and the radix- 
16 FFT algorithms over GF(F„) is shown in Table 1. 
From this table, one can see that the radix 4 and tin* 
radix 16 FFT algorithms require 30% and 70% fewer 
multiplications, respectively, than the more usual radix 2 
FFT algorithm. 


In the above example, it was shown that one can find 
a power of 3 for y such that 

y' 1 ' s- 2 mod F„ 


For tli is y, one has 

y' 1 *' n,! \'2 mod F„ 


From Ref. 5, 


\/o 2 ^(2-. d mod F tl 

Hence multiplications involving integral powers of y 1 - n,! 
can he accomplished either by circular shifts or a circular 
shiit followed by a subtraction, depending on whether 
an even or an odd power of \/2 is involved. As a conse- 
quence, a high radix FFT up to 2“'- also could bo devel- 
oped. For example, the 256-point FFT over GF(F) could 
be computed using a mixed radix FFT of radix 32 and 
radix 8. 


In light of the above discussion, when transforming 
long integer sequences in GF( F„), it is desirable to per- 
form as many high-radix FFT iterations as possible in 
reduce the required multiplications. 
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Table 1. Number of operations required to transform 
d - 256 points FFT over GF(F„), where n - 3, 4. 


Algorithm 

Mod F n 
multiplications 

Mod F„ 

additions 

Circular 

shifts 

Radix 2 
(d = 2 8 ) 

769 

2048 

0 

Radix 4 

Id = (*' '*■ 

513 

2048 

256 

Radix 16 

id - m-i 

225 

2048 

544 


SPACE NETWORK PROGRESS REPORT 42*36 





N77-14058 


DSN Telemetry System Performance With 
Convolutionally Coded Data Using 
Operational Maximum-Likelihood 
Convolutional Decoders 


B. Benjauthrit and B. D. L. Mulhall 

TDA Engineering Office 

B. D. Madsen 

Telecommunications Systems Section 

M. E. Alberda 

DSN Data Systems Section 


This article describes the DSN telemetry system performance with convolution- 
ally coded data using the operational maximum-likelihood convolutional decodei 
(MCD) being implemented in the Network. The report covers data rates from 80 
bps to 115.2 kbps and both S- and X-band receivers. The results of both one - and 
two-way radio losses are included. 


I. Introduction 

DSN telemetry system performance with convolution- 
ally coded data at low data rates was presented in Ref. 1. 
In that article, real-time telemetry data obtained from the 
Compatibility Test Area (CTA 21) were processed by a 
simulated Viterbi decoding computer program to measure 
the system performance. The results covered the data 
rates ranging from 8 to 2048 bps, which will be useful in 
the design of telecommunication links for future space- 
craft such as the Mariner Jupiter-Saturn (MJS), to be 
launched in 1977. Since such spacecraft will also employ 
higher data rates, the work described in this article cover- 
ing data rates from 80 bps to 115.2 kbps was undertaken. 
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The article has a threefold purpose. The first is to give 
overall MCD acceptance test data received from the 
manufacturer. The second and third are to describe the 
test purposes, test setups, and test results from CTA 21 
and from the Telecommunication Development Labora- 
tory (TDL). The tests concentrate mainly on the Block IV 
receiver, high-data-rate, X-band carrier frequency with 
some S-band results to allow S- and X-band performance 
comparison. The effects of using different receiver band- 
widths and subcarrier demodulator assembly (SDA) band- 
widths are also investigated. 

All 32 MCD production units were acceptance-tested, 
and their required performance was verified at the highest 
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data rate, 250 kbps, before delivery. They were also 
checked with the telemetry processing assembly (TPA) 
for proper interface prior to installation. 

Since the actual operational software for the TPA was 
not ready when this project was underway, a special real- 
time computer program, called MCD Performance Evalu- 
ation Program (MCDPEP), was written in assembly 
language for collecting decoded data and formatting 
tli cm into an appropriate compressed form on disc. The 
MCD status and symbol error information from the sym- 
bol synchronizer assembly (SSA) was also recorded. The 
recorded data were then processed by an off-line MCD 
Data Analysis Program (MDAP), also written in assem- 
bly language, to provide bit error rate (BEK), symbol 
error rate (SER), burst error statistics, and MCD normal- 
ization rate (NR). The burst error statistics are expressed 
in terms of error-free runs ( EFR ) and ourst length sample 
distributions. 

The two-way radio losses of the DSN telemetry system 
were also examined for various values of carrier loop 
signal-to-noise ratio margin. These latter test results were 
obtained at TDL. 

At present, only certain tests at data rates 44.8, 67,2 and 
115.2 kbps are available from CTA 21. However, the test 
results from TDL cover a larger data range from 80 bps 
to 115.2 kbps. Additional tests will be accomplished in 
CTA 21 and in Spain, at DSS 62, by DSN Madrid Engi- 
neering. This will allow station-to-station comparison. 

II. MCD Acceptance Test Data 

Before delivery, each production MCD was verified to 
satisfy the functional requirements. The telemetry test 
data used were assumed corrupted by additive white 
aussian noise, and the tests performed at the highest 
allowable data rate, 250 kbps. 

The acceptance tests consisted in part of BER per- 
formance, MCD symbol synchronization recovery inter- 
val, and channel error rate performance. These were 
accomplished for both code rates % and %. Table 1 
summarizes an average performance in terms of bit signal- 
to-noise ratio (ST d /N„), BER, and channel error rate. The 
table includes the corresponding standard deviation fig- 
ures which are noticeably small. The functional require- 
ments, lower and upper bounds, are also given in Table 1. 

The verification of MCD symbol ^synchronization re- 
covery interval was done at SV B /N 0 = 3 dB. For each 
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MCD, ten observations were made for each code rate. 
This gives rise to 320 observations in all. The histogram 
plots of both code rates are depicted in Fig. 1. The 
figure indicates that the MCD time in symbols to re- 
synchronization remains well below the functional re- 
quirements, i.e., 2000 symbols for rate Vz and 3000 
symbols for rate Vb. 1 The corresponding cumulative dis- 
tribution plots of Fig. 1 are given in Fig. 2. 

III. Performance Tests at CTA 21 

The purpose of the test investigation at CTA 21 is to 
obtain telemetry system performance for medium to high 
data rates under the Mark IIM977 DSN telemetry system 
configuration. Emphasis is placed on the X-band Block 
IV receiver with corresponding results for the S-band, 
Block III receiver. Effects of receiver and SDA band- 
width are also examined. Another objective is to determine 
an optimum modulation index for each data rate. 

The equipment used is a typical telemetry string (re- 
ceiver, SDA, SSA, MCD, and TPA), along with symbols 
provided by the SCA via a radio frequency link. The SCA 
piovides a cyclical repetitive 2047-bit pseudo random 
noise (PN) sequence encoded by a 7:Yz (or 3 / 3 ) convolu- 
tional encoder (with alternate symbols inverted) at the 
desired symbol rate and modulation index value. The 
MCD receives the noisy quantized (5-bit parallel) sym- 
bols from the SSA, performs its decoding function, and 
outputs decoded data to the TPA. The SSA, after proper 
initialization from the TPA via commands transmitted by 
the 920 emulator, periodically (once every 10° symbols) 
outputs symbol error count (SEC) to the TPA via the 
SSA/TPA coupler. Thus, MCD decoding performance is 
obtained from output SSA SER and output MCD BER. 

The MCDPEP (resident in the TPA) stores and eval- 
uates all MCD output data in real-time. Data are received 
m 8-bit bytes and sequentially stored in two alternate 
memory blocks (200 16-bit words/block), using the ex- 
ternal direct memory processor (DMP) in the data chain- 
ing mode. Data are processed from one block while the 
other block is being accumulated, Data processing starts 
by searching the stored data to find frame synchroniza- 
tion within the known 2i47-bit PN sequence, consisting 
of 11 consecutive Is. Normal data, and inverted data if 
necessary, arc searched. After finding frame sync, every 
stored data bit is compared with the corresponding known 


'This assumes a nominal bit transition density ( 1/40 - 1/2 at 3 dB 
and greater). For low bit transition densities, the result muv not 
be valid; see Ref, 2. 
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PN bit, being clone 8 bits (one byte) at a time. Every byte 
containing one or more erroneous bits is stored in an 
output error list (OEL), and periodically transferred from 
core to disc in 50-word blocks, 32 bits/word. 

The exclusive-OR byte pattern is stored, with I s indi- 
cating error bits and 0’s indicating correct bits. The corre- 
sponding good byte count (GBC) is also stored in the 
same entry, providing the number of consecutive correct 
bytes since the last previous error byte. 

The 32-bit word format used for each OEL entry is 
shown below. 


GBC (16 bits) XOR (8 bits) 

NR (7 bits) 

1 bit 

Good Byte 
Count 

Byte Error 
Pattern 

MCD 

Average NR 

MCD 

Sync Change 


An entry is made in the OEL whenever one of the 
following events occurs: 

(1) GBC Full Count. An entry is made any time this 
counter reaches full-scale value, which is 8 X 2 ltl or 
524,288 consecutive correct bits. 

(2) MCD BEK, Once every DMP block (or 3200 bits) 
the MCD is commanded to output its 4-bit NR word, 
This number is derived from the MCD internal 
normalization rate. It is related to the input SER 
and consequently the output BER. These NR 
outputs are stored and averaged (over 31 outputs), 
and the average value is entered in the OEL in the 
7 bits provided (4-bit whole number plus 3-bit 
fractional number) once every 99,200 data bits. For 
these NR entries, the GBC value is shown as 0, 
although the actual value is retained and continues 
to be accumulated. 

(3) MCD Node Sync Change. When a node sync 
change occurs, the MCD outputs its own status and 
an entry is made in the OEL, with the least signifi- 
cant bit equal to 1. If the MCD had previously 
achieved symbol sync, and then a sync change 
occurs; this means that the MCD output data are 
temporarily not valid and the program reverts to 
its frame sync search mode. When frame sync is 
again found, this implies that the MCD has re- 
gained symbol sync and the output data are again 
valid. 
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The MCD generates a node sync change when its input 
symbols are uncorrelated. To resync, the MCD attempts 
to find an acceptable level of symbol correlation by re- 
grouping the incoming symbols. 

The OEL continues to accumulate and is stored on disc 
for the test duration, nominally 10 7 bits. At the same time, 
another list is output on the Terminet printer containing 
every other average NR value, or once every 198,400 bits. 
When this list has reached the desired length, the test is 
manually terminated by actuating a switch on the TP A 
console, after which the stored SEC list is printed on the 
Terminet, with one entry every 10“ symbols (or any 
other power-of-ten number of symbols available ) . 

MDAP is an off-line program for use to complete the 
MCD data analysis phase. The program converts the OEL 
entries produced and recorded on disc by MCDPEP into 
a sequence of EFR’s, classifies them, and generates their 
sample distributions. It also generates average BER and 
average NR at every 99,200 decoded bits. At the end of 
each run, the following are printed: 

Total number of decoded bits 
Total number of bit errors 
Average bit error rate 
Total number of EFRs 
Total number of bit bursts 
Total number of normalization rates 
Average bit errors per burst 
Average burst length in bits 
EFR sample distribution 
Number of distinct EFRs 
Burst length sample distribution 
Number of distinct burst lengths 
Symbol error rate 

When a sync change bit is detected, MDAP outputs the 
entire data statistics, reinitializes program parameters, 
and restarts to process the remaining data again. 

Figure 3 gives an overall view of the test configuration. 
(A more detailed configuration is given in Fig. 7 which 
will be discussed later). The ST n /N 0 used in describing 
the system performance is obtained from the average SER 
accumulated for each run via an uncoded performance 
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curve. The MCD error rate calibration curve (MCD input 
STb/No in dB vs MCD normalization rate indication) for 
each rate is also obtained from the program printout, as 
shown in Fig. 4, 


A. Optimization of Modulation Index 


The signal energy per bit (which determines the bit 
error probability of a telemetry channel) of any modula- 
tion system is related to the total received power P? as 
follows: 


ST, _ P T P D 

~-K:j7" nln '' ,,s 

where 


( 1 ) 


Pn = sideband power available for data 

= single-sided noise spectral density, W/Hz 
T ,) = bit time in seconds 

ywDL ~ subcarrier waveform distortion loss 
■q, = system efficiency 

= ytu. + i}sm. 4- 'j/ix/,/.. dB, where 


If the vaiues of P r /N, h and BPS are known, the plot 
of STn/N lt vs 0 can be obtained from Eq. (2). Such plots 
for various values of P r / N„ from 51.45 to 54.58 dB are 
given in Fig. 5. They are for 10- and 30-1 Iz Block IV 
receiver and wide bandwidth Block III SDA. From these 
curves, one can determine the optimum modulation index 
for each P r /N„. To derive this value algebraically, one 
first expresses Eq. (2) as 

II 

A T n Vs 

BPS — —fTr , — sin-’tf 


The maximum data rate, BPS maj , can now be computed 
by setting the derivative of BPS with respect to 0 equal to 
zero, yielding 


0o,„ = tan- 1 dvs > 

f j ^ = Qopt 

This indicates that the optimum modulation index for 
a single-channel system is totally dependent on the proper 
characterization of the system losses. 


yitL = radio loss, dB 

t/n/j/, = subcarrier demodulation loss, dB 

yii.sai. = bit sync and detection loss, dB 

One way to optimize the performance of the telemetry 
system is to maximize the data rate. To achieve this, the 
values of P r /N „ and S7',/N u must be specified. Then the 
parameters P n /Pr and v* arc adjusted to obtain the maxi- 
mum data rate. Another way is to specify the values of 
the data rate, STn/N„, and adjust the parameters Pn/Pr 
and v s to minimize the required Pr/N,,. 

The easiest parameter to adjust in the telemetry system 
is the modulation index. This parameter determines the 
allocation of power in the various sidebands and strongly 
affects the magnitude of the system losses. 

For a single channel, a square-wave subcarrier, Pn/Pr, 
is sin-'#, where 6 denotes the system modulation index. 
Since fjwnr. is usually small (?ni /a, — 0 dB), Eq, (1) may be 
rewritten as 

ST, P T Vg 

v si* 3 *, BPS = 1/T„ (2) 


B. Test Strategies and Preliminary Results 

The end product of the project is to obtain DSN telem- 
etry system performance using convolutional! y coded data 
in terms of hit error probability (BER) vs ST lt /N ( , and the 
corresponding burst characteristics. Toward this end, we 
first proceed to determine the optimum modulation index. 
This is to be used to corroborate the theoretical value of 
the true optimum modulation index. Once this parameter 
is determined, the telemetry system performance for it is 
then obtained. 

The results obtained so far have been for data rates 
44.8, 67.2, 96.8, and 115.2 kbps. Figure 6 describes the 
telemetry system performance at data rate 67.2 kbps. The 
test conditions are given in Table 2. The true optimum 
modulation index for this rate appears in Fig. 6(a) to be 
about 80 deg, which is in agreement with the theoretical 
value. The curves BER vs ST ,/N,, corresponding to modu- 
lation indices 79 and SO deg are given in Fig. 6(b). The 
sample distribution plots of the hurst lengths are given in 
Fig. 6(e)~(f). The average burst length vs modulation 
index and the average hurst error vs modulation index are 
given in Fig. 6(g). The average burst length and average 
hurst error vs ST„/N n arc given in Fig. 6(h). And finally, 
the MCD error rate calibration for this data is depicted 
in Fig. 4. 
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IV. Performance Tests at TDL 

One purpose of the tests conducted ut TDL is to mea- 
sure system two-way radio losses using an operational 
MCD over various data rates for both S- and X-band 
receiver. Another purpose is to obtain telemetry system 
performance in terms of BER vs S7 »/No. 

Equipment used for the tests includes the Block III 
exciter. Mariner Venus/Mercury transponder, X-band 
transmitter, Block III receiver (S-band), Block IV receiver 
(S- and X-band), Block LU SDA, Block ill SSA, and a 
production model MCD. Test data patterns include 31-bit, 
127-bit, and 2047-bit PN codes. The detailed test equip- 
ment setup is shown in Fig. 7. 

Figure 8 describes system two-way radio loss curves 
derived from measurements taken in TDL during January 
1976, for code rate 7 : Vs, 3 bit quantization. The figure pro- 
vides BER vs ST a/ No for various values of the uplink mar- 
gin, M c = P r / N 0 2 B lAh where P c is the carrier power and 
B l o is the single-sided design point phase-lock-loop band- 
width. Tliis parameter is selectable for each test from the 
dial at the receiver panel. Here the theoretically optimum 
downlink modulation indices, discussed in Section 111, 
were used for all tests. This means that the ground receiver 
loop SNR was not held constant but was allowed to vary 
with signal level as it will be in an operational situation. 
These results are not exactly the same as will be obtained 
with the MJS’77 transponder, due to differences in receiver 
loop parameters, but are close enough to provide a basis 
for prediction of the MJS’77 losses. 

The most obvious conclusion to be drawn from the data 
is that the X-band two-way radio losses are so great that 
they will impose operational limitations on the MJS 77 
mission. In order to limit the loss to about 1 dB, the 
uplink carrier margin must be 35 dB or greater. This is 
equivalent to an uplink carrier power of *-120 dBm, 
which is just about what can be achieved at Saturn, range ~ 
10 AU, with a 100-kW, 64-m DSS. Operation at greater 
distances or with a 26-m DSS will require noncoherent 
operation of the X-band downlink. 

The operational thresholds arc also shown in Fig. 8. 
Further, the X-hand telemetry threshold as a function of 
uplink margin is described in Fig. 9. 

The difference between the X-hand and S-band per- 
formance is believed due entirely to the greater turn- 


around phase jitter on the X-band signal due to the VCO 
phase noise being multiplied by a greater ratio (11/3 
times). The magnitude of the problem can be reduced by 
increasing the frequency of the uplink signal (thus re- 
ducing the multiplication ratio) or by using a feed- 
forward technique for removing most of the phase jitter 
from the downlink. Using a feed-forward system has been 
demonstrated in the laboratory and resulted in almost no 
two-way degradation at uplink margins as low as 5 dB. 

It is strongly recommended that both of the above solu- 
tions be investigated. It is probably too late to do anything 
about MJS, hut the standard transponder should cer- 
tainly include the feed-forward features. 

The TDL MCD error rate calibration curve was pre- 
viously given in Fig. 4. The telemetry performance curves 
in terms of BER vs STn/No for rates 7.2 and 115.2 kbps 
are given in Fig. 10. 1 he performance of code rates 1/2 
and 1/3 at data rates 115.2 and 76.8 kbps, respectively, is 
given in Fig. 11. For BER ’. s of 10 to 10 4 , the rate 1/3 code 
is 0.3 dB better. For a BER of 10 r \ the improven is 
1.25 dB. 


V. Telemetry Performance for Radio Frequency 
Subsystem 

Telemetry performance for the Radio Frequency Sub- 
system (RFS) using both the TDL telemetry simulator and 
the Mark III Data System (MDS) as sources is listed in 
Table 3 and plotted in Fig. 12. Telemetry performance is 
similar to that obtained previously during developmental 
tests and no anomalies were found. 

The high rate (115.2 kbps) two-way telemetry loss at 
X-band was found to be slightly worse than that mea- 
sured on the MVM’73 RFS. This was expected since the 
MJS 77 RFS has a wider bandwidth. The measured losses 
for both transponders are shown in Fig. 13. The effective 
telemetry loss at a threshold BER of (5.0 X 10 ') as a 
function of uplink carrier margin is shown in Fig. 14. This 
shows that an uplink carrier margin must he at least 35 dB 
whenever high-rate X-band telemetry is being transmitted 
in the two-way coherent mode. With the present DSN 
capability, this characteristic precludes use of the two-way 
coherent mode with high-rate telemetry for ranges greater 
than 10,0 AU for a 100-kW, 64-m station or 2.0 AU for a 
20-kW, 26-m station. 
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Table 1. JPL MCD «cceptanc- test data 
(averaged over 32 units) 


Code rati’ SIji/N,, 


MCD-Tt’A 1/2 3-0 

■ 1.0 
5.0 

1/3 3.0 

3.5 
■1.5 

MCD-DDA 1/3 3.0 


Bit ertor rate 


No, bits 

■1 ,096 X l() n 
•1,096 X 10 7 
6.553 X 10 7 
•1.096 X 10" 
•1,096 X 10 T 
4.096 x 10" 
4.096 X 10" 


Mean 

7.41 X 10 ' 
3.95 X l()- 6 
1.36 X 10'* 
2.90 X 10- 1 
0.92 X 10-* 
2,83 X 10 n 
2.88 X 10 1 


Std. dev. 

1.82 X 10 1 
0.24 X 10 * 
7.95 X 10 » 
2,47 X 10 7 ’ 
2.77 X 10 " 
5.11 X 10 7 
2.41 X 10- s 


Upper 

bound 


Channel Enor rate Lower 
mean std. dev. bound 


Upper 

bound 


9.0 X 10 1 

11.5 

0.51 

6.0 X 10 71 

7.9 

0.62 

2,5 X 10 

5.0 

0.33 

3.6 X 10 ' 

11. 1 

0.45 

1 .0 X 10 ' 

9.7 

0,44 

7.0 X 10 11 

7.1 

0.24 


10 13 

7 9 

4 3 

lO 12 

9 H 

6 8 


3.6 X 10 1 


Test 

ID 


DAI 

2 

3 

4 

5 

6 

7 

8 

DB1 

2 

3 

4 

5 

6 


8 

DC 1 
o 

.3 

4 

5 

6 
i 

8 

DD1 

o 

3 

4 

5 
0 

7 

8 


Table 2 


, R 7 2-kbns Bloch IV receiver, S-band 

Test conditions and results for 67.z sop 


Total bits Mod index. 
(X LO 7 ) lle & 


UVC IV 
BW, Hz 


1.51 

77 

1,42 

78 

1.21 

80 

1.38 

79 

1.81 

79 

1.41 

79 

1.08 

79 

Cl 

1.10 

( w 

1,35 

77 

1,30 

78 

1,35 

80 

1.33 

79 

1.58 

79 

1.06 

79 

1.00 

79 

1,57 

79 

1.24 

77 

1.30 

78 

1,37 

79 

1.47 

80 

1.58 

80 

1.04 

80 

1 .05 

80 

1.30 

80 

1,30 

77 

1.27 

78 

1,33 

79 

1,31 

80 

1.54 

80 

1.12 

80 

1.15 

80 

1.10 

80 


SDA III 
BW, Hz 


Medium 


Narrow 


Medium 


Narrow 


AVE 

BEK 


1.24 X 10 1 
1.48 X 10 1 
1,12 X 10-' 

1.25 X 10 ' 

1.55 X 10 < 
8.00 X 10 r> 
2,70 X 10 7 

0 

1.31 X 10 * 

1.10 x io 1 

6.56 X 10 » 
1.39 X 10 1 
6.55 X 10 • 
l .00 X 10 ' 

8.10 X 10 7 
0 

1.35 X 10 ' 
1.59 X 10 ' 
1.53 x 10 ' 

1.82 X 1(H 
4.98 X 10 r ' 
6.48 X 10- 7 ’ 
2.70 X 10 7 

0 

7.57 X 10 
1.13 X 10 ' 
8.72 X 10 r ’ 
1.32 X 10- * 
6.44 x 10 5 

9.82 X 10-' 
0 

0 


AVE 
SER, % 


5.97 
6.07 
5.60 
6.05 

5.98 
5.90 
3.02 

1.46 

5.99 
0.03 

5.46 
0.25 
5.76 
5.91 
2.97 
1.44 

0.42 

0.41 

0.41 

6.42 

5.75 

5.84 

2.95 

1,32 


6.20 

5.89 

5.74 

6.03 

5.74 

5.92 

2.97 

1,35 


MCD input 
ST,,/ N„. dB 


3.85 

3.80 
4.03 

4.81 

3.85 

3.85 
5.48 
6.77 


3.84 

3.82 

4.09 

3.72 

3.95 

3.88 

5.51 

0.80 


3.64 

3.65 
3.65 
3.04 
3.95 
3.91 
5.52 
0.93 


3.74 

3.89 
3.96 
3,82 
3.96 
3.87 
5.51 

6.90 


AVE 

MCD NR 


7.75 

8.00 

7.75 

8.13 
7.88 
7.50 

3.13 

1.13 

7.50 

7.75 

7.25 

8.25 
7.38 
7.75 

3.13 

1.00 

8.38 

8.38 

8.63 

9.00 

7.13 

8.00 

3.13 

1 .00 

8.00 

7.00 
7.75 
8.25 
7.75 
7.88 
3.13 

1.00 
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Table 3. Measured telemetry performance 


Band/test 

M.I., deg 

Bit rate, 
kbps 

STn/No, 

dB 

HER 

Loop 
SNR, dB 

XI 

80.0 

115.2 

3.84 

4.7 x in i 

15.4 

X2 

80.0 

1 15.2 

2.84 

5.0 X KH 

14.4 

X3 

80.0 

1 15.2 

5.20 

2.1 X 10 3 

15.2 

X4 

80.0 

44.8 

4.84 

1.8 X 10* 

13.1 

X5 n 

79.5 

1 15.2 

3.08 

2.0 X 10 - 3 

15.2 

X6« 

79.5 

44.8 

3,13 

U X 10- 3 

12.1 

X7“ 

79.5 

29.9 

3.09 

1.5 X 10- 3 

ll.fi 

X8" 

79.5 

21.fi 

3.00 

2.0 X 10- ; ‘ 

10,6 

X9» 

72.0 

7.2 

3.06 

1.4 X 10- 3 

1.9 

X10" 

72.5 

2.5C 

3.42 

2.2 X 10 -a 

8.5 

XU“ 

05.0 

1.20 

2,78 

3.5 X lO- 3 

8.5 

SI 

45.8 

2.56 

3.07 

1.1 X lO- 3 

16.9 

"These tests performed with MDS/RFS 
formed using TLM simulator for source. 

combined; 

others per- 
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Fig. 1. Sample distribution plots of MCD symbol 
synchronization recovery for both code rates 
1/2 and 1/3, ST*/N 0 - 3 dB, 250 kbps 
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Fig. 3. MJS telemetry system test configuration at CTA 21 
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Fig. 6 (contd) 
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Fig, 7. MJS telemetry system configuration at TDL 


94 


JPL DEEP SPACE NETWORK PROGRESS REPORT 42-36 



















JPLMCD \ 
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MARGIN = 30 dB 
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PLANNED DATA MODE 
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ENERGY PER BIT/NOISE SPECTRAL DENSITY, ST^/Ng (dB), BY Y-FACTOR 

Fig. 8. X- and S-band two-way radio losses with various uplink carrier margins for data 
rates from 80 bps to 115.2 kbps, Block IV, 30 Hz 
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Fig. 8 (contd) 
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Fig. 11. Bit error probability vs ST n /N u for code rate 1/2 and 1/3 


Fig. 12. Measured telemetry performance 
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Fig. 14. X-band two-way radio loss vs uplink carrier margin 
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Modification of Hipotronics Discontinuity 
Enamel Wire Test for Wire Screening 

W. D. Schreiner 

Quality Assurance DSN and Mechanical Hardware Section 


The Quality Assurance DSN and Mechanical Hardware Section and the Fabri- 
cation Section have redesigned the wire spindle section of the wire screening test 
equipment to prevent breakage of small gage magnet wire during testing 
operations. 


t. Introduction 

In the last year, use of magnet wire at JPL has increased 
dramatically. With emphasis being directed toward 
smaller and lighter packaging of electronics, the use of 
lighter gage wire has risen sharply. 

Prior to its use, the wire is screened by Quality 
Assurance for defects. In attempting to perform screening 
tests with our existing test equipment on the smaller gage 
magnet wire, problems were encountered with repeated 
breakage of wire. After investigation and analysis, a 
deficiency was discovered in the design of the spindle that 
induced improper load and wire tension during testing 
operations. 

II. Description of Equipment 

The equipment used by Quality Assurance to perform 
the screening test is a Hipotronics Discontinuity Fnamel 
Wire Tester Model DT-1 (Fig. 1). This unit is designed to 
spark-test the insulation coating ot magnet wire and to 
record the number of defects. 
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The procedure is to set a spool of magnet wire on the 
left side spindle, remove a small amount of insulation from 
die end of the wire, feed the wire through the guide 
rollers and through the mercury bath, and then attach the 
end of the wire to the ground terminal of the takeup 
drum. The voltage range is selected and the machine is 
started, pulling the wire through the mercury bath. Any 
defect in the enamel coating completes the circuit 
through the mercury to ground and activates the defect 
counter. A footage counter is preset to run 30.5 meters for 
each test. 


III. Problem 

The problem experienced was frequent breakage of 
wire in the smaller gages (38 and under). This was caused 
by the manufacturer’s design of the feed spindle (Fig. 2). 
The spindle was designed with radiai bearings only, with 
no means for supporting the weight oir the spool. (A spool 
of wire weighs from 0.45 to 4.54 kg.) To adjust tension, a 
screw on the bottom required tightening or loosening. 
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I. 


This action exerted a compression load on bearings not 
designed for this purpose. Erratic operation was the result, 
causing constant breakage of wire. 


IV. Solution 

A redesign of the spindle was accomplished by the 
Quality Assurance DSN and Mechanical Hardware 


Section, and the Fabrication Section. The new design (Fig. 
3) incorporates a compression bearing below the spindle 
to carry the load. In the main body there are two radial 
bearings. This combination provides smooth, free-running 
action regardless of weight. Tension is accomplished by 
use of a spring-loaded nylon button that presses against 
the main shaft; it is adjusted by a set screw. The result has 
been an extremely smooth operation with no breakage or 
loss of wire during test. 
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Fig. 1. Hipotronics Discontinuity Enamel Wire Tester Model DT-1 


Fig. 3. Feed spindle, new design 
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Fig. 2. Feed spindle, manufacturer’s design 
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Two MBASIC Programs That Write 
Application Programs for Accessing a 

Database 

R. M. Smith 

DSN Facility Operations Section 

A method was desired to relieve the tedium of writing and testing application 
programs, Two utility programs were developed to produce application 
programs that perform relational operations on data. No coding is performed by 
the user. 


t. Introduction different approaches to the use of specialized application 

programs were described in Refs. 2 and 3. This article 

The nature of the MBASIC Processor facilitates its use describes tlle result ° f some preliminary efforts to design a 

by people who are not primarily programmers and who sim P le method of P roducin g an MBASIC application 

may never have had any previous programming experi- P™? 3 " 1 while isolating the user from the task of writing 

ence. Many managers, and other people needing to make code ‘ 

use of management data, fall into this category of user. 

However, for many data users, programming is so much Application programs that extract data from a database 
less important than data usage that it would be beneficial may be data dependent (knowledge of data organization is 

to reduce or eliminate programming in their data built into the application program, making the program 

activities. One approach to minimizing programming time sensitive to changes in data organization) or may possess 

is to make use of a generalized application program varying degrees of data independence (immunity of an 

similar to the one described in a previous article (Ref. 1). application program to changes in data organization). The 

A generalized program, once written, allows the user to greater the degree of data independence, the less the 

concentrate upon data usage rather than program writing. effect of changes in the database. The programs described 

However, a generalized program requires that the user in this article are preprocessors that accept relational 

spell out an access strategy each time the program is used. statements from a user, convert the statements to 

If a specific application is required for repeated use, then MBASIC code, and store the statements as an application 

a specialized application program is most desirable. Two program for later use. 
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II. Description of the Program ‘WRITER* 

‘WRITER' produces specialized application programs 
with a moderate amount of data independence. User input 
to the program is fully prompted and makes use of 
relational operations (Refs. 4. 5) to specify database access. 
The relational operations may be invoked in any order 
repeatedly or used singly. Figure 1 illustrates a session 
invoking a restriction and a projection in sequence. Figure 

2 illustrates the same process in relational notation. Figure 

3 is a simplified flow diagram of the program showing the 
iterative nature of the main program in accepting user 
input and choosing a specified subprogram to assemble 
program statements. 

In a typical sequence of events, the user 

1. Enters the name of the new program to be written. 

2. Enters the type of relational operation. 

3. Enters the parameters that describe the selected 
relational operation. 

‘WRITER’ then assembles a program by the following 
process: 

1. Copies a set of generalized code lines to a 
temporary program file. 

2. Appends code lines (to the temporary program file) 
that are created by ‘WRITER’ and are specific to 
the user's application. 

3. Appends standardized subroutines (stored for thk 
purpose) for each relational operation involved. 

4. Copies the temporary program to a file named by 
the user. 


Figure 4 is a copy of the program produced by 
‘WRITER’ using the process depicted in Fig. 1. Line 100 
and lines 902 through 918 are the generalized code lines 
mentioned previously. Lines 2000 through 2320 and 3000 
through 3310 are standardized subroutines for a projection 
and restriction, respectively. (Lines 5000 through 5090 


update a temporary directory relation that describes data 
files accessed by the application program.) 

III. Description of the Program ‘WRITPR* 

WRITPR' produces specialized application programs 
with no data independence. User input is fully prompted 
and makes use of relational operations (Refs. 4, 5) to 
specify database access. Figure 5 illustrates a session 
invoking a restriction and projection (see relational 
notation in Fig. 2). Figure 6 is a simplified flow diagram of 
the program. The user prompting sequence for ‘WRITPR’ 
is similar to that of ‘WRITER’ but ‘WRITPR’ differs in its 
approach to writing the application program (compare 
Figs. 3 and 6). ‘WRITPR’ assembles a program by writing 
code lines (on the program named by the user) that are 
specific to the user’s application. Figure 7 is a copy of the 
program produced by ‘WRITPR’ using the process 
depicted in Fig. 5. 

IV. Miscellaneous Information 

To produce code that is specific to the user’s applica- 
tion, both programs use “WRITE ON” statements that 
incorporate variables and “counters” into a completed 
statement for the application program. Examples of this 
process are shown in Figs. 8, 9, and 10 and are taken from 
the program ‘WRITPR’. Figure 8 shows the lines of code 
that produce lines 110 and 120 of the program presented 
in Fig. 7. The code lines in Fig. 9 produce line 130 of Fig 
7, and the code lines in Fig. 10 produce line 150 of Fig. 7. 

Each of the sessions (illustrated in Figs. 1 and 5) 
requires approximately 3 to 4 min of terminal time and 
produces programs that are fully functional, requiring no 
testing of the MBASIC code. Both provide the user with a 
uniform, extremely simple process for data access Figure 
11 illustrates the data output produced by the application 
programs written by ‘WRITER’ and ‘WRITPR’. 

Data files accessed by these programs must be in, at 
least, first normal form. The file used in this article (Fig, 
12) is in third normal form (Ref. 4). 
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RUN 

ENTER NAME QF NEW PROGRAM*. I HEW 


ENTER RELATIONAL OPERA T ION <PrP* J) OR "STOP'* : R 
SOURCE RELATION: RMS^EGPT .ASSIGNMENT 
DOMAIN NAME: OWNER 
LOGICAL OPERATOR: = 

DOMAIN VALUE • lc 
PROJECT TO: ♦PILE 

ENTER RELATIONAL OPERATION <‘P*F*J> OR “STOP": P 
SOURCE RELATION: ♦FILE 
QUANTITY OF DOMAINS: 3 

TARGET RELATION: TERM INAL 

DOMAIN NAMES: CON ,OPSTAT .LOCATION 

ENTER RELATIONAL OPERATION <F iF<J) OR ’‘STOP”: STOP 
TEST I NEW 


'RMS *EQPT. ASSIGNMENT’ I OWNER - 12 
TT TERMINAL (CON, OPSTAT, LOCATION) 


Fig. 2. Relational notation describing the 
process illustrated in Figs. 1 and 5 


Fig. 1. Illustration of a session at a terminal using WRITER’ to create an 
application program (User responses are to the right of each colon.) 



Fig. 3. Simplified flow diagram of ‘WRITER’ 
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COPY "1NEW" TO TERMINAL 

100 JOIN*1000>PROJ=2000>REST*3000»NN«44>TARG»5000 

110 COPY "RMS^ECPT.RREL" TO "♦INEVI" 

13 0 FL?S*^R«S*EePT.rtSSIGWeriT' *HXS=" OWNER" »LOS*"*' »DVS* 12-' 

140 PES='>FILE" 

150 GOSUB PEST 
160 GOSUE TARG 

170 RENAME "♦HOLB1 AS "♦FILE 
180 FL1S*' ♦FILE " »WD* 3 

190 STRING BPfCCDJ'' 'CON " *"OPSTAT" »" LOCATION " 

£10 NN*3c 

££0 GQSUB PRO.J 

902 STRING RCB<1> 

904 OPEN " ♦HOLE 1 " » INPUT » 1 

906 NT ENEF ILE< 1 > GO TO 914 

908 INPUT FRCN 1 USING "f»>"!RCIKl> 

910 PRINT RCDU> 

912 GO TO 908 
914 CLOSE 1 
916 REMOVE DIRS 

918 STOP CHAR < 13 >,♦ "END OF RUN'* 

£000 OPEN "♦STOl " »OUTPUT »3 
2010 FLGG=0 

£020 DIM Fta> »D<1J» »TC 1> »D*<1 > 

£030 OPEN D IRS ♦ INPUT »4 

£040 AT ENDF ILE<4> GO TO 2140 

£050 INPUT FROM 4 :FS<1 > *D< 1 > »T < 1 > *DS<1 > 

2060 STRING DOMadJO ^ . , __ ct _ c ,- n Tn *, flC( 

2070 IF FS<l>=FLl$ THEN WRITE ON 3:D<1>S'** :DS<.1> ELSE GO TO clOO 

£080 FLG&*FLGe+l 
£090 DOM<FLGQ>*DS< 1 > 

£100 IF FS< 1 >-FL 1$ THEN 60 TO £050 

£110 n=d<:ix 

2120 IF T < 1 > =D < l > THEN GO TO £140 
£130 GO TO £050 
£140 CLOSE 3 »4 
2150 Z=T < 1 > 

2160 HIM DN<OD> »D<1> *DS<1 > »KSCZ> 

£ 1 7 0 OPEN •' ♦STO 1 " » 1 NPUT » 3 
£180 AT ENDF ILE(3 J GO TO £240 
2190 INPUT FROM 3 slK t > »DS< 1 > 

£200 FOR 1=1 UNTIL I=QD?1 

2210 IF DS< 1 >=DPS< I > THEN DN<I>=D<1> 

££2 0 NEXT I 

££30 GO TO £190 

££40 CLOSE 3 

££50 OPEN FL1S » INPUT *3 

£260 OPEN "♦HOLD1" * OUTPUT »4 

££70 AT ENDFILE<3> 60 TO £310 

££80 INPUT FROM 3:K$<J> FOR J«t TO Z 

££82 FOR 1=1 UNTIL I*OD+l 

££84 FOR J*1 UNTIL J=Z*1 

£285 IF I COD THEN CH-NN ELSE CH-13 

£286 IF J=DN< I > THEN WRITE ON 4 USING " < " :KS(. J> iUHHR <CH> 

2288 NEXT J 
2290 NEXT I 
£300 GO TO £280 
£310 CLOSE 3 »4 
2320 RETURNS 


Fig. 4. Data independent application program produced by using ‘WRITER’ 
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3 0 0 0 OPEN ' TO 1 ' > OUTPUT » 3 
3002 FLGG-0 

3005 DIM F$<1> ,IKO *T<1> 

3010 CPEM B I PS, INPUT, 4 

302 0 RT ENBFILE<;4> GO TO 30*0 

3030 INPUT FROM 4 1 > ,B< l > n T < 1 > ,BS< 1 ) 

3 035 STP I N6 BON < T < 1 > '< 

3040 IF FS< 1 >=FL1S THEN WRITE ON 3 *D< 1 > : ' . ' :DS<1 > ELSE GO TC 3050 

3045 FLGG-FLGO+ 1 

3 046 BOM < FLGG "> » BS< 1 > 

3 05 0 IF F* < 1 > sPL 1 S THEN GO TO 3 03 0 
3060 N-B< 1 > 

3070 IF TC1>«BU> THEN GG TO 3030 

3080 GO TO 3030 

3090 CLOSE 3 <4 

31 OC BIM D1 ( 1> ,B2S<' 1 > 

3110 STRING RS<N> 

312 0 OPEN ''♦S'TOl " » INPUT >3 
3130 RT ENBFILE<3> GO TO 3170 
3 14 0 I NPUT FROM 3 : B 1 < 1 > , B2S < 1 > 

3150 IF D£*<n«DX* THEN GO TO 3.170 
3160 GO TO 3140 
3170 CLOSE 3 
3180 Q=B 1 < 1 > 

3190 IF LOS* THEN P-0 ELSE GO TO 3210 
3 200 GO TO 3240 

3210 IF LOS- '■ > " THEN P-1 ELSE GO TO 3230 
3220 GO TO 3240 

323 0 IF LOS -CHAR <6 0> THEN P—1 
8240 OPEN ' ♦HOLBl ' » OUTPUT »5 
3250 OPEN FL IS » INPUT »4 
3260 RT ENBFILEC4;- GO TO 3300 
3270 INPUT FPCM 4:R*<I.-' FOP 1 = 1 TO N 
3280 IF CDMP<flS<Q>fDVS^»P THEN ‘WRITE ON 5 USING 
FlS < I ) ! CHRP < NN ;> FOP 1*1 TO N- 1 \RS < N > : CHRP OS'* 

3290 GO TO 3270 
3300 CLOSE 4.5 
3310 RETURNS 

500 0 OPEN ' ♦TPRNSE .OUTPUT ♦ 1 

5U10 WRITE ON 1 :PES i ^ : I : » • :FLGO : ' . ' :BOM< I '< FOP 1*1 TO FL6C 

5020 CLOSE 1 

5030 APPEND ''♦TRRNSF'' TO BIPS 
5040 RETURNS 




LDKD 'WRITPR'' 

>RUN 

ENTER NKNE OF HEW RROGRF IN? PNEW 
SELECT K RELKT IONHL DRERKT I ON <REST >PROJ * JO IN > 
SOURCE RELATION* RNS^EQF'T .KSSIGNHENT 
" DOMAIN NKNE? OWNER 
L06ICKL OPERATOR: = 

DONKIN VhLUE ? 12 
PROJECT TO ? TERM 
PROJECT ALL DONKINS 7 N 
ENTER CUKHTITY OF DONKINS? 3 
FNTPP OF DONKINS? CON »OPi TH jLOCK I 

Ie^IctTrELRT I SSfiL OP«W I OH '■ PEST .FFOJ . JD IH ) 
TEST DNEW 


DR ENTER "STOP“? REST 


OR 


ENTER "STOP"? STOP 


END OF PUN 


Fig. 5. Illustration of a session at a terminal using 'WRITPR' to create an application program 
(User responses are to the right of each colon.) 



Fig. 6. simplified flow diagram of ‘WRITPR' 
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COPY 'BNEW TO TERM I NHL 

100 STRING CONI »0WNER1 *LCCAT1 *0PSTA1 »RECDT1 
110 DPEN 'RMS*€©PT. ASSIGNMENT', INPUT, 1 
120 AT ENDF ILE< 1 ;■ GU TO 1G0 

130 INPUT FROM 1 :C0N1 » OWNER 1 »LOC AT 1 ,DPSTA 1 ,PECDT 1 

140 IF OWNER 1*' 12' THEN PRINT CONI »0PSTA1 fLOCATl 

150 GO TO 130 

160 CLOSE I FOP 1=1 TO 1 

170 STOP CHAR< 13>+- END OF RUN 


Fig. 7. Data dependent application program produced by using 'WRITPR' 


170 LINE -LINE* 10 

ISO WRITE ON 1 :STR<L INE> : ' OPEN ' : '<' v ;SR1*< JJ> s '<, ' >' : ' , INPUT *' i i 
Jl90 LINE-LINE+10 

3200 WRITE ON 1 *STR<t INE> s ' AT ENDFILE<. ' :JJ : ■' > GO TO 
LINE+30+1 0“R 


Fig. 8. Sample code lines from 'WRITPR* (lines 110 and 120 of Fig. 7) 


3210 LINE«LINE<MO 

3220 WRITE ON 1 USING :STR<L INE >* ' INPUT FROM 

NDOMCn :STR<JJ> :CHAR<44> FOR 1 = 1 TO N-l\DOMt:N>: 
STRC J.J> :CHAR < 13 > 


Fig. 9. Sample code lines from 'WRITPR* (line 130 of Fig. 7) 


3250 LINE “LINE 0 

32G0 WRITE ON I s STRC LINE > s' GO TO ' :LINE-< 1 0*1 G*R) 




Fig. 10. Sample code lines from 'WRITPR* (line 150 of Fig. 7) 
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LORI) 'BNEU' 
>RUN 

BB5C11 OL 

12 

CR6B12 US 

1Y 

END OF RUN 


> 


LORD 'INEU' 
>RUN 

BB5C11 OL 12 
CR6B12 US 1Y 


END OF RUN 

> 



Fig. 11. 'lata printout resulting 
fror* running ‘WRITPR* 
and 'WRITER' 


>CDPY 'RHS^EQPT.R I6NMENT / TO TERNINRL 

RR3R12 » 1Y * 1Y >0L » 01 0675 

RB6C21 * IX * IX *DL> 170273 

CD3B15 » IX > IX *SP f 121 072 

BB5C12 * 14* 14 »0L» 091274 

RX3B09-.il tlY>ER> 190176 

DR4C 1 £ » 1 X » 1 X » 0L > 07 1 07 1 

BB5C1 1 >12 » 12 iOL » 091274 

RR7U15 > 1 Y » IX >US » 250276 

BB5C13 >11 >11 >SP> 091274 

CC7C 02 1 1 X > 1 X >US * 15 1 1 72 

CX5B13 > IY > 1 Y »0L *151071 

CR6B12 * 12 »1Y*US >110376 


Fig. 12. Structure and content of the relation (data file) 
used in this article 
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Differential Range Validation: A New Technique for 
Near- Real -Time Validation of Multistation 
Ranging System Data 

A. L. Berman 

DSN Network Operations Section 


Sear-real-time validation of ranging system data is currently restricted to 
multiple range acquisitions during single passes (Pseudo-DRVID). This article 
describes a new technique (“Differential Range Validation’) which utilizes pre- 
dicted range and doppler pseudo-residuals to validate two-station, contiguous-pass 
range acquisitions down to the 10-meter level. 


I. introduction 

During the first part of this decade, near-real-time vali- 
dation of ranging data was only sporadically attempted 
and far less often successful. In March of 1975, this 
author introduced the “Pseudo-DRVID” Technique (see 
Ref. 1), which gave the Deep Space Network (DSN) the 
first viable near-real-time ranging validation capability. 
Since then, “Pseudo-DRVID” has enjoyed considerable 
success in validating multiple range acquisitions during 
single passes. However, there still remained a persistent 
and unfulfilled desire to be able to validate range acquisi- 
tions between two or more Deep Space Stations (DSSs). 
Responding to tins need, this report presents a technique 
to validate two-station, contiguous-pass range acquisi- 
tions — hereafter to be referred to as “Differential Rang. 
Validation 
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II. The Differential Concept 

The reasons why Planetar)' Ranging Assembly (PRA) 
range acquisitions cannot be simply compared to predic- 
tions 1 are thoroughly explored in Ref, 1. Suffice it to say 
here that the range ambiguity values are (or should be!) 
frugally chosen to be only slightly larger than prediction 
(or orbit determination) errors, and this straightforwardly 
causes a direct comparison of PRA measurements to pre- 
dictions to be without meaning in the conventional sense 
of a “residual.” No matter that the absolute range error 
in predictions may be kilometers or tens of kilometers — it 

i Here considered to be output from the PREDIK program. PREDIK 
is the Network Operations Control Center (XOCC) Sigma 5 track- 
ing prediction program. Inputs to PREDIK me the Simulation 
Output Program (oOP) and Fast Pin-Factor Generation Program 
( FPCP). 
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is, in general, true that tlu* [•roteth of predicted-range 
error (excepting orbital and encounter phases!) is fre- 
quently of tlu* magnitude*: 


-0.1 Hz 

This would map into an increase in predicted-range error 
over (say) five* hours of: 


0.1 cycle f .1 motor 
second [_ 15.3 cycles 


118,000 


seconds] 


.1 18 meters 


At the same time, the absolute range change over tlu 
same 5-hour period may easily be 100,000 km or more. 
This then motivates tlu* central concept of the Differ- 
ential Range Validation Technique: 

FRA (or Mu II for that matter) range measure- 
ments between two different DSSs can be 
validated with a high degree of confidence 
by comparing tlu* differenced range acquisi- 
tions to the differenced predicted range. 

Even better, one is not constrained by the accuracy in 
predicted-range error change over several hours. In near- 
real-time, one is automatically given via the Xetwork 
Operations Control Center (\OCC) Pseudo Residual 
Program a frequent measurement of the rate of growth 
of the predicted -range error — the d op pier pxuedo- 
rcsidual. Obv iously, one can easily and substantial l\ cor- 
rect differenced predicted range by simply adding a term 
computed directly from the observed doppler pseudo- 
residual. Hu* specifies of the Differential Range Valida- 
tion technique are presented in the following section. 


Hi. The Differential Range Validation 
Algorithm 

One starts by noting the relationship between range 
and the output of the PR A: 

mi) K[M(t)} ■ RPRA(I) 

0 < RPRA(t) < K 


where 


/?{ /) round trip range at time I 

K - ambiguity resolution factor: a quantized 
input in units of R(f) 


M(t) integer, determined from independent or- 
bital knowledge 

RPRA{t) — “scaled” output of the PRA, i.c., in the same 
units as ft(/) 

Additionally, the following parameters arc required for 
a complete description of the algorithm: 

Xj = parameter X applicable to DSS, 

TSFi = track synthesizer frequency, Hz 
N-, ~ number of components 

K, = K,(TSF>,N.) - ■ 1 3^p ~ 

= ambiguity 

C — speed of light, m/s 

R< t (t) — actual round trio range, m 

= predicted round trip range, m 

AMO = [R,i(0 - lln(t) modulo K,]/Ki 

= [R,,i(t) - 2tyj(/) modulo KiJ/K) 

PRTR i -■ PRA range measurement, range units 
(RU) 

RPRAi(t) - [ PRTR , ] , m 

TO PRA acquisition time 
Biasj — station range bias RU 

Now consider a PRA acquisition at DSS 1 with a 
TO -■ t , and a subsequent PRA acquisition at DSS 2 
with a TO = t.y. 


R, 1 (/ 1 )aK,[AI„(f l )] -i- RPRAdQ 
+ RPRA,(t s ) 

and 

Aft 8 =sR*,(f 8 ) - &„(/,) 

^RPRAM - RPRAM 

+ iUMaM] - K,[AM*.)3} 

Now in Pseudo-DR VI D, one could reasonably assume 
the same number of components and the same track 
synthesizer frequency for each range acquisition during 
the same pass, and hence all relevant quantities would 
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Iifwe the same ambiguity as a modulus. For range acqui- 
sitions at different DSSs, however, one might find differ- 
ent numbers of components: 

and one can certainly expect the track synthesizer fre- 
quencies to be different: 

rSF„ =£ TSF \ 

If the number of components only were different, all 
quantities could easily be operated on by the smaller 
modulus (ambiguity), since 

K, , 

_ _ 2 (A '-*'•) = integer, Nj > N t 
Howeecr, even minor changes in TSF have a dramatic 

£.i“^p,r l m n — 

Ni = N,~ 10 components 
Am — fl«j = 3 X 10 11 m 


the smaller ambiguity (K,). This is accomplished by 
writing: 1 7 

^[A/ <J2 (^)j = K,L + e; L ~ integer 


or 

6 = (F 2 [A/ M (tu)]) modulo Ki 
follow*" "° w ° perate on AR “ ' vith ** ambi S uit y K. 

AR, as RPRA,(t.) - RPRAJt,) + e 
+ (K,L - 


as 


and 


with 


A RPR A (AH,) modulo 

= RPRA s (td - RPRA,(t,) + e 

6 “ (R;[AR;(f^)]) modulo iCj 


Now let 


so that 


and 


Unfortunately, one does not have A# a8 (f .)• however if 
assumes 


one 


TSF, = 22000000 Hz 


K., — 297684.8679 m 


RPH^ - R n modulo K, ~ 36905 


then 




~ (fCo[A/ ; j 2 (f 2 )]) modulo K, 


m 


One can easily see that even if predictions were in error 
by several times the ambiguity 


Similarly, 


TSF, = 22000010 Hz 
K, ~ 297684.7325 m 
RPRA, - R a , modulo K, = 173289 m 

with a difference in PRA measurements of: 

ARPRA = RPRA 2 - RPRA, = 136384 m 


To compensate for the different ambiguities, the ran 
acquisition performed with the larger ambiguity (say h 
* transformed to an “equivalent” range acquisition 
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Mazih/ — M p ,(t ,) ~ /; / = small integer 

the error in e would still be extremely small: 
assume A TSF ~ TSF, ~ TSF, 

I ~ arbitrary integer 
€ + Ae = (K 2 [Af rt2 (t 2 ) -/]) modulo K, 

Ac = (K 2 A/ fl2 (t 2 ) - jk,) modulo K, 

~ { [KaA/«a(f a )] modulo K, - e] 

+ modulo K, + IK, 

“ (~JKs) modulo K, + IK , 
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so that 


- (— /K. 1 modulo K, + IK, 

= (— /K,2'- v “ Yi) ) modulo K, 

+ (/K, 2<»«-»*> [-^ 7 ]) modulo K, + IK, 

= ^/K 1 2 ( - v *" Vl) [^jr]) modulo K, + IK, 


but since 


ATSF 

JK&x^'O—- <«K i 


Ac = £15^ + / Kl 


(Ac) modulo Ki = /K, 2 (AViVl> 


ATSF 


Finally, one obtains the difference in predicted range: 
AR P Rp2(ta) — Rpi(ii) 


A RPRA,, — (A Rp) modulo K t 

= (Rpn(tu) — Rpi^)) modulo K, 

One now has the differential quantity of interest: 

A RPRA - A RPRA„ = RFRAo(tu) - RPRA,(«,) *i* e 

— (R p n(tn) ~ R,n(t\)) modulo K, 


€ = (K., [M pi (t 2 )]) modulo K, 


A RPRA - A RPRA,, = {RPRAn(t 2 ) ~ RPRA^f,) 

+ K.[M*(fe)] + R in (t ,) 

“ R/^fa)} modulo 

IV. Correction of Predicted Range via Use 
of the Doppler Pseudo-Residual 

As mentioned in Section II, the accuracy of the range 
prediction used in differential range validation can be 
substantially improved by merely utilizing the already 
automatically provided doppler pseudo-residual. To 
facilitate the discussion, define the notation 

&X(t) = X(t 2 ) — X(#j) 
SX(f)=X fl (f)-X,,(f) 


where 


Thus, 


= "actual” quantity at time t 
X„(t) — “predicted” quantity at time t 

SR(t) = R«(f)-R„(t) 

= predicted-range error at time t 


Let one now assume a prediction error over a short time 
period (several hours) of the form: 


and hence: 


8 R(t) zz A + Bt -i- Ct‘ 


777 [8R(0] « B + 2 Ct 


D 2 (t) = 96— TSF„ - 96|^TSF,(l ~ ~~)+ Bias,, 


where 


132(f) = two-way doppler 
TSF fc — received track synthesizer frequency 
TSF r = transmitted track synthesizer frequency 
Bias, i = doppler bias 
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so that (with TSF r === TSF ) 


in a fonn as follows: 


SD2 (t) = 



240 /dR, dRA 
221 



ARJ, — A R,, 4- ASR 
“ Rp\i(t •>) R;,l(tl) 


"c 221 1 

(SD2(t,) + 8D2(t..)'\ "j 

96 240 TSF 

2 j] 


= value of doppler pseudo-residual 

Thus one has 

SD2(( a ) = (MrSF(B + 2Ct a } 

( Qfi\ 940 

T)i2r rsF{B + 2C **> 


V. Final Expression for Differential 
Range Validation 

Incorporating the station range delays and the dop- 
pler pseudo-residual correction, one arrives at the final 
expression: 

A RPRA - aRFRA /( - 


and 


(SD2(fi) + 8D2(t,)} 

/ 96 \ 240 

= ) oof ^SF{B + 2C*i + B + 2Ct 2 } 

/ 96 \ 240 

= V“rj22r TSF{2(35 ' i C[ ' 1 + *‘ 3) > 


Now the quantity one is interested in is the range error 
change, A SR, from ^ to £»: 

A8R - 8R(tJ - 8R(ty) 


= A + Bt, + Ctl - (A + Bt , + C*f) 


= B(t a - £,) + C(f2 - £?) 

= (f 2 - f,) [B + C(f s 4- f-J] 

* /221\ 1 fSD2(tO + 8D2(f 3 ) 
l - L 96 \24°y TSF \ 2 


}] 


RPRA,(t,) - RPRA,(£,) + K, 
c f Bias, Bias, 1 


48 1_TSF 2 TSF, J 


A Rjn(ti) R;i 2 (^ 2 ) 




" c 221 /1\ (8D2(f,) 8D2(f..))l| 

_ 96 240 \ 2 / ( TSF, + TSF,, /J) 


modulo K, 


VI. Preliminary Results of Differential 
Range Validation 

Twelve cases of Viking two-station, contiguous-pass 
range acquisitions were compared with the differential 
range validation technique; the results are presented in 
Table I. The twelve cases presented in Table 1 produced 
the following composite result: 

I A RPRA - aRFRA;, I avB - 9.5 m 

It is noteworthy that these results were obtained via the 
exclusive use of routine tracking predictions. 


which is to say that one can incorporate almost exactly 
predicted-range error growth up to second order: 

8R(t) = R,(f) - R„(t) 

~ A Bt t Ct- 

by simply using the observed doppler pseudo residuals 
at t , and t,\ 

SD2(f,), 8D2(£A 


An HP9810 program containing the differential range 
validation algorithm (see Ref. 2) lias been delivered to 
die Network Analysis Team, Tracking (NAT Track), and 
the technique is considered operational. 

VII. Summary 

Attempts to validate ranging system data in near-real- 
time prior to 1975 were generally unsuccessful. In 1975, 
the Pseudo-DRVID technique was introduced, and it 
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proved quite successful in the near-real-time validation 
of multiple range acquisitions during single passes. How- 
ever, there still existed a need to be able to validate 
ranging data between separate DSSs. The differential 


range validation technique presented in this report 
answers this need by providing a practical method of 
validating two-station, contiguous-pass range acquisitions 
down to the 10-meter level during cruise phases. 
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The 1976 Helios and Pioneer Solar Conjunctions- 
Continuing Corroboration of the Link Between 
Doppler Noise and Integrated Signal 
Path Electron Density 

A. L. Berman, J. A. Wackley, and S.T. Rockwell 

DSN Network Operations Section 


Observed doppler noise (rtns phase jitter) from the 1976 solar conjunctions of 
the Helios 1 and 2 and the Pioneer 10 and 11 spacecraft was processed with a 
recently developed doppler noise model “ISEDB. Good agreement is obtained 
between the observed data and the model . Correlation is shown between devia- 
tions from the ISEDB model and sunspot activity , but it is insufficient to be 
modeled . Correlation is also shown between ISEDB model deviations for (space- 
craft) signal paths on the same side of the sun . 


i. Introduction 

Utilizing observed doppler noise data from the 1975 
Helios 1, Pioneer 10 and Pioneer 11 Solar Conjunctions, 
the authors have constructed a geometric model (the 
“ISED” series) for plasma-induced doppler noise (rms 
phase jitter). The functional form of the model was de- 
veloped by assuming that doppler noise was proportional 
to the integrated signal path electron density, with the 
electron density assumed to be: 



r — heliocentric distance 


Reference 1 describes a modification to the model 
(“ISEDB”) which accounts for heliocentric latitude; Refs. 
2 through 5 document the development of the basic 
model. The most recent model — ISEDB — is defined as 
follows : 

ISEDB = 

KMWM ‘ A ' 10 y 

+ (0.0015) 2 T* 
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where 


F(a,/?) = 1 - 0.05 i— — t g) * ~ ( g ~ Zl^L 


( 


: } 


- 0.00275 f (/? ~ t/2 + g) ° ~ ( g - z!R 

l p 

a = Sun-Earth-Probe angle (SEP), rad 
/? = Earth-Sun-Probe angle (ESP), rad 


and 


~ heliographic latitude, deg 
~ sin-' [cot a(- cos sin a w sin e + sin 8 d cose)] 
a r ,i = right ascension 
8,i = declination 

e = the obliquity of the ecliptic (23.445 deg) 

with 


A 0 = 9.65 X 10“ 
A, = 5 X 10-“' 
A s - 9 X 10-' 


doppler sample rates of either 10 or 1 second. In these 
cases, the following relationship was used to transform 
the doppler noise from the actual sample rate to an 
“equivalent” 60-s sample rate: 

/ SR2 Y'-- 85 

noises , { , = noises, t-A Sfli J 


where 


SRI, 2 = sample rate 1, 2 

The collected noise data was fit with the ISEDB model 
and a new A„ was selected to minimize the standard 
deviation of the residuals: 

A 0 = 11.21 X io-> 

Scatter diagrams for the data (observed doppler noise 
vs the ISEDB model) are seen as follows: 

Fig. 1 — Pioneer 10 
Fig. 2 — Helios 1 
Fig. 3 — Helios 2 
Fig. 4 — Combined 


It will be the intent of this article to ascertain if the 
Helios 1 and 2 and the Pioneer 10 and 11 observed 
doppler noise continues to be well represented bv the 
ISEDB model. 


II. Data Description 

The data consists of pass average, good, two-way 
doppler noise collected during the following time inter- 
vals (DOY = Day of Year, 1976): 


Spacecraft 

Begin (DOY) 

End (DOY) 

Pioneer 10 

51 

175 

Pioneer 11 

79 

215 

Helios 1 

75 

293 

Helios 2 

90 

293 


All doppler data was taken at a 60-scc sample rate with 
die exception of some Helios 2 passes, which utilized 


while the following figures present observed doppler 
noise and the ISEDB model as a function of DOY: 

Fig. 5 — Pioneer 10 
Fig. 6 — Pioneer 11 
Fig. 7 — Helios 1 
Fig. 8 — Helios 2 

Examination of these figures leads one to conclude in a 
qualitative sense that the (ISEDB) model continues to be 
in good agreement with the doppler noise observations 
during the 1976 Pioneer 10 and 11 and Helios 1 and 9 
solar conjunctions. The statistics of the observed doppler 
noise as fit to the ISEDB model were as follows with 
dB - 10 Iog, 0 (observed noise/ISEDB): 



Pioneer 

10 

Pioneer 

11 

Helios 

1 

Helios 

2 

All 

o - , dB 

2.04 

1.77 

2.47 

2.65 

2,28 

Bias, dB 

— 0.37 

-0.01 

-0.63 

11.36 

-1-0.02 
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HI. Correlation of Observed Doppler Noise 
With Sunspot Activity 

In Ref. 2, the authors attempted to correlate observed 
doppler noise witli sunspot activity (the “RISED” model). 
Although a quantitative correlation with sunspot activity 
(as measured by R* — Zurich (daily) sunspot index) 
was demonstrated, it was not (at that time) considered 
optimistically. The data presented in this report was 
processed with the RISED model (now “RISEDB ), and 
a new A„ to minimize the standard deviation was com- 
puted to be 


A {> - 10.16 X 10- ' 

Basically, the RISEDB performance (as compared to 
the ISEDB modeling) can be summarized by spacecraft 
as follows: 

Pioneer 10 — significant degradation 
Pioneer 11 — minor improvement 
Helios 1 — significant improvement 

Helios 2 — -unchanged 
All — minor improvement 

The statistics are as follows: 



Pioneer 

10 

Pioneer 

11 

Helios 

1 

Helios 

2 

All 

a, dB 

2.25 

1.72 

2.20 

2.65 

2.22 

Bias, dB 

-0.43 

+ 0.30 

-0.53 

4- 1.06 

4- 0.05 


ISEDB residuals, “smoothed/phased” sunspots (--XR-; 
see Ref. 2), and RISEDB residuals are seen in the fol- 
lowing figures: 


Fig. 9 — Pioneer 11 


Fig. 10— Helios 1 
Fig. 11 — Helios 2 


Regions which appear to show evidence of correlation 
between observed doppler noise and sunspot activity 
are as follows: 


Pioneer 11, DOY S0-145 (Fig. 9) 

Helios 1. DOY 75-175 ( Fig. 10) 

Helios 2, DOY 90-135 (Fig. 11) 

However, these regions must be balanced against regions 
of little or even negative correlation, as follows: 

Pioneer 11, DOY 175-195 (Fig. 9) 

Helios 1, DOY 1S5-255 (Fig. 10) 

Helios 2, DOY 140-170 (Fig. 11) 


IV. Multispacecraft Correlation 

One would expect to see good correlation between 
ISEDB residuals from different spacecraft, as long as the 
signal paths are on the same side of the sun. Figure 12 
shows ISEDB residuals for Pioneer 10, Pioneer 11, and 
Helios 2 when these signal paths were east (left) of the 
sun, while Fig. 13 shows ISEDB residuals for Pioneer 10, 
Pioneer 11, Helios 1, and Helios 2 when these signal 
paths were west (right) of the sun. Regions of correlation 
between various spacecraft are as follows: 

HE 1/HE 2, DOY 90-134 (Fig. 13) 

PN 11/HE 1, DOY ISO-215 (Fig. 13) 

PN 10/PN 11, DOY 80-105 (Fig. 12) 


V. Summary 

Observed doppler noise data from the 1976 solar con- 
junction phases of the Pioneer 10 and 11 and the Helios l 
and 2 spacecraft was fit to the ISED model (“ISEDB”) 
and subsequently shown to he* in good agreement with 
the model — thus continuing to corroborate the direct link 
between observed doppler noise and integrated signal 
path electron density. Some evidence of correlation be- 
tween sunspot activity and doppler noise is seen (the 
RISED model), but the relationship is insufficient to 
model in a quantitative fashion. Finally, correlation 
between doppler noise deviations from the ISEDB model 
continue to he seen for (spacecraft) signal paths on the 
same side of the sun. 
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Fig. 4. Early 1976 solar conjunction mission composite, observed doppler noise vs the ISEDB model 
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Fig. 5. Pioneer 10 1976 conjunction, observed doppler noise and the ISEDB model vs day of year 
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DAY OF YEAR, 1976 

Ffg. 6. Pioneer 11 1976 solar conjunction, observed doppler noise and the ISEDB model vs day of year 
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Fig. 7. Helios 1 1976 solar conjunction, observed doppler noise and the ISEDB model vs day of year 
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Fig. 9. Pioneer 11 1976 solar conjunction, ISEDB residuals, RISEDB residuals and (smoothed/phased) sunspots vs day of year 
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DSS Range Delay Calibrations: Current 
Performance Level 

G. L. Spradlin 

DSN Network Operations Section 


U is the intent of this report to describe a task undertaken approximately 18 
months ago, under the auspices of the flanging Accuracy Team, to develop a 
means for evaluating Deep Space Station (DSS) range delay calibration 
performance, and through which inconsistencies frequently noted in these data 
could be resolved. Development of the DSS range delay data base is described. 
The data base is presented with comments regarding apparent discontinuities. 
Data regarding the exciter frequency dependence of the delay values are 
presented, and the report concludes by noting the improvement observed in the 
consistency of current DSS range delay calibration data over the performance 
previously observed. 


I. Introduction 

1 he measurement of the distance, or range, to a 
spacecraft, as accomplished by the DSN Ranging Subsys- 
tem, is essentially a precision measurement of the time 
interval between transmission and receipt of a binary- 
coded signal called the range code. Unfortunately, the 
time interval (from transmission to receipt of the range 
code) measured by the Ranging Subsystem is not the 
desired measurement. For the purposes of orbit determi- 
nation, the measurement needed must he from the defined 
station location to a reference point on the spacecraft and 
back, and therefore, the delay actually measured bv the 
Ranging Subsystem must be modified (corrected) to 
account for delays of the range code that are not part of 
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this refercnce-point-to-refercnee-point measurement. 
Among other correction factors that must be applied to 
the measured value to obtain the desired value is the delay 
experienced by the range code as it passes through the 
station equipment, or the station range delay. 

Historically, the station range delay has been considered 
one of the largest contributors to inaccuracies in 
spacecraft range measurements despite the measurement 
(calibration) of the station's range delay prior to and 
immediately after each spacecraft track during which 
range data were taken. Inconsistencies and apparent errors 
in the station range delay calibration could be expected to 
contribute at least several meters of uncertainty to the 
spacecraft range measurement. The inconsistencies of 
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station range delay calibrations were recognized by the 
Ranging Accuracy Team as lieing incompatible with the 
stringent range accuracy requirements of future missions, 
and hence, one of the initial tasks undertaken by the team 
was to develop a data base from which current perform- 
ance could be ascertained, and improvement monitored. 
This report describes this effort, and delineates the current 
station range delay calibration performance. 

II. Establishing The Data Base 

It was decided in early June 1975 that the data base 
would not include any currently available data, nor would 
it include station range delay calibrations as reported via 
the station posttrack reports. Instead, each station capable 
of ranging was requested to perform special range delay 
calibration measurements. 

To provide control over the conditions under which the 
measurements were to be made, each station was 
requested to use a specific set of parameters including 
specifications of range code component integration times, 
level of carrier suppression, downlink ;.ignal level, antenna 
pointing angles, and uplink power level. The report trom 
the station was to include, in addition to the range delay 
calibration, five postrange acquisition Differenced Range 
Versus Integrated Doppler (DRVID) values so that the 
range delay calibration could be corrected for waveform 
distortion, exciter frequency, the number of the receiver 
used, and which maser was used. 

It was originally thought that a one-month period of 
data collection would provide a sufficient data base; 
however, due to a lack of available time at the stations to 
perform unscheduled activities, far too few measurements 
were made available to constitute a data base. Time for 
the stations to take the needed measurements was then 
made available through the Deep Space Network (DSN) 
Scheduling Office. The measurements were to be made on 
a noninterference basis during or immediately after 
normal station posttrack countdown procedures. One-half 
hour was added to the allotted time tor the normal 
posttrack procedures to help alleviate any time bind that 
might be caused by making the additional measurements. 
Measurements were scheduled to be made at each station 
as frequently as once per day with no more than seven 
days between measurements, depending upon station time 
availability, rdso, the time span over which data would be 
collected extended from the original on^ month 

period to a six-month period. 

Coincident with the effort to cstabli'n a data base, a 
task was undertaken to reevaluate and update the standard 


pre- and posttrack range delay •alibration procedures, 
including the reporting fo/mat (Post Track Report). 
Provisions were made in the new procedures to incorpo- 
rate the controlled range delay calibration into the 
standard station procedures. The Post Track Report 
format was redesigned to provide detailed information 
regarding configuration, and to provide the needed data to 
ascertain the health of the ranging subsystem. The obvious 
advantage of the aforementioned incorporation is that the 
data base would now receive almost continuous inputs as a 
product of normal station operations, and station range 
delay calibration performance could l>e more closely 
monitored. 

After a confidence-building period of approximately one 
month, the specially scheduled calibrations were discon- 
tinued, and the new Post Track Report was adopted as the 
data source for the data base. A copy of the portion of the 
new Post Track Report format containing range calibra- 
tion information is shown as Fig. 1. Item A.3 contains the ' 
station range delay calibration data, and the necessary 
configuration information. The information used as input 
to the data base is found in items 3.(a).l and .2, and 3.(b). 1 
and .2. The numbers reported here by the stations are 
converted from the frequency dependent unit of measure- 
ment (Range Unit, RU) output by the ranging subsystem, 
after adjustment by the DRVID correction, to meters 
(one-way) for ease of comparison with flight project 
navigation and radio science requirements, and for entry 
into the data base. 

III. The Data Base 

Station range delay calibration data extracted from the 
data base and plotted in meters (one-way) versus day of 
year are presented in Figs. 2 through 10. Specific station 
performance will not be detailed at this point. Instead, it 
is intended simply to present the data base, and to explain 
peculiarities, events, and characteristics that can be seen 
in the plotted data. A more detailed analysis of the data 
will be presented in a subsequent section of this report. 

A. Gotdstone, Pioneer, DSS 11 

DSS 11 is the station to have most recently acquired a 
ranging capability, this event occurring approximately 
sixty days after data collection began from other stations. 
Early data from DSS 11 indicated that, for whatever 
reason, this station did not reflect the inconsistency of 
range calibrations noted in the early data from other 
stations. Jn fact, it l)ecame an interim goal to try to 
improve the performance of other stations to that level 
already being achieved at DSS 11. As indicated on this 
station’s plot, an Engineering Change Order (ECO) that 


JPL DEEP SPACE NETWORK PROGRESS REPORT 42-36 


139 


altered the station range delay calibration path was 
installed temporarily, day 63 through day 70, and 
permanently on day 158. All calibrations indicated on this 
plot were made using Maser 1 and Receiver 1. 

B. Goldstone, Mars, DSS 14 

Figure 3 is a plot of DSS 14 S-band calibrations. All data 
were taken using the S-Band Polarization Diversity (SPD) 
maser and primarily Receiver 3 of the Block IV Receiver- 
Exciter Subsystem, although there is a scattering of points 
taken using Receiver 4. Also, all calibrations were made 
using the Block IV doppler translator instead of a dish 
mounted Zero Delay Device (ZDD); both devices perform 
the function of simulating a spacecraft by converting a 
transmit-level frequency to a receive-tevel frequency. 

Prior to day 90, DSS 14, like other stations in the 
Network, utilized the Planetary Ranging Assembly (PRA) 
to conduct spacecraft ranging. On day 90 an R&D ranging 
machine, the MU II, was installed to provide range data 
enhancement for the Helios project at superior conjunc- 
tion. The MU II has been retained at DSS 14, and is 
currently used for both Viking and Helios ranging 
operations. 

Figure 4 presents DSS 14 X-band range calibrations. 
The family of points seven to eight meters below the 
major point grouping represents a periodically occurring 
anomaly that recently has been determined to be the 
result of a procedural error (incorrect switch position) 
during the range delay measurement. 

C. Australia, Weemala, DSS 42 

DSS 42 range calibrations presented in Fig. 5 were 
taken using Receiver 5, Maser 1, and a dish-mounted 
ZDD. 

As mentioned earlier, DSS 11 did not indicate the 
inconsistency in range calibrations noted elsewhere in the 
Network. As data collection progressed, it became quite 
clear that DSS 42 was unable to provide consistent 
calibrations, and that the problem was not due to an 
equipment failure, nor was it procedural in nature. It also 
became clear that DSS 42 range calibrations were highly 
exciter reference frequency dependent, and that the 
calibrations at DSS 11 were not. Further, it was found that 
the locations of the dish-mounted ZDDs at these two 
stations were radically different, and that a poor location 
of the ZDD could result in the calibration inconsistency 
being observed, the apparent result of signal multipath 
effects. A pragmatic decision was made to relocate the 
DSS 42 ZDD to a position on the antenna identical to that 
where the ZDD was located on the DSS 11 antenna. This 
was done on day 114 with an obvious and dramatic 
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reduction in range calibration inconsistency. An additional 
small step in the plotted calibrations can be seen to occur 
on day 140 coincident with a klystron change. 

D. Australia, Ballima, DSS 43 

DSS 43 calibrations presented in Fig. 6 reflect a 
multiplicity of configurations. The majority of data 
collected prior to day 30 of 1976 were taken using the 
Block III Receiver-Exciter Subsystem, Receiver 1, the 
SPD Maser, and a dish-mounted ZDD. After day 30, the 
Block IV Receiver-Exciter Subsystem became the prime 
data source, and data from Receivers 3 and 4 in 
conjunction with the SPD maser are represented. As with 
DSS 14, all Block IV measurements utilized the Block IV 
doppler translator. 

Calibrations made at DSS 43, while using the dish- 
mounted ZDD, again indicate a strong exciter-frequency 
dependence. The apparent increase in the calibration 
measurement after day 260 of 1975 resulted from work 
done in the tricone area of the antenna, but did not 
appear to significantly alter the apparent frequency 
dependence of the dish-mounted ZDD calibrations. 

The cause of the dramatic improvement indicated after 
day 140 is not clearly understood, but it is felt to be the 
result of possibly two things: (1) stringent configuration 
control exercised per Viking project requirements, and (2) 
a significant decrease in the frequency dependence of 
range calibrations made using the doppler translator 
(possibly due to an inadvertent but fortunate elimination 
of a multipath or leakage problem). 

Figure 7 presents DSS 43 X-band range calibration data. 
Again a distinct change in the plot occurs at day 140. 
With the exception of the anomalous calibrations reported 
on day 142 and day 162, all calibrations after day 140 fall 
into one of two distinct families of points. Tests conducted 
at DSS 43 indicate that while there is only a minor change 
in the S-band range calibration value as exciter frequency 
is changed, the change in the X-band calibration value is 
quite significant Thn* iiio two families are not unex- 
pected, as the upper grouping of points represents 
calibrations extracted from Viking Orbiter 1 Post Track 
Reports (i.e., the calibrations were made in the region of 
frequency channel nine), and the lower point grouping 
from Viking Orbiter 2 Post Track Reports (calibrations 
made in the region of frequency channel 20). 

E. Spain, Robledo, DSS 61 

DSS 61 range delay calibrations shown in Fig. 8 were 
made using Receiver 5, Maser 1, and a dish-mounted 
ZDD. 
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The small grouping of points between days 261 and 282 
of 1975 at approximately the 360-meter level reflects the 
temporary insertion of an additional cable in the range 
delay calibration path during station troubleshooting of a 
ranging subsystem problem. There is no available 
explanation for a similar family of points between days 43 
and 76 of 1976. 

Although it is not obvious by inspection of DSS 61 
range delay calibrations presented on Fig. 8, tests 
conducted at DSS 61 indicated that this station’s range 
calibration value was also exciter-frequency dependent 
and, like DSS 42, the DSS 61 ZDD was relocated to a 
position similar to the DSS 11 ZDD location on that 
antenna. The ZDD relocation occurred on day 77 of 1976, 
after which a clear improvement was observed in the 
range calibration consistency. An additional discontinuity 
and noted increase in frequency dependence is observed 
on day 139, coincident with an S-band mixer change. 

F. Spain, Robledo, DSS 63 

Like DSS 43, DSS 63 range calibrations presented in 
Fig. 9 indicate that a number of configurations were used 
for calibrations at this station through day 125 of 1976. 
Data prior to day 300 of 1975 reflect use of the Block III 
Receiver 1, the SPD maser, and a dish-mounted ZDD, 
Measurements after day 4 and before day 112 reflect 
mixed use of Block IV Receivers 3 and 4, the SPD maser, 
and the Block IV doppler translator. With the exception 
of four points, all measurements after day 112 were made 
using Receiver 3, the SPD maser, and the Block IV 
doppler translator. 

Use of the high -power klystron at DSS 14 and DSS 43 
has not been mentioned, as use of this device has not had a 
significant impact on the reported range delay calibration. 
This statement cannot be made of calibrations reported at 
DSS 63. For an unknown reason the use of the high-power 
klystron at DSS 63 does result in a significantly different 
range delay. The family of points between 648 and 651 
meters reflects its use. 

The increased point-to-point consistency of calibrations 
after day 135 is again attributed to tight station 
configuration control, a result of Viking Project support. 

Figure 10 presents DSS 63 X-band range calibration 
data. The data were taken using primarily Receiver 4; 
however, a scattering of points representing Receiver 3 is 
also present. 

The family of points between roughly 624 meters and 
628 meters is representative of use of the high-power 
klystron (denoted by the “H,” above the plotted point), 
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and an anomaly recently traced to a faulty relay in the 
Block IV X-band doppler translator. 

Again there is a general reduction in the scatter of the 
calibrations after day 130, and as with the S-band 
calibrations, it is felt this phenomenon is attributable to 
configuration control. 


IV. Calibration Performance Analysis 

The frequency dependence of the station range delay 
calibrations has been mentioned numerous times in the 
preceding section of this report. The cause of the 
frequency dependence of the range delay is not clearly 
understood as of this writing. In fact, the extent to which 
it is a problem is also not clear. In an effort to understand 
how the range delay varies as a function of frequency, 
each station capable of ranging has been requested to 
conduct a “Range Delay Versus Frequency Test.” 
Essentially, the test simply provides calibration data at 
channel center frequencies, channels 5 through 26. 

It was previously mentioned that DSS 11 did not reflect 
the frequency dependence noted in data from DSS 42 and 
DSS 61, and that the ZDDs at DSS 42 and DSS 61 were 
relocated to agree positionally with the DSS 11 ZDD 
location, in an apparently successful attempt to reduce the 
DSS 42 and 61 frequency dependencies. Figure 11 
presents the results of a range delay versus frequency test 
conducted at DSS 11 on day 252 of 1975. Clearly, there is 
a variation of the range delay with change in frequency; 
however, the peak*to-peak change is only about 1.5 
meters. Compared to the approximate 8-meter and 
7-meter changes seen on Figs. 12 and 13 (DSS 42 and DSS 
61, respectively) DSS 11 appears relatively insensitive to 
changes of frequency. 

Figures 14 and 15 present DSS 42 and DSS 61 range 
delay versus frequency data after relocation of their 
ZDDs, The reduction in frequency dependence is 
noticeably significant with both stations showing a range 
delay variation with frequency change of roughly 2 to 3 
meters peak-to-peak. 

Tests run at DSS 14 and DSS 43 also indicate a range 
delay variation with frequency change. The results of the 
DSS 14 test are shown in Fig. 16 (S-band) and Fig. 17 (X- 
band) with indicated peak-to-peak changes of approxi- 
mately 6 meters and 3 meters respectively. The DSS 43 
results shown in Fig. 18 (S-band) and Fig. 19 (X-band) 
indicate an S-band peak-to-peak change of about 3 meters, 
and an X-band peak-to-peak change of roughly 6 meters, 
and for an unknown reason, a significantly different profile. 
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Test data from DSS 63 are not yet available; however, 
inspection of range delay information in the data base 
indicates that neither the S-band nor X-band range delay 
calibrations vary si" lificantly as a function of frccjuency. 

Each spacecraft tracked by the DSN is assigned and 
communicates within a particular frequency channel. All 
station range delay calibrations made and reported via the 
Post Traci: Report are measured at a frequency within the 
frequent*' channel assignment of the spacecraft whose 
track is being reported upon. With the information 
presented above (a station’s range delay calibration is 
frequency dependent), the consistency or repeatability of 
a station’s range delay calibration cannot clearly be 
ascertained without consideration of the spacecraft {or 
spacecraft frequency channel assignment) as an indepen- 
dent variable. That is, the determination of' the consist- 
ency of a station’s range delay calibration should be based 
upon a consistent spacecraft (i.e., consistent frequency), as 
much as it would upon use of a consistent station 
hardware configuration (same receiver, same maser, etc.). 

The data presented in the data base consist primarily of 
range delay calibrations extracted from Post Track 
Reports for Viking Orbiter 1 (VO-1) and Viking Orbiter 2 
(VO-2). The information presented in Table 1 describing 
station range delay calibration performance was con- 
structed considering the spacecraft as an independent 
variable. The mean and standard deviation (1 cr) were 
computed using the last N range delay calibrations 
(excluding known anomalous points mentioned earlier) 
reported from each station. The mean for each configura- 
tion has been indicated on the data base plots (Figs. 2 
through 10). Only those station hardware configurations 
recently and frequently used have been considered. 

It was mentioned at the beginning of this report that 
station range delay calibrations have historically been 
considered a major source of inaccuracies in the measure- 
ment of spacecraft range, and that the inconsistencies in 
the station calibrations were suspected to have contrib- 
uted up to ‘-veral meters of uncertainty to the range 
measurement. 


Data presented in Figs. 2 through 10 clearly indicate 
that unexplained point-to-point inconsistencies of from 5 
to 10 meters were commonplace in the calibrations 
recorded during the first several months of data collection. 
The data base also delineates that by relocation of the 
ZDDs at DSS 42 and DSS 61, stringent configuration 
control at all stations, an increasing knowledge of the 
sensitivity of the range delay to many variables, and 
continuous monitoring and feedback regarding calibration 
performance to the stations, the consistency of the 
calibrations has been greatly increased, and point-to-point 
consistency of less than one meter (1 cr) is the current 
level of performance. 


V. Summary 

Station range delay calibrations have long been 
considered one of the major contributors to inconsistencies 
observed in spacecraft range measurements. Recognition 
of this plus knowledge of future mission needs for highly 
accurate range data led to an effort undertaken by the 
Ranging Accuracy Team, to acquire knowledge regarding 
the current ability of the Network to provide consistent 
range delay calibrations, and to develop a data base 
through which improvements in calibration performance 
could be monitored. 

It has been the intent of this report to describe the data 
base development, graphically present the data base, and 
to pro* Id's some explanation of trends or characteristics 
that are observed in the graphical presentation. A 
substantial amount of data is presented regarding one 
easily discernible characteristic-the frequency depen- 
dence of the range delay. 

Table 1 is presented at the conclusion of the analysis 
section, and provides a summary of station range delay 
calibration performance. The mean and standard devia- 
tions for each frequently used station configuration 
(receiver, maser, and spacecraft) are presented. It can be 
seen that station calibrations are now quite consistent, 
with the standard deviations being generally less than one 
meter. This represents a significant improvement over the 
5- to 10-meter variations observed early during the data 
base development. 
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Table 1. Station range delay calibration performance 


DSS 

Coitii guration 

N 

MEAN, m 

a, m 

11 

RCV R 1/TWM 1/ VO- 1 

33 

307.0 

0.4 


RCVR1 /TWM l/VO-2 

62 

306.2 

0.4 

14 

RCVR3/SPD/VO-1 

57 

493.5 

0.5 


RCVR3/SPD/VO-2 

30 

495.6 

0.4 


RCV R4/XRO/VO- 1 

53 

481.7 

0.6 


RCVR4/XRO/VO-2 

23 

482,5 

0.3 

42 

RCVR5/TWM1/VO-1 

8 

381.7 

0.6 


RCVR5/TVVM l/VO-2 

33 

382.1 

0.5 

43 

RCVR3/SPD/VO-1 

43 

509.1 

0.3 


RCV R3/SPD/ VO-2 

22 

508.7 

0.3 


RCVR4/XRO/VO-1 

41 

498.0 

0.7 


RCV R4/XRO/ VO-2 

20 

494.7 

0.6 

61 

RCVRS/TWM 1/VO- 1 

12 

353.5 

0.9 


RCVR5/TWM l/VO-2 

35 

350.9 

0.7 

63 

RCVR3/SPD/VO-1 

36 

656.7 

0.5 


RCVR3/SPD/VO-2 

29 

657.3 

0.5 


RCVR4/XRO/VC-1 

31 

634.1 

0.5 


RCVR4/XRO/VO-2 

29 

635.0 

0.5 
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A. MISSION INDEPENDENT. 

1 • DSS ; DAY ; PASS : MISSION ; S/C ID 

2. ACQUISITION_ _Z ; END OF TRACK Z. 

3. RANGE DELAY: 


(a) S-B AND PRETRACK MEASUREMENTS AT -150 DBM (1, 2) AND AT 
EXPECTED SIGNAL LEVEL (3, 4): 


1 . 

2 . 

4. 


RANGE 
RANGE " 
RANGE 

RANGE 

XMTR PWR 

EXC FREQ 


RU ; DRVID 

RU; S/L -150 

DBM; 

RDA ATTR 

DB 

RU ; DRVID 

RU 




RU ; DRVID 

RU; S/L 

.DBM; 

RDA ATTR 

DB 

RU ; DRVID 

RU 




KW; EXC NO. 

; TWM 

t 

RCVR NO. 



HZ; CARRIER SUPPRESSION (3 and 4) DB 


(b) X-BAND PRETRACK MEASUREMENTS AT -150 DBM (1, 2) AND AT 
EXPECTED SIGNAL LEVEL (3, 4): 


1. 

2. 

RANGE 

RU ; DRVID 

RU; S/L 

-150 DBM; RDA ATTR 

DB 

RAN G E 

RU ; DRVID 

RU 

— 


3. 

4. 

RANGE 

RU ; DRVID 

RU; S/L 

DBM: RDA ATTR 

DB 

RANGE 

RU ; DRVID 

RU 

' 



TWM ; RCVR NO. 


(c) S-BAND POST TRACK MEASUREMENTS AT ACTUAL SIGNAL LEVEL OF THE 
PASS: 


RANGE RU ; DRVID RU; S/L DBM; RDA ATTR 

RANGE RU; DRVID RU * 

XMTR PWR K W ; EXC NO, ; TWM : RCVR NO 

EXC FREQ HZ; CARRIER SUPPRESSION DB 

XMTR PWR DURING PASS KW; KLYS NO. 

(d) X-BAND POST TRACK MEASUREMENTS AT ACTUAL SIGNAL LEVEL OF THE 
PASS: 

^ ANG ? HU ; DRVID RU; S/L DBM: RDA ATTR 

RANGE RU; DRVID RU * 

TWM ; RCVR NO. 

(e) REMARKS: * 

4. STATION TIME OFFSET WITH RESPECT TO DSN MASTER: 

(a) OFFSET_ UNCERTAINTY REFERENCED TO ** 

(b) TIMING SYSTEM RESET SINCE LAST POST TRACK REPORT: 

MAGNITUDE DIRECTION . 

5. SYSTEM NOISE TEMPERATURE DEG. K 

6. EQUIPMENT FAILURES OR ANOMALIES (indicate start and end times or start and 
estimated time for return to operation. Reference TFR(s) and DR(s).) 


Fig. 1. Post Track Report format 
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Fig. 6. Station range delay calibration data, DSS 43 S-band 
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Fig. 7. Station range delay calibration data, DSS 43, X-band 
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Fig. 8. Station range delay calibration data, DSS 61 


.1 

i 



JPL DEEP SPACE NETWORK PROGRESS REPORT 42-36 149 


666f 


664t 


6621 


660f 


658f 


656^ 


654t 


652 \ 


650 r 


648h 


O 


°Q,§> 

8 


car 

cP % 


°°$o 
6 




o~o 

o 


cA 

cfS> 

□ 


CpS 03 ° 


A RCVR 1 
O RCVR 3 
□ RCVR 4 


high-power 

KLYSTRON 


646E_L 


180 


200 


_L 


220 240 260 280 300 320 340 360 


20 


40 


60 


B0 


100 120 140 


160 


180 200 220 240 260 280 300 


1975 


1976 


DAY OF YEAR 


Fig. 9. Station range delay calibration data, DSS 63 S-band 






Fig. 10. Station range delay calibration data, DSS 63 X-band 



Fig. 11. DSS 11 station range delay versus frequency 
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DSN Research and Technology Support 

E. B. Jackson 

Radio Frequency and Microwave Subsystems Section 


1* activitiesof the Venus Station^ IdUnd the 
( MTF ) during the period June 14 through vctoo 
progress noted. 

IpeJon of the Venus Sjajionj * reported* on* u>tth 

extensive ^teldaZ Station (DSS 14), and stability and 

OSS 13 

Z^sZonty, Pulsar notion IZ'Zd 

Zldlo fcZsyZlronixationtransrn^ 

noted as 39 transmissions were made to Australia, DS , P 


l. Development Support Group 

The activities of the Development Support Group, in 
operating the Venus Station (DSS 13) and the Microwave 
Test Facility (MTF) during the period June 14 through 
October 10, 1976, are discussed. 


II. Station Automation 

In support of RTOP 68 (Station and Network Monitor 
and Control Technology Development) and RTOP 69 


•racking Station Systems Technology), DSS 13 will be the 
jmonstration station with which remote^ operated, 
attended tracking operation is demonstrated. 

The first phase, remote controlled observation of pulsars 
as been demonstrated on several occasions. During is 
eriod, 35 station support hours, of which 6 hours were 
jmote controlled tracking, were utilized m ^PP ort ° 
lis project. During the remote c ^ tr oiled Pulsar obse - 
ions, pulsars 0031-07, 0329 + 54, 0355 + 54 0628-28 

736-40, and 0833-45 were successfully observed withou 
ntervention by on-station personnel. The rest of e 
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support time was devoted to troubleshooting and perform- 
ance measurement. 

III. Pulsar Observations 

In support of the Radio Science Experiment “Pulsar 
Rotation Constancy” (OSS 188-41-51-09), DSS 13 provided 
130-3/4 hours of observations, during which the emissions 
from the pulsars tabulated in Table 1 were recorded. 
These data, recorded at 2388 MHz, left circular polariza- 
tion (LCP), are used to determine precise pulse-to-pulse 
spacing, changes in this spacing, pulse shape, and pulse 
power content of the signals received from these pulsars. 

IV. Maser-Receiver-NAR Reliability-Stability 
Testing 

Reliability and stability testing of the DSS 13 total 
receiving system is conducted during nonmanned station 
periods. The 26-m antenna is prepositioned to a fixed 
azimuth and elevation and the Noise Adding Radiometer 
(NAR) data collection system automatically records total 
receiving system temperature as the 28-m antenna beam is 
swept across the sky by the rotation of Earth. During this 
period, the antenna was fixed in azimuth at 360 degrees 
and progressively positioned from 49.9 to 49.1 degrees 
elevation, and 496-/2 hours of testing were automatically 
performed. This testing is performed at 2295 MHz, using 
right circular polarization (RCP) on the 26-m antenna. 

V. Unattended Operation, DSS 13 

In preparation for the planned unattended operation of 
DSS 13, the Polarization Diversity S-Band (PDS) feedcone 
removed from DSS 14 has been modified into a feedcone 
having all the capabilities of an S-Band Polarization 
Diversity (SPD) feedcone, except simultaneous reception 
of two polarizations, but using only two waveguide 
switches instead of three. Two Block III Receiver-Exciters, 
obtained surplus from the Spaceflight Tracking and Data 
Network (STDN), have been installed in the Operations 
and Data Processing building at DSS 13 and completely 
checked out with the STDN and DSN frequency 
converters. The first “milestone” demonstration is planned 
to be an automated telemetry reception track from a 
spacecraft in March, 1977. 

VI. 26-m Antenna Measurements and 
Evaluation 

In order to evaluate the effects of several * sars of 
experimental operation on the stability of the surface 
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panel adjustments, careful measurements of the surface 
deviation from the desired parabola were made on the 26- 
m antenna. Measurements were made at three different 
elevation angles: near zenith, near horizon, and at 33 
degrees. The subreflector was removed and returned to 
JPL for precision measurements, and installation of a 
removable vertex plate. While the feedcone was off the 
antenna, antenna waveguide run measurements were 
made in preparation for installation of the modified 
feedcone described above, designated SVU (S-band, 
Venus, Unattended). 

Replacement of the oil seal on one elevation ball screw 
was performed, and an examination of the gearbox was 
made to ensure that no hidden damage existed. The oil 
seal replacement solved the oil leakage problem, and, 
although brass shavings and a piece of broken gear were 
found in the gearcase, no immediate problem is foreseen. 

After reinstallation of the subreflector, new vertex plate, 
and feedcone, pointing evaluation tracks were performed, 
using an automated program called SCOUR (SCan and 
COrrect Using Receiver), which utilizes the SDS-930 to 
provide antenna control. Using Virgo A, 3C144, and 
3C273 as sources, 21-1/2 hours of pointing error 
evaluation were performed. 

Several problems were experienced with the S-band 
maser refrigerator with which the 26-m antenna is 
equipped. During this period, the displacer crosshead 
drive has been overhauled and replaced twice. The system 
is working correctly at this time, although some audible 
roughness exists in the drive mechanism. 

VII. X Band Radar, 8495 MiSz, 400 W 

During a series of radar observations of the comet 
D’ Arrest, the VA949J klystrons used as final amplifiers 
failed and were removed. In preparation for restoring the 
system to service, the traveling-wave tube amplifiers 
(TWTA) used as klystron drivers were tested, and the 
TWT power supplies (Logimeuics Corp. Model A300) 
were modified, with a new HV cable connector for 
improved reliability. This modification required fabrica- 
tion of a complete new back panel assembly and 
rearrangement of connectors. After modification, one of 
the three power supplies developed an intermittent fault 
and was repaired. At this time, the TWTAs are 
operational. 

VIII. Energy Conservation 

The Operations and Data Processing (G-51) and 
Laboratory and Office (G-60) buildings at the Venus 
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Station have each been equipped with a time clock with 
which the lighting {interior and exterior) and air 
conditioning are operated. The air conditioning is on a 
multiday time clock, while the office and laboratory 
lighting is controlled by personnel-set timers, with a total 
time of 10 hours available. Excluding the transmitting 
systems, the air conditioning systems at Venus are the 
largest users of electricity, and automatic turnoff of this 
usage when the station is not manned will reduce 
consumption. 

IX. Deep Space Network High-Power 
Transmitter Maintenance Facility (DSN 
HPTMF) 

The X-3070 klystron loaned to Arecibo Radio Observa- 
tory has been returned in apparent good condition as 
evidenced by a vacuum check. Arecibo did not find it 
necessary to use this klystron. 

At the request of the Transmitter Cognizant Operating 
Engineer (COE), a klystron socket tank, klystron focusing 
magnet, magnet adapter and 500-kW RF water load were 
loaned to Varian Corporation to be used in testing a 
repaired DSN 100-kW klystron. These items were 
returned after testing was completed. 

In continuing support of the DSN mission, four each, 
20-kW klystrons, yarian Model 5K70SG, were tested to 
assure usability in case of need. (These klystrons had been 
obtained from the STDN.) Two of the four klystrons meet 
ail applicable specifications. Three arc detectors, fabri- 
cated by a vendor, were checked, wiring errors corrected, 
repairs made where necessary, and operation and response 
time verified. Additionally, a light source and power 
supply for the crowbar on the HV power supply at DSS 
14 were repaired. 

Routine scheduled maintenance was performed on the 
20-kW test transmitters, and the DSS 13 transmitter 
control system wa» converted to the same solid state 
configuration used in the DSN. During this maintenance 
period, a documentation update program was started to 
ease the planned configc -ation of the system for unat- 
tended operation for uplink transmissions to spacecraft, 
During this maintenance period, the X-3070 m?gnet used 
in the Advanced Systems Demonstration transmitter failed 
due to inadequate coolant flow, and the system is 
currently inoperative. 

As part of a program to improve the spectral purity of 
the uplink transmissions, a high-resolution spectrum 
analysis was made of the carrier emitted by the 20-k\V 


test transmitters and also of the HV power supply used 
with these transmitters. The same spectral components, at 
similar relative amplitudes, were found in both the RF and 
HV dc spectra. Detailed analysis is being performed. 


X. Planetary Radio Astronomy 

In support of the Radio Science experiment “Planetary 
Radio Astronomy” (OSS 196-41-73-01), DSS 13 measures 
and records the radiation received at 2295 MHz from the 
planet Jupiter and various radio calibration sources. These 
measurements use the 26-m antenna, the S-band station 
receiving system, and the NAR. During this period, 109- 
3/4 hours of observations were made, during which the 
radiation received from Jupiter and the calibration sources 
tabulated in Table 2 were measured and recorded. 


XI. Interstellar Microwave Spectroscopy 

In support of the Radio Science E^eriment “Interstel- 
lar Microwave Spectroscopy” (OSS 188-41-52-12) DSS 13 
provided 23 hours of station support, including 17-3/4 
hours of tracking reception at 2273.5 and 2321.729 MHz. 
Sources used for this experiment included 3C58, NGC 
7822, and 118 4*48. 


XII. VLBI Support 

In support of the several programs that utilize Very 
Long Baseline Interferom "try (VLBI) techniques, DSS 13 
devoted 6-3/4 hours of tracking to check out new 
equipment utilized for this purpose. Sources observed 
were 3C84 and 3C273. 


XIII. Clock Synchronization System 

The interval between transmission of clock synchroniza- 
tion signals was shortened at the request of the Viking 
Project. This increased transmission schedule, program- 
med on five-day centers, was instituted in support of the 
encounter, orbit, and landing phases of Vikings I and II 
spacecraft with Mars. 

Although Spain, DSS 63, had some difficulty with their 
receiving system, 18 transmissions were made to them as 
scheduled by DSN Scheduling. Another 21 transmissions 
were made to Australia DSS 43, for a total transmission 
time of 40-3/4 hours. Both stations are now reporting 
good reception and correlation with their local station 
clocks. 
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Table 1. Pulsars observed at DSS 13, June 14 through 
October 10, 1976 


0031-07 

0833-45 

1818-04 

0329+54 

1133 -MO 

1911-04 

0355 + 54 

1237 + 25 

1929+10 

0525 + 21 

1604-00 

1933+16 

0628-28 

1642-03 

2021+51 

0736-40 

1706-16 

2045-16 

082.3 + 26 

1749-28 

2L11+46 



2218 + 47 


Table 2. Radio calibration sources observed at DSS 13, 
June 14 through October 10, 1976 


3C17 

3C273 

3C418 

3C48 

3C274 

NRAO530 

3C84 

3C279 

OJ287 

3C123 

3C286 

PKS0237-23 

3C138 

3C309.1 

FKS2 134-00 

3045 

3C345 

Virgo A 

3C147 

3C348 


3C218 

3C353 
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